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CHAPTER  1 
INTRODUCTION 


1.1  Electronics  in  Aviation 

Since  the  1930’s  aircraft  of  all  varieties  began  adopting  systems  that  required 
electrical  power  for  their  operation  [1],  Over  the  decades,  this  trend  has  only  increased. 
From  avionics  and  flight  control  to  environmental  control,  electronic  systems  have  grown 
and  evolved  into  more  important  and  mission  critical  roles.  As  modern  aircraft  become 
more  reliant  on  electronic  systems,  a  greater  burden  is  being  placed  on  the  electric  power 
generation  to  operate  these  systems.  The  power  demands  are  being  met  with  larger  and 
more  powerful  electrical  generation  systems,  further  enabling  the  capability  of  on  board 
electronics.  However,  with  higher  power  demands,  the  thermal  load  created  by  the  sum 
of  these  systems  has  been  increasing  proportionately. 

Figure  1.1  [1]  and  Figure  1.2  [2]  show  the  increasing  electrical  demands  and 
subsequent  rise  in  the  aircraft  heat  load.  Traditionally  the  primary  sink  for  thermal  energy 
on  an  aircraft  is  the  on  board  fuel.  However,  only  a  finite  amount  of  energy  can  be  stored 
this  way  and  modern  aircraft  are  quickly  approaching  this  limit.  Additionally,  the 
incorporation  of  structural  composite  materials  of  low  thennal  conductivity  and 
emissivity  has  created  the  adverse  secondary  effect  of  insulating  the  aircraft.  Altogether, 
these  effects  have  forced  a  greater  burden  on  thermal  management  systems  (TMS). 
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Figure  1.1  Aircraft  Electric  Power  Demands  over  Time  [1] 
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Figure  1.2  Aircraft  Avionics  Heat  Load  Progression  over  Time  [2] 


The  majority  of  the  TMS  developed  have  mainly  focused  on  moving  and 
distributing  the  heat  loads  away  from  critical  components.  However,  this  approach  is 
limited.  Modem  aircraft  electronics  are  also  continually  becoming  lighter  and  more 
compact,  thus  increasing  the  power  density.  This  means  that  thennal  management  will 
need  to  be  more  efficient  in  removing  large  transient  heat  loads  and  need  to  be  collocated 
with  the  heat  source.  “5th  generation  aircraft  have  3-5  times  the  heat  load  of  legacy 
platfonns  while  being  limited  in  their  ability  to  reject  heat  to  the  environment.”  [3] 


1.2  Research  Concentration 

A  relatively  new  challenge  for  thennal  management  that  is  emerging  has  to  do  with 
the  presence  of  electromechanical  actuation  systems  (EMAS)  being  developed  for  aircraft 
primary  flight  control  surfaces.  The  emergence  of  EMAS  has  brought  about  measurable 
advantages  in  weight  savings,  maintainability,  complexity,  and  efficiency.  However,  they 
are  not  without  drawbacks.  Prior  solutions  for  flight  control  actuation  used  hydraulic 
actuators.  One  of  the  advantages  these  systems  had  over  current  proposed  EMAS 
replacement  was  the  presence  of  hydraulic  fluid  acting  secondarily  as  a  coolant.  The 
removal  of  hydraulic  actuation  systems  in  favor  of  EMAS  means  the  removal  of  its  very 
effective  thennal  management  process,  thus  creating  a  new  challenge  [4]  [5]  [6]. 

Some  specific  solutions  considered  for  the  thennal  management  of  EMAS 
incorporate  the  following: 

•  Thennal  energy  storage  (on  aircraft) 

o  Phase  Change  Material  (PCM)  for  transient  heat  absorption 
o  Heat  pipe,  loop  heat  pipe,  or  other  high  conductive  materials  to  conduct 
the  heat  from  PCM  to  airframe  or  fuel 

•  Thermal  energy  removal  (to  ambient) 

o  Convection,  air  jet  cooling 
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o  Convection,  liquid  to  air  heat  exchanger 
o  Convection,  air  cooling  with  fan 


Each  of  these  methods  has  advantages  and  disadvantages,  whether  it  is 
perfonnance,  integration,  or  cost. 

In  this  study,  the  suitability  of  forced  convection  with  air  will  be  examined. 
Forced  convection  is  a  very  common  and  effective  method  used  for  thermal  management, 
particularly  in  cooling  electronic  equipment.  There  exist  multiple  approaches  to 
accomplish  forced  convection  cooling  using  air.  One  common  approach  is  to  use  an 
enhanced  surface  or  heat  exchanger  to  remove  the  thennal  energy  away  from  critical 
components  and  convect  it  to  an  airstream.  The  most  common  source  of  the  air  flow 
would  be  an  electric  fan  due  to  its  low  cost,  light  weight,  small  size,  and  low  energy 
consumption.  However,  on  an  aircraft  there  exist  other  sources,  such  as:  the  aircraft 
environmental  control  system  (ECS)  air,  engine  bypass  air,  or  channeling  external  air 
passing  over  the  aircraft  surface  (ram  air).  Although,  all  methods  have  considerable 
systematic  differences,  they  all  involve  transmission  of  thermal  energy  to  the  atmosphere 
surrounding  the  aircraft.  Therefore,  an  investigation  into  properties  governing  this 
process  will  be  applicable  to  all  similar  solutions. 

For  the  case  of  the  electromechanical  actuation  systems  for  flight  control,  the 
electromechanical  actuator  (EMA)  will  be  housed  in  what  are  known  as  wing  bays, 
unpressurized  enclosures  that  are  isolated  from  the  ECS  air  flow  but  exposed  to  changes 
in  the  atmospheric  conditions.  In  other  words,  for  most  aircraft  the  internal  conditions  of 
components  within  a  wing  bay  will  be  static  air  resembling  the  environmental  conditions 
of  the  atmosphere  for  a  given  altitude.  Physical  characteristics  of  the  atmosphere  vary 
relative  to  altitude  (e.g.,  temperature,  pressure,  humidity,  vapor  pressure,  speed  of  sound, 
concentration  of  gases,  etc.).  These  properties  likely  have  a  drastic  effect  on  the 
perfonnance  of  forced  convective  heat  transfer  when  trying  to  dissipate  thermal  energy  to 
the  atmosphere. 

Figure  1.3  [7]  [8]  shows  the  physical  characteristics  of  atmospheric  air  for 
altitudes  from  sea  level  to  100  km.  The  right  side  of  the  chart  shows  the  different  layers 
of  the  atmosphere  as  well  as  representative  altitudes  of  various  man-made  aircraft  and 
natural  occurring  phenomena.  Figure  1.3  also  shows  that  a  typical  flight  ceiling  for 
commercial  and  military  air  breathing  aircraft  is  about  40,000  feet  (12km)  above  mean 
sea  level  (ASL). 
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Figure  1.3  Air  Pressure,  Density,  Temperature  and  Speed  of  Sound  Relative  to  Altitude  [7] 

18] 


Of  the  four  properties  shown  in  Figure  1.3  (pressure,  density,  temperature,  and 
speed  of  sound),  three  are  interconnected  through  the  ideal  gas  law,  ( PV  —  nRT )  and  can 
significantly  influence  heat  transfer.  Changes  to  ambient  temperature  will  mostly  affect 
the  temperature  gradient  between  the  sources  to  the  sink.  The  changes  to  this  property 
can  be  studied  and  characterized  through  experimentation  at  sea  level.  However,  this  is 
not  the  case  for  changes  to  the  pressure  and  density.  Although,  changes  to  density  and 
pressure  are  well  understood  analytically,  there  lacks  experimental  data  to  support  this 
understanding  as  it  applies  to  predictive  methods  of  cooling  for  EMAS.  From  sea  level  to 
40,000  feet  the  pressure  and  density  of  the  air  decrease  by  roughly  80%.  The  effects  this 
can  have  on  the  performance  of  convective  heat  transfer  could  be  profound. 

When  designing  forced  convective  heat  transfer  solutions  for  all  varieties  of 
systems  engineers  use  several  methods  for  predicting  the  performance  of  the  designs. 
Through  observation,  theoretical  analysis,  and  experimentation,  certain  relationships  have 
been  developed  describing  the  processes  of  convective  heat  transfer.  These  relationships 
have  been  quantified  as  correlations.  These  correlations  can  be  used  to  describe  such 
characteristics  as  the  heat  transfer  coefficient  (HTC)  and  friction  factor,  which  are  both 
important  theoretical  properties  when  defining  the  processes  of  heat  and  mass  transfer. 
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Correlations  that  would  be  relevant  to  heat  transfer  perfonnance  for  electromechanical 
actuator  (EMA)  TMS  have  been  developed,  supported  in  theory,  and  observed 
experimentally.  However,  they  have  not  been  extensively  vetted  in  atmospheric 
conditions  experienced  by  commercial  and  tactical  military  aircraft. 


1.3  Purpose 

The  purpose  of  this  thesis  was  to  experimentally  find  the  forced  convection  heat 
transfer  coefficient  and  pressure  drop  relationships  in  a  low  pressure  environment  under 
varying  conditions  of  pressure  and  velocity,  to  study  methods  of  heat  transfer 
enhancement,  and  to  validate  experimentally  existing  heat  transfer  and  pressure  drop 
correlations  relevant  to  EMAS  cooling  at  low  pressure. 
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CHAPTER  2 
BACKGROUND 


2.1  Literature  Review 

In  order  to  begin  planning  the  experimental  investigation  into  the  topic  of  forced 
convective  heat  transfer  at  high  altitude,  a  comprehensive  understanding  of  relevant 
publications  must  first  be  established.  There  has  been  a  gap  in  published  experimental 
evidence  regarding  high  altitude  (40,000  feet  above  mean  sea  level  (ASL)),  low  Reynolds 
number  convective  heat  transfer,  relevant  to  the  thennal  management  of  electronic 
systems.  However,  that  is  not  to  say  that  there  is  no  supporting  literature  related  to  this 
topic  at  altitudes  greater  than  sea  level. 

The  literature  review  on  the  works  in  this  field  begins  in  the  late  1980’s  from  an 
investigation  into  heat  transfer  process  of  electronics  at  high  altitude.  In  1989,  Bar- 
Shalom  published  a  thesis  titled  Altitude  Effects  on  Heat  Transfer  Processes  in  Aircraft 
Electronic  Equipment  Cooling  [2],  This  work  took  an  analytical  approach  to 
characterizing  the  changes  in  environment  conditions  and  their  impact  on  cooling 
electronics  in  wing  bays.  It  also  took  steps  to  detail  the  contributing  factors  and  issues  for 
heat  transfer  processes  throughout  a  range  of  altitude  from  0  -  70,000  feet  ASL  and 
offered  insight  on  how  to  optimize  cooling  for  these  conditions.  One  conclusion  of  this 
study  was  when  considering  only  pressure  changes,  the  heat  transfer  thermal  resistance 
increases  exponentially  with  increasing  altitude.  This  is  due  to  a  reduced  heat  transfer 
coefficient  (HTC)  at  low  pressures.  With  a  lower  HTC,  the  temperature  difference 
between  the  source  and  fluid  becomes  larger  as  a  result  of  the  higher  thermal  resistance. 
Figure  2.1  [2]  shows  the  relationship  that  the  decreasing  convective  HTC  has  with  the 
temperature  rise  between  the  electronics  and  ambient,  with  changing  altitude. 


Figure  2.1  Temperature  Rise  per  Heat  Flux  vs.  Convective  HTC  |2]. 
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With  a  lower  HTC  the  electronics  will  have  a  higher  operating  temperature  since 
it  requires  a  higher  temperature  difference  to  remove  the  same  amount  of  heat.  This 
drives  system  requirements  for  electronics  cooling,  since  having  a  higher  temperature 
difference  means  that  the  electronics  are  operating  closer  to  unsafe  temperatures.  This 
becomes  a  key  element  in  aircraft  electronic  TMSs,  which  can  minimize  the  working 
temperature  of  components  at  altitude  by  designing  systems  to  have  the  largest  possible 
HTC. 

Another  important  point  from  Bar-Shalom’s  work  was  regarding  the  contribution 
of  radiation  to  the  thennal  management  process  at  high  altitude.  Before  this  study  the 
practice  was  to  ignore  radiation  heat  transfer.  However,  with  the  diminished  HTC  at 
altitude,  the  temperature  of  the  heat  source  or  the  electronics  increases  and  the  effect  of 
radiation  becomes  a  larger  percentage  of  the  overall  heat  transfer.  Further,  this  thesis 
suggested  ways  to  design  optimal  cooling  methods,  but  it  falls  short  with  regard  to 
experimental  evidence  to  support  its  conclusions. 

Belady  [9]  sought  to  provide  an  overview  on  the  impact  altitude  has  on  the  fan 
cooling  of  electronics.  The  majority  of  thennal  management  solutions  for  electronics 
favors  the  use  of  electronic  fans  to  provide  fluid  motion.  This  is  because  they  have  low 
cost  and  weight,  reasonable  reliability,  and  are  relatively  efficient.  Perfonnance  of  the  fan 
is  also  affected  by  changes  in  air  pressure  and  density.  Therefore,  when  studying  fan 
cooling  above  mean  sea  level  (ASL)  this  factor  must  be  accounted  for.  Belady  found  that 
at  10,000  feet  ASL  a  fan  rotating  at  the  same  speed  produces  a  decreased  dynamic  head 
that  is  proportional  to  the  decrease  in  air  density,  while  the  power  consumed  by  the  fan  is 
found  to  be  1/3  less.  In  other  words,  the  ability  of  a  constant  volume  axial  fan  to  drive  a 
pressure  head  decreases  due  to  the  air  having  a  lower  density.  This  exacerbate  the  effect 
on  thennal  management  due  to  the  influence  density  has  on  the  heat  transfer  coefficient 
in  both  laminar  and  turbulent  flow  conditions.  Another  observation  made  by  Belady  is 
that  from  sea  level  to  10,000  feet  ASL  the  temperature  rise  in  the  air  exiting  the  heat  sink 
increased  by  50%.  This  is  an  indication  of  the  heat  sink  surface  temperature  rising,  as 
there  are  fewer  air  molecules  to  remove  the  heat,  resulting  in  a  higher  average 
temperature  of  the  air.  This  is  a  direct  result  of  the  negative  impact  low  density  air  has  on 
the  heat  transfer  coefficient.  Belady  uses  the  following  basic  and  well  established 
correlations  to  obtain  results.  Equation  2-1  [9]  and  Equation  2-2  [9]  are  the  Nusselt 
number  correlations  Belady  used  for  both  laminar  and  turbulent  flow  regimes. 

Equation  2-1  [9] 

NM-Laminar  ~  Chi’  <  2,000 

Equation  2-2  [9] 

NuTurbulent  =  0.023 Re°-8Pr0A)  Re  >  10,000 


where  Chl  is  a  constant. 

In  summary  Belady  found: 
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•  “For  laminar  flow,  pressure  drop  is  independent  of  altitude  and  is  a  linear  function 
of  volumetric  flow.  Thus,  the  system  impedance  curve  is  the  same  for  all 
altitudes”  [9]. 

•  “For  laminar  flow,  the  heat  transfer  coefficient  is  independent  of  altitude  and 
volumetric  flow.  It  is  constant  for  a  given  geometry  as  long  as  the  flow  is  laminar. 
This  implies  that  the  temperature  rise  to  the  local  air  temperature  is  also  constant” 
[9]. 

•  “For  turbulent  flow,  pressure  drop  is  a  function  of  altitude  and  the  square  of  the 
volumetric  flow.  Thus,  the  system  impedance  curve  decreases  by  the  density  ratio 
as  the  altitude  increases”  [9], 

•  “For  turbulent  flow,  the  heat  transfer  coefficient  is  a  function  of  density  and 
volumetric  flow  both  to  the  0.8  power  for  circular  ducts.  This  implies  that  the 
local  air  temperature  rise  is  inversely  proportional  to  density  and  volumetric  flow 
both  to  the  0.8  power”  [9]. 


Belady  also  produced  a  table  (See  Figure  2.2)  expressing  the  behavioral  changes  of 
the  axial  fan  as  well  as  the  heat  transfer  over  a  heat  sink  in  relation  to  changes  in  altitude, 
or  density. 


Description 

Parameter 

Relation 

Fan  Performance 

Fan  Air  Flow 

.  AP  cc  p 

Fan  Power 

HP  cc  p 

Heat  Capacity 

Temperature  Rise 

AT  oc  l/(p  G) 

Laminar  Flow 

Pressure  Drop 

AP  oc  G 

Heat  Transfer  Coeff 

h  =  constant 

Temperature  Delta 

AT  -  constant 

Turbulent  Flow 

Pressure  Drop 

AP  oc  p  G 2 

Heat  Transfer  Coeff 

h<xp°-8G0-8 

Temperature  Delta 

ATcc  Up08  G0-8 

Figure  2.2  Summary  of  Results  from  Belady  [9] 


Belady’s  investigation  supports  the  applicability  of  current  methods  used  to 
predict  the  heat  transfer  coefficient  and  pressure  drop  at  altitudes  from  sea  level  to  10,000 
feet  ASL.  However,  there  is  a  significant  difference  between  10,000  and  40,000  feet 
ASL.  The  ambient  air  used  in  Belady’s  experiment  never  fell  below  67%  of  atmosphere 
while  modem  aircraft  operate  closer  to  20%  atmosphere  and  in  some  cases  lower.  This 
margin  is  too  great  to  ignore. 

In  1999,  Wong  [10]  [11]  conducted  an  experiment  on  the  effects  of  altitude  on 
electronics  cooling  perfonnance.  The  difference  from  Belady’s  study  is  that  Wong 
managed  to  gather  experimental  data  by  simulating  altitudes  from  0  -  16,400  feet  ASL  (5 
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km).  At  16,400  feet  ASL  the  atmospheric  pressure  is  53%  of  that  at  sea  level.  This  is 
closer  to  the  20%  conditions,  or  40,000  feet  MLS,  but  still  falls  short.  Wong’s 
conclusions  were  meant  to  aid  in  the  TMS  design  for  specific  electronics  operating  at 
high  altitudes.  However,  there  were  still  some  good  takeaways  from  this  study  that  apply 
to  the  use  of  analytical  methods  for  predicting  heat  transfer.  The  primary  conclusions 
gave  support  for  two  methods  of  predicting  heat  transfer  perfonnance,  the  superposition 
principle  for  predicting  the  surface  temperature  of  a  flat  plate  heater  array  and  the  thennal 
wake  function  of  a  flat  plate.  This  approach  was  better  suited  for  the  prediction  of  surface 
temperatures  based  on  a  reference  temperature,  rather  than  explicit  prediction  of  the 
HTC.  Although,  it  did  offer  some  insight  into  how  heat  transfer  is  influenced  by  altitude. 
For  instance,  Wong  confirmed  that  the  steady  state  surface  temperatures  at  higher 
altitudes  have  a  much  higher  value  than  they  would  at  sea  level,  given  constant  velocity 
air  flows  and  that  this  effect  is  more  prevalent  at  higher  airspeeds.  Overall,  this 
experiment  falls  short  in  completing  a  direct  comparison  of  the  changes  in  the  heat 
transfer  coefficient  due  to  altitude. 

Xu  [12]  explored  ability  of  using  correction  factors  to  predict  forced  convective 
heat  transfer  at  high  altitude,  the  non-dimensional  pressure  drop  and  thermal  resistance. 
Xu  created  correction  factors  derived  from  non-dimensional  heat  transfer  and  fluid  flow 
parameters  (i.e.  Reynolds  number)  to  predict  the  change  in  perfonnance  of  a  complex 
processor  cooling  geometry  at  an  altitude  of  3,000  meters  (9,800  feet)  ASL.  The 
advantage  of  this  approach  was  that  it  produced  a  more  general  use  method  for 
approximating  the  HTC  for  any  flat  plate  geometry  and  flow  regime.  The  results  of  this 
experiment  support  the  use  of  analytical  methods  for  predicting  heat  transfer  and  pressure 
drop,  to  0.7  atm  (9,800  feet  ASL). 

One  industrial  report  done  by  EPAC  Software  Inc.  [13]  explored  fan  cooling  at 
altitudes  higher  than  16,400  feet  ASL.  They  sought  to  improve  their  model  for  predicting 
the  perfonnance  of  an  aircraft  ECU  at  altitudes  up  to  80,000  feet  ASL.  The  investigation 
was  launched  when  EPAC  found  a  large  discrepancy  in  the  predicted  temperatures  of  the 
ECU  coldplate  at  altitudes  above  30,000  feet  ASL.  The  method  of  predicting  the 
temperature  was  based  on  a  proprietary  software  setup  for  their  specific  application. 
However,  the  study  revealed  the  problem  to  be  primarily  with  the  perfonnance  of  the  fan. 
It  was  found  that  the  method  for  predicting  dynamic  pressure  head  of  a  fan  cannot  be 
used  above  30,000  feet  ASL  due  to  a  disconnection  that  happens  where  constant 
volumetric  flow  rate  (VFR)  is  no  longer  proportional  to  rotational  speed.  Instead,  a  mass 
flow  rate  approach  must  be  used.  This  study  provides  valuable  insight  into  the  field  but 
also  serves  to  raise  more  questions  about  the  validity  of  some  very  standard  methods  of 
predicting  fan  cooling  in  high  altitude  conditions,  such  as  the  validity  of  fan  impedance 
curves. 


2.2  Synopsis 

The  results  of  the  literature  review  uncovered  a  gap  in  experimental  evidence 
supporting  analytical  methods  and  empirical  correlations  of  predicting  convective  fan 
cooling  performance  relevant  to  the  thennal  management  of  EMAS,  at  altitudes  above 
16,000  feet  ASL.  Further,  the  review  also  conveyed  the  importance  of  having  an  accurate 
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and  validated  theoretical  approach  to  predicting  heat  transfer  in  the  low  pressure 
environment.  At  altitudes  up  to  and  above  40,000  feet  ASL  the  lower  density  air  has  a 
significant  impact  on  the  perfonnance  of  both  fan  perfonnance  and  the  convective  heat 
transfer  cooling.  This,  in  addition  to  the  industry  thrust  to  adopt  more  high  power 
electronics  for  mission  critical  systems,  drives  the  need  to  conduct  further  research  into 
validating  current  methods  of  heat  transfer  prediction  at  high  altitude.  The  implication  of 
providing  some  insight  into  this  research  area  would  be  to  further  advance  the 
understanding  of  high  altitude  convective  cooling  systems.  This  could  be  in  the  form  of 
validating  the  applicability  of  existing  methods  at  higher  altitude,  offering  a  better 
analytical  approach,  such  as  a  new  or  enhanced  correlation,  or  providing  insight  for  better 
optimized  system  elements,  such  as  heat  sink  or  fan  design. 
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CHAPTER  3 
EXPERIMENT 


3.1  Design  Concept 

To  study  forced  convective  heat  transfer  at  high  altitude  conditions  an  experiment 
needed  to  be  constructed  that  allows  for  the  control  of  properties  of  convective  heat 
transfer  (air  velocity  and  temperature  difference)  as  well  as  the  environment  in  which  it 
occurs  (atmospheric  pressure).  It  was  determined  that  the  best  means  to  accomplish  this 
was  through  the  construction  of  an  open  loop  wind  tunnel,  with  a  fan  driving  the  flow  at 
the  beginning  of  the  tunnel  and  a  heater/heat  sink  at  the  end  of  the  tunnel  before  the  exit. 
A  high  speed  fan  was  used  to  drive  the  airflow  within  the  wind  tunnel  and  allowing  for  a 
wide  range  of  Reynolds  number  within  the  test  section.  Airflow  was  measured  using  a 
differential  pressure  sensor  and  Pitot  tube  within  the  duct  channel.  Lastly  the  wind  tunnel 
was  placed  within  an  environmental  chamber  to  simulate  the  conditions  of  high  altitude. 
Figure  3.1  shows  a  conceptual  diagram  of  the  layout  for  the  open  loop  wind  tunnel. 


Figure  3.1  Conceptual  Design  of  Experimental  Test  Setup 


To  give  greater  control  over  the  flow,  a  ball  valve  was  placed  at  the  entrance  of  the 
wind  tunnel  to  provide  a  means  to  restrict  flow  thus  allowing  for  control  of  the  velocity. 
The  heat  sink  was  attached  to  a  heater  that  would  provide  thermal  energy  that  was 
transferred  to  the  fluid.  Measurements  for  the  accurate  characterization  of  this  process 
such  as  temperature,  absolute,  dynamic  and  static  pressures  and  power  dissipation  were 
taken. 

Since  the  experiment  sought  to  study  forced  convection  at  high  altitude,  a  means 
was  obtained  to  lower  the  ambient  pressure  in  which  this  open  loop  wind  tunnel  operates. 
The  primary  atmospheric  difference  concerning  the  forced  convective  heat  transfer 
coefficient  from  sea  level  to  50,000  feet  ASL  is  the  density  of  air.  Temperature  is  also  a 
key  difference  within  this  altitude  range,  but  it  does  not  affect  the  FITC.  Therefore,  the 
objective  would  be  to  obtain  a  means  to  lower  the  air  density  through  the  control  of 
ambient  pressure.  This  was  achieved  by  obtaining  an  environmental  chamber.  The 
chamber  contains  an  airtight  space  that  utilizes  a  vacuum  pump  to  decrease  the  internal 
pressure,  simulating  atmospheric  conditions  at  altitudes  above  sea  level. 
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3.2  Environmental  Chamber 


For  the  purpose  of  simulating  atmospheric  conditions  above  sea  level,  a  Tenney 
Versa  3  environmental  chamber  was  obtained.  The  chamber,  originally  built  in  the 
1950’s,  was  designed  to  control  pressure  and  temperature  for  a  rectangular  volume  of  21” 
x  21”  and  25”  in  height.  Figure  3.2  shows  the  outside  of  the  chamber  as  well  as  the  sealed 
workspace  available  within. 


Figure  3.2  Tenney  Versa  3  Environmental  Chamber  (Left)  with  Working  Volume  (Right). 


The  chamber  was  recently  equipped  with  an  updated  control  unit  for  active  control 
of  chamber  pressure  and  temperature.  The  pressure  is  reduced  via  a  rotary  vane  vacuum 
pump  powered  by  an  electronic  motor.  The  pump  itself  can  achieve  a  vacuum  of  about 
0.080  TORR  (211,000  feet  ASL).  However,  due  to  imperfect  sealing  of  the  chamber  the 
maximum  vacuum  that  was  measured  to  be  about  100  TORR  and  later  80  TORR  (52,000 
feet  ASL)  was  achieved  through  the  addition  of  a  secondary  vacuum  pump.  The  heating 
and  cooling  of  the  air  within  the  chamber  is  done  with  an  electric  coil  heater  and  vapor- 
compression  cycle  refrigeration  system.  Although  these  systems  were  available, 
temperature  control  was  not  used  for  this  experiment.  The  primary  reason  was  negative 
effect  that  the  refrigeration  system  had  on  the  absolute  pressure  control,  by  causing  (±0.2 
PSI)  fluctuations  during  operation.  Secondly,  there  was  also  little  value  in  varying  the 
chamber’s  temperature  since  the  HTC  measurement  is  independent  of  the  temperature 
gradient  between  the  heat  sink  and  ambient  air  temperature. 


3.3  Wind  Tunnel  Hardware 

The  four  major  components  responsible  for  the  execution  of  the  experiment  are  the 
fan,  valve,  heater/heat  sink,  and  data  acquisition  (DAQ)  System.  The  design  and  use  of 
all  four  will  be  discussed  in  detail  within  this  section. 
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3.3.1  Fan 


The  Fan  that  was  chosen  is  an  AMETEK  Rotron  Propimax  2L  3715SF.  This  is  a 
very  high  perfonnance  axial  fan  about  3.0  inches  in  diameter  capable  of  spinning  in 
excess  of  25,000  RPM.  This  fan  proves  to  be  a  good  fit  for  the  experiment  because  of  its 
compact  size  and  high  volumetric  flow  rate.  The  fan  is  powered  by  a  DC  brushless  motor 
that  rotates  a  2-blade  cast  aluminum  rotor.  The  fan  requires  a  DC  voltage  between  10-28 
VDC  to  operate  and  will  draw  1.3  Amps  (at  28VDC)  while  driving  air  at  1.202  kg/m  . 
However,  this  power  draw  will  decrease  as  the  density  of  air  being  driven  decreases  [14]. 
The  fan  will  not  run  at  voltages  lower  than  -10VDC  as  this  appears  to  be  a  minimum 
threshold  voltage  needed  to  drive  the  motor  and  electronics.  This  also  means  that  there  is 
a  minimum  flow  rate  that  the  fan  can  reach  when  airflow  is  being  controlled  with  fan 
voltage.  This  is  the  reason  a  control  valve  was  also  for  controlling  the  airflow.  For  more 
information  on  the  Propimax  2L  fan,  a  data  sheet  is  available  in  APPENDIX  A. 

A  fan  curve  is  a  metric  for  evaluating  the  performance  of  any  given  fan  against 
various  loads.  Figure  3.3  shows  the  fan  curve  for  the  Propimax  2L  3715SF  as  provided 
by  the  manufacturer.  This  shows  for  a  given  voltage  input  the  fan  will  drive  a  certain 
volume  of  air  of  0-60  liters  per  second  against  a  static  pressure  load  of  0-4  inches  of 
water  (inH20).  Therefore,  to  determine  the  airflow  that  this  fan  will  produce  in  a  given 
system  the  pressure  drop  of  this  system  must  be  determined. 


The  fan  curve  in  Figure  3.3  shows  the  perfonnance  of  the  fan  at  an  air  density  of 
1.202  kg/nr.  Since  this  experiment  would  operate  at  reduced  densities  of  air,  the  amount 
of  pressure  the  fan  can  drive  would  also  be  reduced.  Therefore,  in  order  to  accurately 
characterize  the  fans  capability  at  low  densities,  fan  curves  were  needed  for  its  operation 
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at  reduced  pressure.  Work  done  by  Lin  [14],  at  the  University  of  Central  Florida,  directly 
contributed  to  characterizing  this  fan  at  reduced  pressures.  Lin  gathered  experimental 
data  on  the  performance  of  the  Propimax  2L  fan  and  constructed  fan  curves  for  pressures 
at  1.0  atm,  0.7  atm,  0.5  atm  and  0.2  atm.  The  plots  from  this  work  can  be  seen  in  Section 
3.4.4. 


3.3.2  Valve 

In  this  experiment  a  ball  valve  is  used  to  control  the  airflow  entering  the  open 
loop  wind  tunnel.  The  valve  works  as  a  sphere  or  ball  with  a  cylindrical  hole  through  it 
rotates.  The  cylindrical  hole  allows  for  the  passage  of  fluid  from  one  side  of  the  valve  to 
the  other.  As  the  valve  rotates  from  0-90°  this  cylindrical  hole  rotates  into  the  side  of  the 
valve  partially  covering  it  and  effectively  changing  the  size  of  the  hole  that  fluid  can  now 
pass  from  one  side  of  the  valve  to  the  other.  By  decreasing  this  orifice  size  the  impedance 
to  the  system  is  increased,  thus  decreasing  the  volumetric  flow  rate  (VFR).  This  allows 
for  active  modulation  of  the  airspeed  with  the  valve. 

Another  requirement  for  this  system  is  that  it  must  be  remotely  operated.  This  is 
because  the  regulation  of  airflow  will  need  to  be  done  actively  and  the  experiment  will  be 
completely  contained  inside  the  environmental  chamber.  Therefore,  an  electronically 
controlled  valve  must  be  used.  The  Plast-O-Matic  2-Way  electronically  actuated  ball 
valve  was  selected  for  this  purpose.  There  are  several  types  of  ball  valve  available  for  this 
product  line.  A  ball  valve  was  chosen  that  fit  a  3.0  inch  diameter  pipe  coupling  and  had 
1 .0  inch  diameter  inner  orifice  size.  The  actuator  that  controlled  the  valve  was  powered 
by  19  VDC  and  responded  to  an  analog  signal  of  0-10  VDC  for  position  control  of  0° 
(closed)  to  90°  (fully  open). 


3.3.3  Heater/Heat  Sink 

The  heater/heat  sink  would  provide  the  source  for  the  thennal  energy  generation 
as  well  as  the  geometry  dissipating  the  thermal  energy  to  the  fluid  stream.  The  purpose  of 
this  experiment  required  that  the  geometry  used  as  the  heat  sink  should  be  able  to 
characterized  using  widely  accepted  methods  of  correlation.  Additionally,  it  should  be  a 
variation  of  an  existing  geometry  type  that  is  widely  adopted  within  industry,  so  as  to 
keep  its  applicability  to  what  might  be  used  to  cool  EMAS.  Therefore,  a  parallel  plate 
sink  was  chosen.  This  type  of  heat  sink  is  very  popular  in  forced  air  cooling  of  electronics 
because  of  its  good  perfonnance  and  low  cost.  The  other  constraint  on  the  geometry  is  the 
size.  Larger  geometries  will  limit  the  maximum  airflow  velocity  over  the  heat  sink,  thus 
limiting  the  range  of  Reynolds  number,  as  well  as  making  the  integration  into  the 
experiment  more  cumbersome  given  the  finite  amount  of  space  within  the  environmental 
chamber.  However,  a  smaller  geometry  is  more  adversely  affected  by  error  due  to 
experimental  measurement  methods,  such  as  thennocouples  imbedded  altering  the 
thennal  conduction  and  obstructions  within  the  air  stream  such  as  thermocouple  (TC) 
probes  and  a  Pitot  tube. 

Figure  3.4  shows  the  heat  sink  that  was  chosen  for  the  experiment.  It  was 
manufactured  by  Heat  Sink  U.S.A.  It  is  a  4.0  inch  long  parallel  plate  extruded  aluminum 
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heat  sink.  The  material  is  6063-T5  unfinished  aluminum  (k  =  209  W/m-K,  s  =  0.04-0.07). 
Its  width  is  1.813  in.  and  height  1.25  in.  with  a  total  of  6  fins.  This  geometry  offered  a 
good  fit  given  the  scale  of  the  experiment. 


Figure  3.4  Heat  Sink  Photograph  (Left)  and  Profile  Dimensions,  in  Inch  (Right) 


A  thin  film  resistive  heater  was  permanently  attached  to  the  base  of  the  heat  sink. 
The  type  chosen  was  an  etched  foil  sandwiched  between  two  layers  of  Polyimide  film. 
The  entire  element  is  0.01  in.  thick.  This  would  allow  for  easy  integration  and  attachment 
to  the  bottom  of  the  heat  sink,  as  seen  in  Figure  3.5.  In  order  to  have  equal  heating  of  the 
heat  sink  over  the  entire  area  of  the  base,  the  heating  element  must  cover  the  entire  1.813 
in.  x  4.0  in.  surface  of  the  base.  This  requires  two  identical  heaters  each  2.0  in.  *  2.0  in. 


Figure  3.5  Thin  Film  Heaters  Installed  on  the  Base  of  Heat  Sink. 
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The  thin  film  heater  was  manufactured  by  Omega,  KMLV-202/10-P.  This  item 
had  a  watt  density  of  10W/in2  at  28  VDC  and  total  power  output  of  40W  each,  yielding  a 
heat  input  range  of  0-80W.  Also,  this  heater  comes  with  pressure  sensitive  adhesive  on 
the  flat  side  of  the  element  for  ease  of  attachment  to  another  flat  surface.  Ultimately  the 
heaters  will  be  wired  in  parallel  so  that  they  both  have  an  equal  amount  of  voltage 
dropped  across  their  resistive  element. 


3.3.4  DAQ 

The  last  critical  piece  of  hardware  needed  for  the  setup  was  the  data  acquisition 
system  (DAQ).  This  item  would  serve  as  the  primary  means  of  control  for  the  experiment 
as  well  as  being  responsible  for  logging  all  of  the  experimental  data.  The  device  chosen 
to  do  this  is  a  National  Instruments  (NI)  Compact  Reconligurablc  I/O  Module  (cRIO). 
The  cRIO-9022  chassis  offered  the  capability  of  both  a  real-time  (RT)  controller  and 
field-programmable  gate  array  (FPGA)  module  to  collect  data  from  up  to  8  unique  NI 
modules  and  transmit  the  information  through  Ethernet  to  a  computer  at  up  to  1000 
Mbits/s.  For  this  experiment  a  total  of  7  cRIO  modules  were  utilized.  The  model  numbers 
for  these  modules,  their  quantities,  the  mode  it  was  operated  in  (RT/FPGA),  and  the 
specific  function  within  the  experiment  are  listed  below. 

•  4  x  NI-92 1 1  Thermocouple  Cards  (RT) 

o  16  x  Input  -  Thennocouple  measurements 

•  lx  NI-9205  Analog  Input  (RT) 

o  4  x  Input  -  Differential  Pressure  Sensors 
o  1  x  Input  -  Absolute  Pressure  Sensor 
o  lx  Input  -  Shunt  Resistor  Voltage 

•  lx  NI-9263  Analog  Output  (RT) 

o  1  x  Output  -  Heater  Voltage 
o  1  x  Output  -  Fan  Voltage 
o  1  x  Output  -  Valve  Position 

o  1  x  Output  -  Power  Supply  Current  Limit  (5V  Constant) 

•  1  x  NI-94 1 1  Digital  I/O  (FPGA) 

o  lx  Input  -  Fan  Encoder  Pulse  Count 

3.4  Wind  Tunnel  Design 

Ideally  the  wind  tunnel  should  have  no  influence  on  the  process  of  forced 
convection  and  should  be  completely  removed  as  a  variable  within  the  experiment.  In 
other  words,  the  wind  tunnel  should  deliver  perfectly  uniform,  constant  velocity  air  over 
the  entire  wetted  area  of  the  heat  sink.  In  practice,  this  is  never  the  case  since  no  system  is 
ideal.  However,  steps  can  be  taken  to  minimize  any  adverse  effects  the  wind  tunnel  has 
on  the  quality  of  the  flow  over  the  heat  sink.  To  do  this  several  guidelines  that  have  been 
developed  and  proven  were  adopted  for  this  experiment.  Many  of  these  guidelines  are 
published  as  standards  in  AMCA:  210-99  [15].  This  reference,  along  with  some  others, 
offers  recommendations  for  duct  settling  lengths  needed  to  create  a  uniform  flow 
necessary  for  the  reliable  and  accurate  dynamic  pressure  measurements  via  a  Pitot  tube. 
Figure  3.6  shows  how  the  minimum  length  requirements  are  used  between  critical 
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components  within  a  duct  to  allow  for  the  accurate  air  velocity  measurement  with  a  Pitot 
tube. 


PL.l  PL.2  PL. 3 


NOTES 

1.  Dotted  lines  on  fan  inlet  indicate  an  inlet  bell  and  one  equivalent  duct  diameter  which  may  be  used  for  inlet  duct 
simulation.  The  duct  friction  shall  not  be  considered. 

2.  Dotted  lines  on  the  outlet  indicate  a  diffuser  cone  which  may  be  used  to  approach  more  nearly  free  delivery. 

Figure  3.6  Design  Length  Strategy  for  Duct  Air  Velocity  Measurement  [15] 


These  length  recommendations  are  based  off  of  multiples  of  the  ducts  diameter  or 
effective  diameter,  known  as  the  hydraulic  diameter  ( DH ).  This  allows  for  the  standards  to 
be  scalable  to  any  shape  of  duct  (e.g.,  rectangular,  square).  The  following  guidelines  [15] 
[16]  for  facilitating  accurate  airflow  measurement  were  used: 

•  Located  the  Pitot  tube  section  providing  5  or  more  duct  diameters 
upstream  and  1-1/2  or  more  diameters  down  stream  of  Pitot  tube  free  of 
elbows,  size  changes  or  obstructions. 

•  Use  a  honeycomb  type  of  flow  straightener  5  duct  diameters  upstream  of 
Pitot  tube 

•  The  distance  between  the  fan  and  flow  straightener  should  be  no  less  than 
3.5  duct  diameters. 

•  Duct  diameter  should  be  at  least  30  times  the  diameter  of  the  Pitot  tube. 


3.4.1  Design  Integration 

Using  the  guidelines  for  airflow  measurement  as  requirements  it  is  possible  to 
begin  the  development  of  putting  physical  dimensions  to  the  wind  tunnel  design.  The  first 
step  was  to  select  the  desired  size  of  the  duct,  since  many  of  the  requirements  for  duct 
length  are  dependent  on  the  hydraulic  diameter  of  the  duct.  There  are  three  factors 
driving  this  portion  of  the  design;  the  heat  sink  geometry,  the  diameter  of  the  fan  and  the 
allotted  space  within  the  environmental  chamber. 


17 

Approved  for  public  release;  distribution  unlimited. 


The  fan  diameter  is  3.0  inches,  any  duct  size  chosen  should  be  less  than  or  equal 
to  this  size.  Having  a  duct  size  greater  than  this  would  limit  the  maximum  air  flow  within 
the  duct  ultimately  lowering  the  total  dynamic  pressure  and  making  an  accurate  reading 
more  difficult.  Additionally,  any  duct  size  chosen  should  be  similar  to  the  overall  profile 
size  of  the  heat  sink  (1.813  in.  x  1.25  in.).  The  reason  for  this  is  that  all  of  the  air  flow 
will  be  directed  through  the  channels  of  the  heat  sink  (no  bypass  air).  This  is  to  insure  a 
more  accurate  approximation  for  air  velocity  over  the  heat  sink  surface,  measured  from 
the  Pitot  tube.  Having  no  bypass  air  means  an  approximation  would  not  need  to  be  made 
of  the  actual  average  velocity  through  the  heat  sink  channels.  In  a  bypass  air  situation  the 
air  over  the  heat  sink  is  slightly  less  than  the  average.  This  is  due  to  added  surface  friction 
and  boundary  layer  forces  dropping  a  greater  pressure  across  the  finned  section  of  the 
duct.  Instead,  with  no  bypass  air,  the  average  air  velocity  over  the  surface  can  easily  be 
determined  based  on  the  change  in  duct  cross-sectional  area.  The  duct  diameter  should 
also  closely  emulate  the  size  of  the  heat  sink  to  minimize  the  size  of  any  expansion  or 
contractions  that  will  need  to  happen  as  a  result  of  a  duct  size  change. 

The  last  and  most  difficult  constraint,  to  the  design  of  the  wind  tunnel,  was  the 
size  of  useable  space  inside  the  environmental  chamber.  A  measurement  of  this  space 
shows  the  base  of  21  in.  by  21  in.  and  25  in.  in  height.  It  is  desirable  that  the  duct  be  on  a 
flat  surface  to  allow  for  easy  construction  and  integration  of  all  other  necessary 
components.  Therefore,  all  of  the  length  requirements  much  fit  within  a  21  in.  x  21  in. 
square  on  the  inside  floor  of  the  environmental  chamber.  To  maximize  the  one¬ 
dimensional  length,  the  duct  was  oriented  from  one  corner  to  the  other  along  the  diagonal 
of  the  square  base.  This  gave  a  one-dimensional  length  of  29.7in.  However,  since  the 
wind  tunnel  was  three-dimensional  the  true  duct  length  was  less  than  this  since  it  would 
impact  the  walls  at  the  comers  in  the  environmental  chamber. 

Dictated  by  the  size  of  the  heat  sink,  a  2.0  in.  x  2.0  in.  square  duct  was  chosen. 
Calculations  of  the  minimum  lengths  governing  the  distance  between  certain  components 
of  the  wind  tunnel  were  then  carried  out  based  on  the  requirements  for  accurate  flow 
measurement.  Table  3-1  shows  the  computed  values  for  the  necessary  section  length  of 
the  duct  based  off  multiples  of  the  hydraulic  diameter,  DK 

The  hydraulic  diameter  is  an  attribute  commonly  used  in  describing  duct 
geometry  in  non-circular  ducts,  such  as  in  this  case.  It  acts  as  an  equivalent  dimension 
simplifying  the  alternate  duct  shape  to  as  if  they  were  circular.  This  allows  one  to  apply 
calculations  normally  reserved  for  circular  flows  to  any  geometry.  In  this  case  the  duct  in 
question  was  a  2.0  in.  square  duct,  where  the  hydraulic  diameter  is  represented  as, 

Equation  3-1 

4  ab 

Dh  = - 

2  (u  +  b) 


Where  a,  is  the  height  and  b,  is  the  width  of  the  duct.  Therefore,  the  hydraulic  diameter 
for  this  duct  is  2.0  inches. 
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Table  3-1  Critical  Requirements  for  Duct  Length 


Length  (in.) 

Element 

7.0 

Fan  to  Honeycomb  exit,  Settling  Distance  ~3.5  x  DH 

1.0 

Honeycomb  Size  (not  in  total  sum) 

10.0 

Honeycomb  to  Pitot  Tube,  Settling  Distance  -5.0  x  DH 

1.5 

Pitot  Tube  Size 

3.0 

Pitot  Tube  to  Heat  Sink,  Settling  Distance  ~1.5  x  DH 

4.0 

Heat  Sink  Size 

2.0 

Exit  Section  1.0  x  Dh 

27.5 

TOTAL 

From  Table  3-1,  it  was  determined  that  a  duct  length  of  27.5  in.  was  needed  to 
allow  for  accurate  airflow  measurement.  However,  this  only  accounts  for  the  wind  tunnel 
length  after  the  fan,  and  not  the  fan  (1.75  in.)  or  the  control  valve  (7.75  in.).  Ultimately, 
the  wind  tunnel  would  need  to  be  37  in.  long  in  order  to  fit  all  the  necessary  elements  of 
the  experiment.  This  is  8  inches  longer  than  the  environmental  chamber  would  allow,  in  a 
straight  line.  Decreasing  the  size  of  the  duct  (i.e.  Dh  )  to  just  the  height  and  width  of  the 
heat  sink  did  help,  but  was  only  able  to  decrease  the  duct  length  by  about  5.0  inches. 
Instead,  a  new  design  was  adopted  that  did  not  require  the  wind  tunnel  to  be  in  a  straight 
line.  Several  design  concepts  were  vetted  that  incorporated  turns  within  the  duct.  Figure 
3.7  shows  a  concept  where  the  wind  tunnel  turns  180°  back  under  itself. 


Figure  3.7  Wind  Tunnel  Design,  180°  Turn  Concept. 


The  valve,  fan  and  settling  length  after  the  fan  were  all  in  the  upper  section  of  the 
wind  tunnel  and  the  honeycomb,  settling  distance  for  the  flow  measurement,  Pitot  tube, 
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heat  sink  and  exit  were  in  the  lower  section.  This  decreased  the  longest  linear  section  of 
the  wind  tunnel  to  20.5  in.  which  worked  although  it  added  some  other  components  such 
as  turns  within  the  duct. 


3.4.2  Detailed  Design 

In  order  to  properly  determine  if  this  concept  would  work,  a  detailed  design  was 
developed  and  integrated  into  the  allotted  space  of  the  environmental  chamber,  using  a 
three  dimensional  (3D)  CAD  model.  By  taking  and  modeling  the  existing  components 
(fan,  valve,  Pitot  tube  and  heat  sink)  a  wind  tunnel  was  designed  around  these  parts  and 
their  settling  lengths.  The  duct  size  was  selected  to  be  a  2.0  in.  square  duct.  Fitting  the 
20.5  in.  test  section  into  the  environmental  chamber  proved  to  be  a  tighter  fit  than  it  did 
on  paper.  The  model  in  Figure  3.8  shows  how  the  integration  of  the  wind  tunnel  into  a  3- 
Dimensional  space  limited  how  close  the  ducting  could  be  to  the  two  opposing  comers  of 
the  environmental  chamber.  This  is  because  the  incorporation  of  the  turning  sections 
(elbows)  accounted  for  remaining  of  the  allocated  space. 


Figure  3.8  Detailed  CAD  Model  of  Wind  Tunnel  Design 


Figure  3.9  shows  the  final  configuration  of  the  wind  tunnel  with  the  length 
dimensions  displayed  for  the  spacing  between  crucial  sections. 
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Figure  3.9  Wind  Tunnel  Sectioned  View  with  Dimensions,  in  Inch. 


It  is  worth  noting  that  the  9  inches  separating  the  honeycomb  (Figure  3.10, 
Section  16)  and  inlet  to  the  Pitot  tube  does  not  meet  the  minimum  recommended  settling 
distance  of  10  in.  Due  to  fitting  other  critical  components  near  the  heat  sink  the  Pitot  tube 
had  to  be  moved  up  by  1  inch.  The  consequence  of  this  was  thought  to  be  minimal  given 
that  the  distance  it  is  moved  reduced  the  settling  distance  by  5%  and  the  5  x  Du  is  an 
approximation  of  the  necessary  settling  length. 


3.4.3  Wind  Tunnel  Section  Design 

Critical  to  the  functionality  of  the  wind  tunnel  was  its  ability  to  maximize  the 
airflow  over  the  heat  sink  test  section,  thus  allowing  for  the  widest  range  of  Reynolds 
numbers  tested.  To  achieve  this,  the  wind  tunnel  was  designed  to  minimize  the  pressure 
drop  experienced  by  the  flow.  The  exercise  of  minimizing  the  pressure  loss  is  a  well- 
documented  field  that  has  a  wide  range  of  academic  and  industrial  references.  Figure  3.10 
shows  the  elements  adopted  by  this  design  in  order  to  minimize  the  pressure  loss  within 
the  wind  tunnel. 
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Figure  3.10  Features  Used  to  Minimize  Pressure  Loss  in  Wind  Tunnel. 


The  majority  of  the  components  are  controlled  contractions  (nozzles)  or  expansions 
(diffusers)  within  the  duct.  The  other  two  major  elements  are  the  two  honeycomb  flow 
straighteners  and  the  two  90°  expanding  comers.  The  pressure  loses  within  a  pipe  flow 
can  be  separated  into  two  categories,  the  major  and  the  minor  losses.  The  total  pressure 
loss  or  system  resistance  ( Pioss )  is  found  by  a  sum  of  both  the  major  ( Vmajor )  and  minor 
losses  (p minor),  as  seen  by  Equation  3-2. 


Equation  3-2 


Pioss  ^  '  Pmajor  T  ^  '  Pmir 


The  major  loss,  or  head  loss,  in  a  pipe  flow  are  due  to  friction  with  the  moving 
fluid  against  the  walls  of  the  duct.  These  losses  are  dependent  on  the  velocity  of  the  flow 
(t/oo),  the  hydraulic  diameter  of  the  pipe  (Du,  Equation  3-1),  the  length  of  the  pipe  (L), 
the  density  of  the  air  (p),  and  the  friction  coefficient.  The  friction  coefficient  (/),  for 
laminar  flow,  was  estimated  by  the  following  equation. 

Equation  3-3  [17] 

64 

/  =  T> — '  Reo  £  3,000 

KeD 


For  turbulent  flow,  the  friction  coefficient  was  found  from  the  following  equation. 

Equation  3-4  [18] 

/  =  0.316 fleD~1/4;  3,000  <  ReD  <  2  x  104 
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The  Reynolds  number  ( ReD )  of  the  flow  is  based  on  the  hydraulic  diameter  within  the 
duct  element. 

Equation  3-5 

Dtf  U0 o  p 

ReD  =  - 

H 


The  major  losses  in  all  sections  of  the  duct  were  calculated  through  the  use  of  Equation 
3-6  [19]. 


Equation  3-6 


Vmajor  loss 


The  minor  losses  are  associated  with  changes  in  the  cross-sectional  area  of  the  duct, 
changes  in  velocity  or  direction  of  airflow  or  obstructions  in  the  air  stream.  For  each  case 
the  non-dimensional  change  in  pressure  is  referred  to  as  the  loss  coefficient,  p.  The  loss 
coefficient  is  related  to  the  pressure  drop  through  the  use  of  the  following  equation. 

Equation  3-7 

/ PU oo2 

Vminor  loss  V  l  n 


For  this  experiment,  the  minor  losses  were  either  from  a  cross-sectional  area 
change  or  change  in  direction  from  a  bend  in  the  wind  tunnel.  The  pressure  loss  from 
both  of  these  cases  is  computed  differently.  A  full  record  of  the  computations  for 
determining  the  minor  losses  for  every  section  within  the  wind  tunnel  is  presented  in 
APPENDIX  B. 

The  final  element  that  was  considered  when  evaluating  the  total  wind  tunnel 
pressure  loss  was  the  heat  sink.  This  was  estimated  through  a  method  developed  by 
Simons  [20]  for  predicting  the  pressure  drop  across  parallel  plate  heat  sinks.  This  method 
and  its  uses  are  further  detailed  in  Section  4.2.2  of  this  thesis. 


3.4.4  Pressure  Loss 

After  all  the  elements  of  the  wind  tunnel  were  designed  and  the  computations  for 
estimating  their  major  and  minor  losses  were  detennined,  it  was  necessary  to  compute  the 
total  system  loss,  or  system  resistance.  There  were  several  reasons  for  computing  the 
predicted  system  resistance  before  completing  construction  of  the  wind  tunnel.  The  first 
reason  was  for  the  selection  of  the  differential  pressure  transducers  for  measuring  flow 
velocity  and  static  pressure  drop.  These  sensors  were  relatively  expensive  and  it  was 
useful  to  know  the  exact  ones  that  were  needed  for  the  experiment  rather  than  using  a 
sensor  that  was  out  of  range  and  needing  to  purchase  more.  The  second  reason  was  to 
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develop  predictions  for  a  comprehensive  testing  matrix.  Knowing  the  capabilities  of  the 
wind  tunnel  (maximum  volumetric  flow  rate)  ahead  of  time  would  give  insight  to  the 
expected  airflow  conditions  over  the  heat  sink.  This  would  allow  for  testing  parameters  to 
be  developed  with  the  goal  of  maximizing  the  range  of  Reynolds  number  tested  over  the 
heat  sink.  Lastly,  this  information  was  used  to  validate  the  wind  tunnel. 

The  process  of  calculating  the  total  system  pressure  loss  was  fairly  simple.  The 
system  pressure  loss  is  a  sum  of  all  the  major  and  minor  losses  for  every  section 
throughout  the  wind  tunnel  (Equation  3-2).  This  will  yield  the  loss  coefficients  for  each 
section.  However  the  pressure  loss  is  also  dependent  on  the  velocity  of  the  flow,  as  seen 
in  Equation  3-7.  The  velocity  is  derived  from  dynamic  pressure  (Poy),  measured  at  the 
Pitot  tube,  shown  in  the  following  equation. 


Equation  3-8 


In  order  to  find  the  pressure  drop  within  the  wind  tunnel  both  the  density  of  the  air 
and  velocity  had  to  be  known.  This  was  achieved  by  plotting  the  wind  tunnel  system 
resistance  against  the  fan  curves  [14]  at  various  ambient  pressures  (i.e.,  air  densities), 
shown  in  Figure  3.11,  Figure  3.12,  Figure  3.13,  and  Figure  3.14.  The  intersection  of  the 
system  resistance  and  fan  flow  rate  were  found.  This  predicted  the  maximum  flow  rate  of 
the  system  for  a  given  fan  RPM. 
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The  values  found  at  the  intersection  point  are  shown  in  Table  3-2. 


Table  3-2  Values  for  the  Intersection  Points  of  the  Fan  Curves  and  System  Resistance 


Pressure 

(atm) 

Density 

(kg/m3) 

Fan  Speed 
(RPM) 

Flow  Rate 
(m3/s) 

Pressure  Head 
(inH20) 

1.00 

1.202 

25,000 

0.024 

2.59 

0.70 

0.918 

25,000 

0.023 

1.89 

0.50 

0.699 

25,000 

0.023 

1.52 

0.20 

0.326 

25,000 

0.011 

0.24 

0.20 

0.326 

30,000 

0.014 

0.36 

From  the  intersection  points  of  the  fan  curves,  both  the  volumetric  flow  rate  and 
pressure  head  were  determined.  Next,  the  local  velocity  at  every  point  within  the  duct  can 
be  determined  by  dividing  the  volumetric  flow  rate  by  the  cross-sectional  area.  A  table 
displaying  all  the  values  for  pressure  loss  as  well  as  the  flow  conditions  at  every  point 
within  the  duct  for  the  case  of  running  the  fan  at  25,000  RPM  at  1.0  atm  can  be  found  in 
APPENDIX  C.  By  knowing  the  velocity  within  the  duct  at  the  location  of  the  Pitot  tube 
and  over  the  heat  sink  it  was  possible  to  accurately  size  the  sensors  used  for  the 
experiment. 
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3.5  Test  Matrix 

Test  matrix  was  derived  by  finding  the  maximum  airflow  conditions  within  the 
duct  and  subdividing  it  into  five  20%  increments  based  on  velocity.  This  was  done  for  all 
cases  where  the  maximum  volumetric  flow  rate  was  detennined  using  the  fan  curve.  Then 
by  linear  interpolation  of  these  points,  estimations  were  made  for  the  maximum  airflow  at 
the  remaining  conditions,  (0.9,  0.8,  0.6,  0.4,  0.3,  and  0.1  atm).  Table  3-3  is  the  derived 
test  matrix  which  gives  estimation  for  the  flow  conditions  such  as  velocity,  Reynolds 
number,  dynamic  pressure,  and  static  pressure  throughout  the  entire  range  of  parameters 
that  this  test  was  interested  in.  The  values  listed  under  ’Atmosphere  and  Altitude’  were 
based  on  the  1976  Standard  Atmosphere  [7].  Since  the  experiment  was  not  conducted  at 
sea  level  the  corresponding  pressure  for  823  feet  ASL  was  used.  The  values  for  the 
dynamic  pressure  in  the  duct  were  used  to  size  the  differential  pressure  sensor  for  the  air 
velocity  measurement.  The  static  pressure  drop  across  the  heat  sink  was  also  used  to  size 
which  sensor  will  be  appropriate  for  this  measurement. 
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Table  3-3  Preliminary  Test  Matrix  for  Desired  Airflow  Conditions 


Atmosphere  and  Altitude 

Duct  Air  Flow 

Heat  Sink  Air  Flow 

Percent  of 

Maximum  airflow 

Atmospheric 

Pressure 

(atm) 

Density 

(kg/m3) 

Altitude 

(feet) 

Atmospheric 
Pressure  (psi) 

Kinematic 

Viscosity 

(m2/s) 

Velocity  (m/s) 

Reynolds 

Number 

Dynamic 
Pressure  (in 
of  H20) 

Velocity  (m/s) 

Reynolds 

number 

Static  Pressure 
Drop  (in  of 
H20) 

97% 

Atmosphere 

100% 

0.97 

1.196 

823 

14.26 

1.52E-05 

9.46 

31719 

0.2148 

26.36 

17481 

1.3878 

80% 

0.97 

1.196 

823 

14.26 

1.52E-05 

7.57 

25375 

0.1375 

21.09 

13985 

0.9317 

60% 

0.97 

1.196 

823 

14.26 

1.52E-05 

5.68 

19031 

0.0773 

15.82 

10489 

0.5603 

40% 

0.97 

1.196 

823 

14.26 

1.52E-05 

3.78 

12687 

0.0344 

10.54 

6993 

0.2766 

20% 

0.97 

1.196 

823 

14.26 

1.52E-05 

1.89 

6344 

0.0086 

5.27 

3496 

0.0855 

90% 

Atmosphere 

100% 

0.90 

1.125 

2885 

13.23 

1.61E-05 

9.33 

29426 

0.1966 

26.00 

16221 

1.2903 

80% 

0.90 

1.125 

2885 

13.23 

1.61E-05 

7.46 

23541 

0.1258 

20.80 

12977 

0.8673 

60% 

0.90 

1.125 

2885 

13.23 

1.61E-05 

5.60 

17656 

0.0708 

15.60 

9732 

0.5224 

40% 

0.90 

1.125 

2885 

13.23 

1.61E-05 

3.73 

11770 

0.0315 

10.40 

6488 

0.2585 

20% 

0.90 

1.125 

2885 

13.23 

1.61E-05 

1.87 

5885 

0.0079 

5.20 

3244 

0.0803 

80% 

Atmosphere 

100% 

0.80 

1.022 

6045 

11.76 

1.77E-05 

9.14 

26192 

0.1714 

25.48 

14442 

1.1545 

80% 

0.80 

1.022 

6045 

11.76 

1.77E-05 

7.31 

20954 

0.1097 

20.39 

11553 

0.7774 

60% 

0.80 

1.022 

6045 

11.76 

1.77E-05 

5.49 

15715 

0.0617 

15.29 

8665 

0.4695 

40% 

0.80 

1.022 

6045 

11.76 

1.77E-05 

3.66 

10477 

0.0274 

10.19 

5777 

0.2333 

20% 

0.80 

1.022 

6045 

11.76 

1.77E-05 

1.83 

5238 

0.0069 

5.10 

2888 

0.0730 

70% 

Atmosphere 

100% 

0.70 

0.918 

9540 

10.29 

1.97E-05 

8.95 

23032 

0.1476 

24.96 

12703 

1.0236 

80% 

0.70 

0.918 

9540 

10.29 

1.97E-05 

7.16 

18425 

0.0945 

19.97 

10162 

0.6907 

60% 

0.70 

0.918 

9540 

10.29 

1.97E-05 

5.37 

13819 

0.0531 

14.98 

7622 

0.4183 

40% 

0.70 

0.918 

9540 

10.29 

1.97E-05 

3.58 

9213 

0.0236 

9.98 

5081 

0.2088 

20% 

0.70 

0.918 

9540 

10.29 

1.97E-05 

1.79 

4606 

0.0059 

4.99 

2541 

0.0659 

60% 

Atmosphere 

100% 

0.60 

0.810 

13470 

8.82 

2.24E-05 

8.93 

20283 

0.1297 

24.91 

11186 

0.9266 

80% 

0.60 

0.810 

13470 

8.82 

2.24E-05 

7.15 

16226 

0.0830 

19.92 

8949 

0.6266 

60% 

0.60 

0.810 

13470 

8.82 

2.24E-05 

5.36 

12170 

0.0467 

14.94 

6712 

0.3806 

40% 

0.60 

0.810 

13470 

8.82 

2.24E-05 

3.57 

8113 

0.0208 

9.96 

4475 

0.1908 

20% 

0.60 

0.810 

13470 

8.82 

2.24E-05 

1.79 

4057 

0.0052 

4.98 

2237 

0.0608 

50% 

Atmosphere 

100% 

0.50 

0.699 

17970 

7.35 

2.59E-05 

8.913 

17460 

0.1114 

24.85 

9629 

0.8254 

80% 

0.50 

0.699 

17970 

7.35 

2.59E-05 

7.13 

13968 

0.0713 

19.88 

7703 

0.5597 

60% 

0.50 

0.699 

17970 

7.35 

2.59E-05 

5.35 

10476 

0.0401 

14.91 

5777 

0.3412 

40% 

0.50 

0.699 

17970 

7.35 

2.59E-05 

3.57 

6984 

0.0178 

9.94 

3852 

0.1720 

20% 

0.50 

0.699 

17970 

7.35 

2.59E-05 

1.78 

3492 

0.0045 

4.97 

1926 

0.0554 

40% 

Atmosphere 

100% 

0.40 

0.583 

23265 

5.88 

3.11E-05 

7.77 

12698 

0.0706 

21.66 

7003 

0.5690 

80% 

0.40 

0.583 

23265 

5.88 

3.11E-05 

6.22 

10158 

0.0452 

17.33 

5603 

0.3881 

60% 

0.40 

0.583 

23265 

5.88 

3.11E-05 

4.66 

7619 

0.0254 

13.00 

4202 

0.2385 

40% 

0.40 

0.583 

23265 

5.88 

3.11E-05 

3.11 

5079 

0.0113 

8.67 

2801 

0.1217 

20% 

0.40 

0.583 

23265 

5.88 

3.11E-05 

1.55 

2540 

0.0028 

4.33 

1401 

0.0402 

30% 

Atmosphere 

100% 

0.30 

0.462 

29775 

4.41 

3.92E-05 

6.57 

8515 

0.0401 

18.33 

4697 

0.3630 

80% 

0.30 

0.462 

29775 

4.41 

3.92E-05 

5.26 

6812 

0.0257 

14.67 

3758 

0.2496 

60% 

0.30 

0.462 

29775 

4.41 

3.92E-05 

3.94 

5109 

0.0144 

11.00 

2818 

0.1551 

40% 

0.30 

0.462 

29775 

4.41 

3.92E-05 

2.63 

3406 

0.0064 

7.33 

1879 

0.0805 

20% 

0.30 

0.462 

29775 

4.41 

3.92E-05 

1.31 

1703 

0.0016 

3.67 

939 

0.0275 

20% 

Atmosphere 

100% 

0.20 

0.326 

38385 

2.94 

5.56E-05 

5.231 

4779 

0.0179 

14.59 

2637 

0.1972 

80% 

0.20 

0.326 

38385 

2.94 

5.56E-05 

4.18 

3823 

0.0115 

11.67 

2109 

0.1374 

60% 

0.20 

0.326 

38385 

2.94 

5.56E-05 

3.14 

2868 

0.0064 

8.75 

1582 

0.0870 

40% 

0.20 

0.326 

38385 

2.94 

5.56E-05 

2.09 

1912 

0.0029 

5.84 

1055 

0.0465 

20% 

0.20 

0.326 

38385 

2.94 

5.56E-05 

1.05 

956 

0.0007 

2.92 

527 

0.0168 

10% 

Atmosphere 

100% 

0.1 

0.163 

52805 

1.47 

1.11E-04 

3.62 

1655 

0.0043 

10.11 

914 

0.0748 

80% 

0.1 

0.163 

52805 

1.47 

1.11E-04 

2.90 

1324 

0.0027 

8.09 

731 

0.0537 

60% 

0.1 

0.163 

52805 

1.47 

1.11E-04 

2.17 

993 

0.0015 

6.06 

548 

0.0354 

40% 

0.1 

0.163 

52805 

1.47 

1.11E-04 

1.45 

662 

0.0007 

4.04 

365 

0.0200 

20% 

0.1 

0.163 

52805 

1.47 

1.11E-04 

0.72 

331 

0.0002 

2.02 

183 

0.0079 

3.6  Experiment  Construction 


The  many  unique  sections  of  the  duct  were  built  using  several  methods.  Figure 
3.15  shows  the  entire  experiment  setup  inside  the  environmental  chamber. 
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Figure  3.15  Wind  Tunnel  inside  Environmental  Chamber. 


A  combination  of  commercial  off-the-shelf  (COTS)  items,  custom,  and  3D 
printed  parts  was  used  for  the  build.  The  constant  diameter  corner  and  a  friction  fitting  to 
attach  it  to  the  actuated  valve  were  purchased.  The  clear  acrylic  straight  square  duct  was 
built  from  a  flat  0.22”  thick  sheet.  The  intricate  components  that  could  not  be  purchased 
or  easily  fabricated  were  made  using  3D-Printing.  A  total  of  7  parts  were  made  using  3D 
printing,  shown  in  red  in  Figure  3.16. 
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Figure  3.16  3D  Printed  Parts  within  Wind  Tunnel 


Of  the  7  parts,  5  were  made  using  ABS  plastic  with  a  tolerance  estimated  to 
±0.035  in.  The  other  two  remaining  parts,  the  Heat  Sink  Housing  and  the  Insulation  Base 
Block,  were  made  using  polycarbonate  plastic  also  with  a  tolerance  of  ±0.035  in.  The 
decision  was  made  to  use  polycarbonate  (PC)  for  these  parts  because  of  its  higher 
working  temperature.  The  glass  temperature  of  PC  is  around  140°C.  This  is  20°C  higher 
than  the  temperature  rating  of  the  heaters.  This  would  ensure  that  heat  sink  will  never  get 
hot  enough  to  permanently  damage  the  housing.  Additionally,  all  of  the  3D  printed  parts 
that  came  in  contact  with  the  airstream  were  smoothed  to  minimize  friction  losses,  with 
the  exception  of  the  two  expanding  corners.  Sanding  the  inner  surface  of  these  parts  was 
made  overly  difficult  given  the  turning  vanes  obstructing  the  most  of  the  surface. 


3.6.1  Sensors:  Types,  Usage  and  Installation 

The  next  challenge  in  the  construction  of  the  wind  tunnel  was  the  installation  of 
all  the  sensors  used  to  operate,  control,  and  conduct  measurement  for  the  experiment. 
Three  different  pressure  measurements  are  needed  for  this  experiment;  dynamic  pressure 
using  the  Pitot  tube,  static  pressure  difference  across  the  heat  sink,  and  ambient  pressure. 
For  these  three  measurements  five  sensors  are  needed.  Both  the  dynamic  and  static 
pressure  drop  measurements  used  two  different  differential  pressure  sensors,  for  a  low 
and  high  range.  Because  this  experiment  covered  such  a  wide  range  of  pressure 
measurements,  the  uncertainty  of  these  measurements  was  minimized  by  using  two 
ranges  for  sensors.  As  seen  in  Table  3-3,  the  dynamic  measurement  could  see  a  pressure 
range  of  0.2  -  0.0007  inH20  and  a  static  pressure  range  of  1.4-0.02  inH20.  For  these  two 
measurements  Setra  Model  264  Differential  pressure  transducers  were  chosen.  The 
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ranges  chosen  for  both  measurements  were  0-0.1  &  0-0.25  inH20  for  the  dynamic 
measurement  and  0-1.0  &  0-2.5  inH20  for  the  static  pressure  measurement.  Switching 
between  either  sensor  was  easily  done  prior  to  an  experiment,  based  on  the  anticipated 
measurement  range.  The  third  pressure  measurement  needed  was  the  absolute  pressure 
within  the  environmental  chamber.  For  this  the  Setra  Model  552  absolute  pressure  sensor 
was  used.  It  was  capable  for  covering  the  entire  range  of  pressures  for  the  experiment,  0- 
15  PSIA.  All  five  of  these  sensors  required  an  excitation  voltage  of  9-30  VDC  and  output 
0-5VDC  for  their  full  scale  range.  They  were  all  mounted  within  the  environmental 
chamber  with  the  exception  of  the  absolute  sensor  that  was  relocated  to  a  fitting  attached 
to  the  outside  of  the  chamber.  This  was  due  to  an  unknown  error  with  the  sensor  output 
that  was  observed  when  operating  it  inside  the  chamber.  Figure  3.17  shows  how  the  four 
differential  sensors  were  mounted  within  the  experimental  setup  inside  the  environmental 
chamber. 


Figure  3.17  Mounting  Locations  of  the  4  Differential  Sensors. 


The  transition  of  changing  from  one  sensor  to  the  other  was  made  by  moving  the 
vinyl  tubing  that  connects  to  the  Pitot  tube  or  static  pressure  taps  from  one  sensor  to  the 
other.  The  location  of  the  measurement  devices  were  also  careful  placed.  As  discussed  in 
Section  3.4  the  location  of  the  Pitot  tube  is  placed  9.0  in.  downstream  of  the  honeycomb 
flow  straightener  at  the  center  point  of  the  duct.  Figure  3.18  shows  the  mounting  location 
for  the  two  static  pressure  taps  used  to  measure  the  difference  in  pressure  before  and  after 
the  heat  sink. 
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Figure  3.18  Location  of  Pressure  Taps  on  Heat  Sink  Housing 


A  total  of  16  Type-T  thermocouples  were  used  in  this  experiment.  Table  3-4  lists 
all  the  thennocouple  designations  and  their  location  within  the  experiment. 

Table  3-4  Thermocouple  Numbering 


TC  NUMBER 

MEASUREMENT  LOCATION 

TC01 

AIR  AT  INLET  1 

TC02 

AIR  AT  INLET  2 

TC03 

HEAT  SINK  BASE  1 

TC04 

HEAT  SINK  BASE  2 

TC05 

HEAT  SINK  BASE  3 

TC06 

HEAT  SINK  BASE  4 

TC07 

HEAT  SINK  BASE  5 

TC08 

HEAT  SINK  TIPI 

TC09 

HEAT  SINK  TIP  2 

TC10 

HEAT  SINK  TIP  3 

TC11 

AIR  AT  EXIT  1 

TC12 

AIR  AT  EXIT  2 

TC13 

HOUSING  EXTERNAL  1 

TC14 

FAN 

TC15 

HOUSING  EXTERNAL  2 

TC16 

AIR  AT  EXIT  3 

A  total  of  5  thermocouples  are  used  for  the  measurement  of  the  temperature  of  the 
heat  sink  base.  The  thermocouples  are  placed  at  the  five  locations  shown  in  Figure  3.19, 
with  TC03  being  at  the  upstream  most  location  and  TC06  being  closest  to  the  exit  of  the 
air  leaving  the  heat  sink. 
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Figure  3.19  Heat  Sink  Base  Thermocouple  Placement  in  Inch 


These  5  measurements  were  used  in  the  direct  calculation  for  the  thermal 
resistance,  which  was  used  to  find  the  experimental  heat  transfer  coefficient.  They  were 
also  used  to  observe  the  symmetry  and/or  uniformity  of  the  heat  source.  If  any  irregular 
or  unexpected  pattern  was  observed  it  would  be  a  concern.  To  ensure  the  quality  of  these 
measurements,  steps  were  taken  to  eliminate  measurement  error.  The  thermocouples  are 
located  on  the  base  of  the  heat  sink  just  above  the  interface  point  between  the  heater  and 
heat  sink.  To  avoid  any  contact  that  the  thermocouple  tips  would  have  with  the  heater, 
tiny  holes  were  drilled  at  each  placement  location.  Thin  channels  were  then  cut  for  the 
thermocouple  wire  to  be  run  to  the  rear  of  the  duct  so  that  when  the  heater  was  installed 
there  was  a  flat  surface  to  surface  contact  between  heater  and  the  aluminum  base.  Lastly, 
thermal  grease  was  used  to  fill  the  holes  at  the  thermocouple  tips  and  Kapton  tape  was 
used  to  cover  the  interface  and  create  a  barrier  between  heater  and  thermocouple,  and  to 
eliminate  air  pockets  that  may  have  existed.  This  method  of  attachment  was  also  used  in 
other  surface  thermocouple  placement  (e.g.,  the  tips  of  the  heat  sink,  the  fan,  and  sides 
and  bottom  of  the  housing). 

Figure  3.20  shows  the  relative  location  of  the  two  thermocouples  used  for 
measuring  the  temperature  of  the  air  entering  the  test  section  prior  to  encountering  the 
heat  sink. 
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Before  the  Heat  Sink 
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Figure  3.20  Thermocouple  Placement  at  the  Entrance  of  the  Test  Section,  Shown  as  Section 

17 


Those  two  stainless  steel  thermocouple  probes  were  placed  at  the  inlet  to  the  heat 
sink  housing,  roughly  1.5  inches  upstream  of  the  heat  sink.  The  air  within  the 
environmental  chamber  was  kept  at  a  very  uniform  temperature  by  means  of  an  internal 
fan  agitating  the  air.  Because  of  this,  very  little  temperature  variation  is  expected  of  the 
air  entering  the  heat  sink  inlet  and  only  two  thermocouples  were  thought  needed  for  this 
measurement. 


Figure  3.21  shows  a  cross-section  for  the  test  section  at  the  midpoint  for  the  heat 
sink,  2.0  inches  from  the  front  and  rear. 


Center  of  HeatSink 


TC15 


Figure  3.21  Thermocouple  Placement  at  the  Center  of  the  Heat  Sink 
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A  total  of  8  thermocouples  are  placed  at  this  location.  Of  these,  the  three  heat  sink 
base  thermocouples  have  been  discussed.  The  other  three  placed  at  the  tip  of  the  heat  sink 
fins,  designated  as  TC08,  TC09,  and  TC10,  were  used  to  give  an  inclination  on  the  heat 
sink  fin  efficiently.  Three  were  chosen  to  observe  any  variation  in  temperature  across  the 
tips.  Unfortunately,  one  of  these  thennocouples  was  damaged  during  installation  (TC08) 
and  repairing  it  was  considered  untenable  given  its  relatively  low  importance  and  high 
potential  for  damaging  other  thermocouples  during  the  process.  The  last  two 
thermocouples,  TC13  and  TC15,  were  later  added  to  the  experiment.  Their  purpose  was 
to  measure  the  housing  temperature  at  the  sides  and  base.  This  information  was  used  to 
determine  the  thermal  losses  by  conduction  into  the  housing. 


The  last  crucial  measurement  within  the  heat  sink  was  the  air  at  the  exit  of  heat 
sink.  This  measurement,  like  the  heat  sink  fin  tip  measurement,  was  not  absolutely 
critical  to  the  experiment  but  was  found  to  be  useful  in  post  experimental  analysis.  The 
placement  for  these  thermocouples  is  shown  in  Figure  3.22. 


Figure  3.22  Thermocouple  Placement  at  the  Exit  of  the  Heat  Sink,  Section  21 


Obtaining  a  measurement  for  the  average  air  temperature  at  exit  was  challenging. 
This  was  because  the  air  temperature  was  likely  not  unifonn,  due  to  surface  temperature 
variation,  thermal  boundary  layer  growth  and  asymmetric  mixing  of  the  air.  For  these 
reasons,  it  was  determined  that  three  thennocouples  centered  at  the  lower,  middle  and 
upper  third  of  a  single  channel  would  give  the  best  results  for  the  average  temperature  of 
air  leaving  the  center  channel  of  heat  sink. 

Lastly,  a  thermocouple  (TC14)  was  placed  on  the  stator  hub  of  the  fan  used  to 
drive  the  air  within  the  wind  tunnel  (Figure  3.23).  This  measurement  was  needed  for  the 
safe  operation  of  the  experiment.  In  some  cases  the  fan  was  run  beyond  its  spec  power 
limits  for  extended  periods  of  time.  It  became  necessary  to  monitor  the  temperature  of  the 
fan  and  implement  safety  stops  within  the  software  to  avoid  thermal  failure  and 
permanent  degradation  of  the  motor. 
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Figure  3.23  Fan  Thermocouple  Placement 


The  last  parameter  needing  recording  for  this  experiment  is  the  power  dissipated 
by  the  heater.  This  was  done  by  measuring  the  current  flow  to  the  heater  via  a  shunt 
resistor.  A  shunt  resistor  is  a  highly  stable  and  dependable  method  for  measuring  current 
for  a  constant  DC  voltage.  By  using  a  small  electrical  resistance  (0.0197Q)  in  the  circuit 
and  measuring  the  voltage  drop  across  the  resistor  the  current  can  be  found.  Equation  3-9 
was  used  to  determine  the  current  by  measuring  the  voltage  drop  across  a  known 
resistance. 


Equation  3-9 

T  _  ^cLrop 

‘ Heater  ^ 

K-e,  Shunt 

The  power  to  the  heater  is  therefore  I  Heater  Repeater,  where  Re, Heater  is  the  heater 
electrical  resistance. 


3.6.2  Wiring,  Power  Distribution  &  Programming 

After  all  the  instrumentation  was  selected  and  placed  properly  within  the  wind 
tunnel  it  became  time  for  the  final  wiring  of  all  sensors  to  the  DAQ,  distributing  power  to 
the  necessary  components,  and  programing  the  remote  operation  and  data  acquisition  of 
the  experiment.  The  wiring  and  methodology  for  power  and  signal  commands  are 
discussed  here  briefly.  More  details  are  available  in  the  fonn  of  a  wiring  diagram  in 
APPENDIX  D.  The  approach  taken  when  wiring  many  components  with  multiple  power 
sources  to  a  single  data  collection  source  was  to  allow  all  sources  to  act  as  isolated  power 
sources  with  a  single  common  ground  terminated  to  earth  at  the  DAQ.  A  grounding  bus 
was  made  to  do  this.  This  was  advantageous  for  a  number  of  reasons.  The  primary  was  to 

36 

Approved  for  public  release;  distribution  unlimited. 


eliminate  any  ground  loops.  Secondly,  this  ensured  that  all  measurements  have  the  same 
common  ground  voltage.  This  guarantees  that  there  was  no  electrical  measurement  offset 
errors. 

The  equipment  used  to  power  the  electronics  was  as  follows.  For  the  heater  and 
fan,  two  BK  Precision  series  1685B  switch-mode  power  supplies  were  used.  The  power 
supplies  were  controlled  by  a  0-5  VDC  source  corresponding  to  a  0-60VDC,  5amp  max 
output.  The  output  will  be  controlled  by  the  NI  (cRIO).  This  allowed  for  active  control  of 
the  fan  and  the  power  to  the  heater.  All  other  electronics  with  the  exception  of  the  cRIO 
were  powered  by  a  common  source,  a  laptop  switching  power  supply  outputting  19.5 
volts,  3.0  amp  maximum.  This  provided  all  5  pressure  transducers  with  the  necessary 
excitation  voltage,  power  for  the  shunt  isolated  amplifier,  and  the  operating  power  for  the 
actuated  valve.  The  cRIO  was  powered  by  a  benchtop  DC  power  supply. 

Since  all  of  these  components  were  being  operated  remotely  within  the  vacuum 
chamber,  wiring  needed  to  pass  into  the  chamber  through  a  small  circular  hole  on  the  side 
while  maintaining  an  airtight  seal.  For  this,  a  custom  seal  was  designed.  The  seal  was 
made  of  a  %”  aluminum  plate  with  holes  cut  for  the  5  wire  bundles  that  passed  through. 
This  was  more  advantageous  than  one  large  hole  as  less  room  temperature  vulcanization 
(RTV)  silicone  was  needed  to  fill  the  air  gap  between  aluminum  and  wire  insulation. 
Figure  3.24  shows  the  design,  attachment,  and  sealing  of  the  wire  pass-through  plate. 


Figure  3.24  Wire  Pass-Through  Seal  Design 


The  last  step  in  the  functionality  of  the  experiment  was  the  software  programing 
for  control  and  data  acquisition.  For  this,  NI  Lab  VIEW  2014  was  used.  The  software 
allowed  for  the  creation  of  a  customized  virtual  instrument  (VI)  for  the  NI  cRIO  chassis 
and  modules  to  execute  desired  commands  based  on  various  settings  and  inputs.  The  goal 
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of  this  programing  was  the  remote  operation  of  experimental  commands  for  airflow  and 
heater  power  as  well  as  recording  of  data.  The  program  was  given  a  target  for  airflow 
(dynamic  pressure)  and  heater  power  (voltage  drop  across  the  shunt),  and  control  loops 
within  the  VI  manipulated  the  fan  speed  (voltage  to  the  fan),  valve  position  (0-10VDC 
signal  voltage)  and  voltage  to  the  heater  (0-60 V). 

The  original  methodology  for  controlling  airflow  was  to  run  the  fan  at  a  high 
speed  and  use  only  the  valve  to  regulate  airflow.  However,  this  concept  was  abandoned  to 
avoid  running  the  fan  in  high  fatigue  conditions  for  longer  than  required  periods  of  time. 
Instead,  the  following  logic,  shown  in  Figure  3.25,  was  used  to  achieve  the  desired  air 
flow  within  the  wind  tunnel. 


Figure  3.25  Air  Flow  Control  Loop  Operational  Logic  Diagram 


In  summary,  this  control  scheme  had  two  operating  modes.  The  control  loops  for 
both  of  these  operating  modes  as  well  as  the  hardware  involved  are  shown  in  Figure  3.26. 
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Figure  3.26  Wind  Tunnel  Control  Diagram 


The  starting  mode  or  default  mode  (valve  control)  used  the  valve  to  regulate 
airflow.  The  fan  was  run  at  16,000  RPM,  and  the  valve  would  open  from  the  closed 
position,  increasing  the  airflow.  If  this  control  scheme  was  able  to  reach  the  airflow  set 
point  then  the  program  would  actively  regulate  it  to  make  sure  it  stays  on  target.  This  was 
needed  since  the  fan  required  an  active  PID  loop  to  keep  it  at  constant  speed,  as  it  would 
lose  efficiency  when  the  temperature  increased.  The  other  control  scheme  (fan  control) 
was  activated  when  the  valve  reached  the  fully  open  position  and  the  airflow  target  was 
not  yet  met.  When  this  condition  took  place,  the  outer  control  loop  would  alter  the  fan 
speed  set  point,  increasing  it,  until  the  target  was  achieved. 

During  any  operation,  the  Lab  VIEW  VI  also  had  safety  measures  in  place.  If  the 
temperature  of  the  heat  sink  reached  120°C,  the  program  cut  power  and  cooled  it  down. 
Also,  if  the  fan  temperature  reached  50°C,  its  voltage  was  decreased  to  a  safe  range  and 
the  valve  was  fully  open  for  cooling. 

In  addition  to  controlling  the  airflow  and  heater  power,  the  Lab  VIEW  VI  was  also 
tasked  with  data  collection.  For  the  purpose  of  this  experiment,  it  was  necessary  to  take 
measurements  in  a  steady  state  condition  where  the  temperatures  of  the  heat  sink  and 
airflow  were  no  longer  changing  under  constant  testing  conditions.  The  approach  to 
achieve  this  was  to  have  the  data  logging  be  initiated  by  a  trigger  that  monitors  the  rate  of 
change  for  certain  inputs  that  were  desired  to  be  steady  (i.e.  the  heater  base  temperature 
TC05)  and  the  flow  velocity  ( Pdyc ))•  This  was  accomplished  by  measuring  the  standard 
deviation  averaged  over  the  last  100  points  over  10  seconds  (100  samples/10  sps  =  10 
seconds).  When  the  standard  deviation  fell  below  a  target  value,  data  recording  was 
started.  The  program  would  record  10  seconds  of  data  at  10  sps.  After  this  was  done,  it 
would  then  cool  the  heat  sink  to  40  °C  before  executing  its  next  test.  It  would  do  this  for 
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as  many  testing  conditions  were  input  and  then  shut  down.  More  infonnation  on 
detennining  what  targets  were  adequate  to  achieve  steady  state  is  in  section,  3.7.5. 


3.7  Validation  and  Calibration 

Before  testing  could  begin,  it  was  necessary  to  calibrate  all  the  instrumentation 
used  in  the  experiment  as  well  as  characterizing  certain  aspects  of  the  test  setup  in  order 
to  ensure  accurate  operation  of  the  wind  tunnel  that  will  yield  reproducible  results. 


3.7.1  Pressure  Sensors 

A  total  of  five  sensors  were  used  in  this  experiment:  two  for  dynamic  pressure, 
two  for  static  pressure  drop,  and  one  for  absolute  pressure.  All  five  sensors  were 
purchased  from  Setra  and  were  calibrated  by  the  manufacturer.  The  results  of  the 
calibration  are  in  APPENDIX  E.  The  results  of  the  calibration  were  used  to  generate 
linear  regression  equations  for  each  of  the  pressure  transducers.  These  equations  were 
used  to  correct  the  voltage  output  from  the  sensor  to  the  true  pressure  on  the  sensor 
diaphragm,  within  the  Lab  VIEW  VI.  Also,  the  standard  error  of  the  regression  equations 
was  documented  and  used  in  the  uncertainty  analysis  for  the  experiment  (APPENDIX  H). 


3.7.2  Thermocouples 

The  thermocouples  are  Type-T  probes  purchased  from  Omega.  They  were 
calibrated  by  a  resistive  temperature  detector  (RTD)  in  a  silicone  oil  bath.  The  calibration 
was  performed  in  the  range  of  temperature  0-150  °C.  The  details  of  this  calibration  are 
in  APPENDIX  F.  The  results  of  the  calibration  were  used  to  generate  linear  regression 
equations  for  each  of  the  16  thennocouples.  These  equations  were  used  to  correct  the 
thennocouple  reading  within  the  Lab  VIEW  VI  to  the  true  temperature.  This  correction 
was  made  during  the  data  recording  process.  The  standard  error  of  the  regression 
equations  was  documented  and  used  in  the  uncertainty  analysis  for  the  experiment 
(APPENDIX  H). 

3.7.3  Duct  Air  Flow  Calibration 

When  using  one  Pitot  tube  for  velocity  measurement,  it  is  important  to  make  sure 
that  the  measured  pressure  at  the  Pitot  tube  location  could  be  used  to  accurately  represent 
the  average  dynamic  pressure  within  the  entire  duct  cross-section.  If  the  Pitot  tube  is 
placed  in  the  center  of  the  duct,  the  measured  velocity  will  be  higher  than  the  actual 
average  velocity  due  to  boundary  layer  effect  on  the  velocity  profde.  Figure  3.27  shows 
that  a  boundary  layer  begins  to  form  decreasing  the  velocity  at  the  edges  and  increasing 
the  velocity  at  the  center.  In  the  experiment,  only  the  velocity  at  the  center  of  the  duct 
will  be  measured,  it  is  important  to  establish  a  relationship  between  the  centerline 
velocity  and  the  average  velocity  over  the  experimental  range  of  Reynolds  number. 
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Figure  3.27  Developing  Flow  between  Two  Plates  [21]. 


By  generating  a  constant  flow  within  the  duct  and  moving  the  Pitot  tube  to 
different  locations  and  recording  the  dynamic  pressure,  it  is  possible  to  know  the  velocity 
profile  for  that  particular  flow  condition,  thus  the  average  velocity  as  well  as  the 
centerline  velocity.  To  do  this  a  9  point,  3x3,  grid  pattern  was  chosen.  Figure  3.28  shows 
the  location,  in  green  dots,  for  the  Pitot  tube  transverse  measurement. 


Figure  3.28.  Location  of  Transverse  Measurement  Points  within  Duct  Cross-Section 


The  next  difficulty  was  how  to  conduct  this  measurement  for  all  the  flow 
conditions  possible  throughout  this  experiment  (i.e.  Reynolds  number  range).  The 
experiment  greatly  varied  in  velocity  and  air  pressure  throughout  testing,  both  of  which 
affect  the  dynamic  pressure  reading.  To  mimic  the  changes  to  testing  conditions  for  the 
transverse  measurement  using  the  vacuum  chamber  was  not  possible,  on  account  of  the 
fact  that  relocating  the  Pitot  tube  after  each  measurement  would  require  the  constant 
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opening  and  closing  of  the  chamber  door,  making  this  process  too  time  consuming. 
Instead,  what  was  decided  was  to  create  a  velocity  profile  for  the  entire  range  of 
Reynolds  numbers  used  in  the  experiment.  This  was  accomplished  by  conducting  the 
calibration  outside  the  chamber  at  ambient  pressure  for  823  feet  ASL  and  only  varying 
the  velocity  within  the  duct  to  achieve  the  desired  Reynolds  number.  One  adverse 
consequence  of  running  this  test  at  ambient  pressure  was  that  the  velocity  had  to  be  very 
low  in  order  to  achieve  the  lower  Reynolds  numbers,  much  lower  than  if  they  were 
reached  at  higher  altitude  conditions.  Because  of  this  the  test  was  limited  by  the 
‘calibration  limit’  of  the  dynamic  pressure  sensors.  This  limit  was  found  to  be  near  Re  ~ 
1,500  for  the  0-0.1  inH20  sensor.  This  was  not  ideal  since  the  experiment  would  involve 
Reynolds  numbers  lower  than  1,500.  However,  the  test  was  still  able  to  capture  the 
transition  from  developing  to  fully  developed  flow.  Since  enough  information  was  taken 
within  the  developing  region  above  1,500,  a  strong  characterization  for  that  region  was 
developed  and  extrapolated  to  cover  lower  Reynolds  number  scenarios.  The  duct  velocity 
corrections  are  shown  in  Figure  3.29. 


Duct  Velocity  Adjustment  Based  on  the  Average  Reynolds  Number 


Figure  3.29  Air  Velocity  Correction  Curve 


This  plot  was  created  by  measuring  how  much  greater  the  center  point  Pitot  tube 
measurement  was  than  the  average  flow  and  plotting  it  as  a  function  of  the  Reynolds 
number.  This  velocity  correction  curve  was  useful  in  understanding  the  condition  of  the 
flow  within  the  duct  for  any  given  Reynolds  number.  More  importantly  it  also  quantifies 
the  difference  in  the  measurement  taken  by  the  Pitot  tube  at  the  center  of  the  duct  and  the 
average  velocity  of  the  airflow.  From  Figure  3.29  two  distinct  curves  immerge.  The 
dotted  red  line  at  Re  =  10,000  divides  the  developing  flow  curve  and  fully  developed 
turbulent  flow  curve.  The  developing  region  was  represented  by  a  power  curve  best  fit  to 
the  five  measurements.  The  fully  developed  turbulent  flow  was  represented  by  a  zero 
slope  line  since  boundary  layer  growth  was  minimized  due  to  turbulence.  The  corrections 
in  Figure  3.29  were  used  to  derive  the  relationship  between  Reynolds  number  at  the 
center  point  (ReDuctcP)  and  the  average  Reynolds  number  (ReDuctAvg),  as  shown  in 
Equation  3-10  (developing)  and  Equation  3-11  (fully  developed  turbulent). 
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Equation  3-10 


Re 


Duct.Avg 


ReDuct,cp  221.69 
099 


ReDuct,CP  <  10,000 


Equation  3-11 

R^Duct  cp  42.65 

Re Duct.Avg  —  1~02  '  Re  Duct, cp  —  10,000 


Equation  3-10  is  mathematical  correction  for  the  center  point  Reynolds  number  to 
the  true  average,  below  Re  =  10,000,  for  developing  flow,  and  Equation  3-11  is  the 
correction  for  the  Reynolds  numbers  above  10,000  in  turbulent  fully  developed  flow.  The 
corrections  were  used  to  correct  the  dynamic  pressure  measurements  after  each  test. 


3.7.4  Heater  Power 

The  next  calibration  necessary  for  the  operation  of  the  experiment  was  for  the 
measurement  of  power  dissipation.  This  experiment  uses  two  resistive  thin  film  heaters.  It 
was  possible  to  measure  the  thermal  energy  by  directly  measuring  the  electrical  energy 
dissipated  through  the  resistors.  The  way  to  measure  the  power  to  the  heater  was  to 
measure  its  resistance  and  record  the  current  in  the  circuit  for  the  heater.  Its  total 
resistance  is  measured  to  be  9.76  Q.  The  current  was  measured  using  a  calibrated  shunt 
resistor,  0.0197  Q.  By  measuring  the  voltage  across  the  shunt,  the  total  current  flowing 
through  the  circuit  was  found.  Since  this  current  will  be  equal  for  all  components  within 
the  closed  circuit  the  power  can  be  found  from  P  —  I2Re.  Doing  a  calibration  of  the 
entire  circuit,  from  the  shunt  through  the  isolated  amplifier  and  to  the  cRIO  was  the  most 
accurate  way  to  characterize  this  measurement.  Figure  3.30  shows  the  circuit  diagram  for 
the  calibration  setup. 


Figure  3.30  Wire  Diagram  for  Shunt  Current  Measurement  Calibration 


The  red  wire  circuit  is  the  simulated  power  to  the  heater.  A  load  resistor  was  used 
rather  than  the  actual  heater.  A  current  controlled  power  supply  provided  the  power  and 
for  reference,  a  Fluke  87  multimeter  was  used  to  monitor  current.  The  blue  wire  was  the 
current  sensing  circuit.  It  detected  the  voltage  drop  across  the  shunt  and  amplified  it  by  a 
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factor  of  100  before  it  was  recorded  by  the  cRIO.  A  Fluke  8846A  multimeter  was  used  to 
measure  the  voltage  differential  being  recorded  by  the  cRIO.  By  comparing  the  reading 
from  the  multimeter  and  the  cRIO  reading  a  correction  was  made  within  the  Lab  VIEW 
VI  so  that  it  reflected  the  actual  current  of  the  load  circuit.  Figure  3.31  shows  the  results 
from  this  calibration. 


0.000  1.000  2.000  3.000  4.000  5.000  6.000 

Applied  Current  (Amps) 


Figure  3.31  Shunt  Current  Calibration  Regression  Curve 


A  linear  regression  was  performed  to  develop  an  equation  of  best  fit,  within  the 
LabVIEW  VI,  that  would  translate  the  acquired  voltage  at  the  cRIO  into  current  to  the 
heater.  Table  3-5  displays  the  results  of  the  linear  regression  that  was  used  within  the 
LabVIEW  VI,  as  well  as  the  standard  error  that  is  used  in  the  uncertainty  analysis 
(APPENDIX  H). 

Table  3-5  Results  from  Linear  Fit  of  Shunt  Calibration 


Curve  Fit  From  Shunt  Current  Calibration 

Slope  (Amps/VDC) 

0.500542 

Intercept  (Amps) 

0.003375 

Standard  Deviation  (Amps) 

0.001534 

Standard  Error  (±Amps) 

0.003065 

3.7.5  Steady  State  Verification 

A  method  was  development  to  verify  that  a  steady  state  condition  has  been 
achieved  during  the  experiment  and  then  to  trigger  data  logging.  As  stated  in  section 
3.6.2,  the  LabVIEW  VI  would  calculate  the  standard  deviation  of  the  last  100  samples 
taken  over  10  seconds  for  both  air  velocity  and  heat  sink  temperature.  The  calculation  for 
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the  running  standard  deviation  is  an  equivalent  metric  to  the  rate  of  change  of  a 
measurement.  If  the  SD  of  air  velocity  and  heat  sink  temperature  fell  below  a  threshold 
then  both  could  be  considered  steady  state.  The  heat  sink  temperature  steady  state 
measurement  used  TC05.  The  air  velocity  steady  state  measurement  used  the  dynamic 
pressure  taken  at  the  Pitot  tube.  Both  standard  deviation  measurements  were  used  to 
quantify  the  steady  state  condition  of  the  system  so  that  the  program  could  repeatedly 
self-determine  when  to  log  data. 

The  next  step  was  to  detennine  what  target  standard  deviation  would  be 
appropriate  to  signify  that  the  system  has  reached  steady  state.  If  the  target  was  too  large, 
then  data  logging  would  be  trigged  before  the  system  had  truly  reached  steady  state.  If  it 
was  too  small,  the  program  may  not  trigger  at  all.  To  find  the  correct  values  for  the  target 
standard  deviation,  the  following  process  was  used.  It  began  by  selecting  reasonable 
targets  for  the  dynamic  pressure  and  temperature,  (Heater  SD  Target:  0.02,  Dynamic 
Pressure  SD  Target:  0.0002).  Four  tests  were  run  with  these  parameters  and  the  results 
are  shown  in  APPENDIX  G.  Two  were  run  at  0.2  atm  with  high  and  low  velocity  and 
two  were  run  at  0.8  atm  with  high  and  low  velocity.  This  process  was  thought  to  give  a 
good  range  representing  the  envelope  of  all  testing  conditions.  The  observations  from 
these  tests  found  that  the  standard  deviation  of  the  airflow  needs  to  be  relaxed  for  low 
altitude  conditions  (Target:  0.00025)  and  can  be  lowered  for  high  altitude  conditions 
(Target:  0.00015).  This  was  due  to  the  use  of  two  different  differential  pressure  sensors 
for  the  dynamic  pressure  measurement.  Overall  the  airflow  was  able  to  reach  and  hold 
steady  state  significantly  quicker  than  the  heater  temperature.  In  addition  to  this  the  heat 
sink  temperature  still  appeared  to  have  a  slight  slope  during  data  recording.  Therefore, 
further  investigation  into  the  heat  sink  temperature  settling  time  is  needed.  To  investigate 
this  further  it  was  necessary  to  study  all  of  the  thennocouple  readings  and  determine  their 
steadiness  relative  to  the  TC05  measurement.  Figure  3.32  shows  the  results  from  an 
entire  test  run  at  0.2  atm,  with  0.0029  inH20  dynamic  Pressure,  40  W  to  the  heater  and 
having  targets  of:  Heater  SD  Target  =  0.015  °C;  Flow  SD  Target  =  0.00015  inH20. 
These  conditions  were  the  most  extreme  when  considering  settling  time  to  reach  steady 
state. 
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The  test  runs  from  the  initial  start  to  when  the  target  standard  deviations  are  met 
and  the  data  is  logged.  Table  3-6  displays  the  results  from  the  logged  data,  showing  the 
100  point  average,  standard  deviation  and  the  slope. 


Table  3-6  Temperature  Measurements  for  Settling  Time 


TC01 

AIRJN 

TC02 

AIRJN 

TC03 

HSB01 

TC04 

HSB02 

TC05 

HSB03 

TC06 

HSB04 

TC07 

HSB05 

TC15 

Bottom 

TC09 

HST02 

TC10 

HST03 

TC11 

AIREXIT01 

TC12 

AIREXIT02 

TC16 

AIREXIT03 

TC13 

Side 

AVG  (°C) 

25.67 

25.69 

88.44 

87.92 

96.07 

87.75 

92.49 

39.54 

85.14 

85.27 

40.65 

35.32 

42.06 

69.11 

STD  DEV(°C) 

0.0167 

0.0164 

0.0167 

0.0187 

0.0159 

0.0220 

0.0211 

0.0369 

0.0263 

0.0261 

0.0862 

0.0707 

0.0523 

0.0232 

Slope  (°C/sec) 

0.0015 

0.0023 

0.0033 

0.0047 

0.0037 

0.0056 

0.0061 

0.0120 

0.0079 

0.0072 

0.0254 

0.0210 

0.0162 

0.0068 

The  slope  was  important  in  determining  if  the  temperature  of  the  heat  sink  was 
still  increasing  or  had  leveled  off.  A  zero  slope  means  it  had  reached  steady  state.  As 
stated  previously,  because  of  noise,  achieving  a  zero  slope  is  not  feasible.  Therefore,  it 
was  desirable  for  the  slope  be  as  small  as  possible. 

Table  3-6  shows  that  at  every  instance  the  slope  was  positive.  This  includes  the 
air  entering  the  heat  sink.  This  means  that  the  temperatures  of  the  heat  sink  and  ambient 
air  were  still  rising.  The  explanation  for  this  is  that  the  environmental  chamber  was  very 
close  to  an  isolated  system;  therefore  any  thermal  energy  tends  to  remain  within  the 
chamber.  Hence  the  temperature  increases.  The  impact  of  this  was  that  the  slope  of  TC05 
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would  never  reach  zero  during  a  test.  However,  the  slope  of  TC05  would  approach  zero 
relative  to  the  increase  in  ambient  temperature,  which  could  also  be  considered  steady 
state.  The  difference  in  the  rate  of  change  between  TC05  and  the  air  at  inlet  will  be 
denoted  as  Ad.  So,  it  was  desirable  to  select  the  best  target  standard  deviation  that  would 
yield  a  minimum  A 6  and  could  be  reasonably  met.  From  Table  3-6,  the  difference  in 
slope  from  the  inlet  air  (average  of  TC01  and  TC02)  and  TC05  was  Ad  —  0.0018  °C / 
sec.  So  to  verify  that  having  a  Heater  SD  Target  of  0.015  °C  would  be  sufficient  one 
final  test  was  run  with  extended  time  to  observe  what,  if  any,  improvements  could  be 
made  to  Ad.  Three  more  tests  were  run  to  investigate  how  the  temperature  difference,  of 
the  heat  sink  base  (TC05)  and  inlet  air  (TC01  &  TC02),  changes  over  a  period  of  time 
after  the  target  standard  deviations  have  been  met.  Each  test  was  run  for  about  200%  of 
its  normal  settling  time  given  the  Heater  SD  Target  of  0.015  °C.  The  results  (APPENDIX 
G)  from  the  extended  time  tests  were  compared  with  data  from  previous  tests  which  used 
the  conventional  triggering  method  (Table  3-7) 

In  the  far-right  column  is  the  recorded  temperature  difference  for  each  test.  By 
averaging  the  standard  runs  and  subtracting  it  from  the  extended  time  run,  a  value  for  the 
increase  in  temperature  difference  over  the  extended  test  was  found.  For  the  two  0.8  atm 
tests,  the  difference  is  negligible.  In  the  case  of  the  0.2  atm  test  the  difference  was 
somewhat  significant.  However,  when  viewed  in  the  context  of  the  additional  time 
needed  to  observe  this  increase,  the  true  impact  on  the  experiment  was  minimal, 
~0.10  °C/min.  Therefore,  it  was  concluded  that  target  standard  deviation  of  0.015  °C  for 
TC05  was  satisfactory  to  use  to  allow  the  heat  sink  to  reach  steady  state  conditions. 


Table  3-7  Comparison  of  the  Temperature  Difference  for  Extended  Time  Testing 


Heater  SD  Target  (°C) 

-  or  - 

Additional  Time  after  SD 
Target  was  met  (min) 

Pressure 

(atm) 

Flow  % 

Mean  Absolute 
Pressure  (PSI) 

Mean  Dynamic 
Pressure 
(inH20) 

Mean  Heater 
Power  (inH20) 

Mean  Delta  T 
(°C) 

0.015 

0.8 

80% 

11.55 

0.10970 

60.05 

25.98 

0.015 

0.8 

80% 

11.47 

0.10980 

59.89 

25.76 

0.015 

0.8 

80% 

11.49 

0.10962 

59.84 

26.18 

(Extended  by  3:25) 

0.8 

80% 

11.50 

0.10986 

59.83 

25.94 

Increase  in  Delta  T  over  Extended  test: 

-0.03 

0.015 

0.8 

40% 

11.57 

0.02707 

59.86 

38.91 

0.015 

0.8 

40% 

11.49 

0.02777 

59.70 

37.54 

0.015 

0.8 

40% 

11.49 

0.02722 

59.66 

38.74 

(Extended  by  6:40) 

0.8 

40% 

11.50 

0.02782 

59.64 

39.70 

Increase  in  Delta  T  over  Extended  test: 

1.30 

0.015 

0.2 

40% 

2.90 

0.00281 

39.44 

68.52 

0.015 

0.2 

40% 

2.95 

0.00295 

39.40 

66.28 

(Extended  by  27:20) 

0.2 

40% 

2.89 

0.00269 

39.39 

70.25 

Increase  in  Delta  T  over  Extended  test: 

2.85 
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By  running  the  tests  for  additional  time  after  the  heater  target  was  met,  the  amount 
of  additional  temperature  rise  could  be  observed.  Since  the  results  show  that  the 
temperature  rise  was  negligible,  it  was  detennined  that  SD  <0.015  °C  was  a  good  target 
to  signify  heat  sink  had  reached  steady  state. 


3.8  Testing 

The  experiments  were  carried  out  closely  following  the  test  matrix  outlined  in 
section  3.5.  Testing  consisted  of  first  verifying  the  experiment  functionality  and  then 
sealing  the  chamber  and  taking  it  to  the  desired  vacuum  (i.e.,  altitude).  Once  the  desired 
pressure  was  reached  and  observed  to  be  stable  (the  chamber  was  not  airtight  and 
required  a  pump  to  be  running  at  all  times)  testing  began.  At  any  given  pressure,  a  total  of 
4-5  tests  were  done  at  varying  heater  power  and  flow  velocity.  Heater  power  inputs  were 
mostly  chosen  conservatively  to  limit  high  temperatures.  Velocities  were  chosen  to  cover 
a  complete  range  of  Reynolds  numbers  for  each  pressure  condition.  This  was  done  by 
using  the  maximum  velocity  and  incrementally  reducing  it  by  20%  for  each  subsequent 
test.  In  practice,  this  plan  worked  relatively  well.  However,  some  changes  were  made  for 
specific  velocities  of  some  measurements.  In  some  cases  the  maximum  velocity  was 
dictated  by  the  fan  speed.  It  was  too  strenuous  for  it  to  maintain  high  speeds,  necessary 
for  the  target  maximum  of  the  test  matrix,  for  long  periods  of  time  during  some  of  the 
higher  chamber  pressures.  This  issue  was  thought  to  be  caused  by  the  added  resistance  at 
the  exit  of  the  wind  tunnel  created  by  its  proximity  to  the  walls  of  the  environmental 
chamber.  For  this  reason,  the  issue  was  not  observed  during  the  airflow  validation 
process.  The  added  resistance  meant  the  fan  would  have  to  run  at  speeds  higher  than 
25,000  RPM  to  hit  its  target  velocity.  Thus,  the  velocity  targets  for  some  cases  were 
adjusted. 
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CHAPTER  4 

EXPERIMENTAL  RESULTS  AND  DISCUSSION 


The  experiment  focused  on  the  heat  sink  pressure  drop,  heat  transfer  coefficient, 
and  predictions  of  these  measurements.  The  predictions  were  derived  from  standard 
correlations.  Additionally,  the  validity  of  these  predictions  was  supported  by  accounting 
for  heat  loss  through  a  thennal  energy  balance  and  experimental  error  through  an 
uncertainty  analysis.  A  discussion  of  the  results  is  provided  in  this  chapter. 


4,1  Heat  Transfer  Coefficient 

The  primary  measurement  for  this  experiment  was  the  heat  transfer  coefficient.  The 
heat  transfer  coefficient  for  this  heat  sink  was  detennined  from  the  following 
relationship. 

Equation  4-1 

u  —  Qconv 

tLava  ~  A  AT 

^base1-*1  base 


where  Qconv  is  the  thennal  energy  transfened  from  the  heat  sink  to  the  air,  Abase  is  the 
area  of  the  heat  sink  base,  and  ATbaseis  the  temperature  difference  between  the  heat  sink 
and  air.  However,  it  is  more  useful  to  study  the  HTC  in  its  non-dimensional  fonn,  the 
Nusselt  number,  Nu  (Equation  4-2).  Presenting  the  heat  transfer  with  the  Nusselt  number 
would  make  the  results  comparable  to  existing  Nusselt  number  correlations. 

Equation  4-2 

L* 

NUavg  —  h-avg  T 

K-fluid 


The  Nusselt  number  is  the  average  heat  transfer  coefficient  multiplied  by  the 
characteristic  length  of  the  heat  sink  base  (L*)  over  the  thennal  conductivity  of  the  fluid 
(k fluid)-  The  thennal  conductivity  was  found  using  known  thermo-physical  properties  of 
air.  For  simplicity,  the  heat  sink  was  assumed  to  be  at  uniform  temperature.  To  estimate 
this  temperature,  an  average  of  five  heat  sink  base  measurements  (TC03,  TC04,  TC05, 
TC06  &  TC07)  was  averaged  with  an  average  of  the  heat  sink  fin  tip  temperatures,  (TC09 
&  TC10).  The  estimate  for  the  fluid  temperature  was  found  by  taking  an  average  of  the 
two  inlet  air  thennocouples  (TC01  &  TC02).  The  calculation  for  QC0nv  required  a 
thennal  energy  balance. 


4,1,1  Convective  Thermal  Energy  Calculation 

A  thermal  energy  balance  was  used  to  calculate  thermal  energy  imparted  to  the 
fluid  through  convection,  by  subtracting  losses  through  conduction  ( Qcond )  ar|d  radiation 
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{QRad)  from  the  electric  power  input  to  the  heater.  This  began  with  the  measurement  for 
electrical  power  dissipation  by  the  heater,  ( Pfieater )•  At  steady  state  conditions,  the 
electrical  power  to  the  heater  was  equal  to  the  total  thermal  energy  dissipated  by  the  heat 
sink.  Thermal  energy  was  dissipated  by  the  heat  sink  through  convection,  conduction, 
and  radiation.  Therefore,  knowing  the  amount  of  energy  transmitted  to  the  fluid  required 
estimates  for  the  losses  in  energy  through  conduction  and  radiation.  Equation  4-3 
represents  the  energy  balance  and  gives  the  convected  thermal  energy. 

Equation  4-3 

Qconv  —  P Heater  ~  ( Qcond  4”  QRad ) 


To  find  the  losses  due  to  conduction,  an  evaluation  of  the  heat  sink  test  section 
was  perfonned.  Figure  4.1  diagrams  the  primary  surface  to  surface  contact  of  the  heat 
sink  and  heater  to  the  polycarbonate  housing.  These  areas  were  along  the  bottom  through 
the  fiberglass  insulation  and  on  the  two  sides  of  the  heat  sink. 


Polycarbonate  3D  Printed  Plastic 


Fiberglass  insulation 


Figure  4.1  Thermal  Energy  Conduction  Losses 


The  conductive  thermal  energy  loss  was  estimated  using  one  dimensional  (ID) 
conduction.  Equation  4-4  is  the  formula  for  estimating  the  ID  conduction  losses  from  the 
heat  sink  into  the  housing.  Measurements  for  the  A Tside  HS  and  A TBottom  HS  were  found 
using  different  averages  for  thermocouples  near  the  surface  in  question.  A TBottom  HS  was 
from  the  difference  of  Tavg  Base  &TC  15.  The  Tavg  Base  is  simply  the  average  of  the  five 
thermocouples  imbedded  in  the  base  of  the  heat  sink  (TC03,  TC04,  TC05,  TC06  & 
TC07).  ATside  HS  was  the  difference  between  T  and  TC 13.  The  T  uc 

was  found  using  an  average  of  the  measurements  of  the  two  base  side  thermocouples 
(TC04  &  TC06).  This  was  then  averaged  with  an  average  of  the  two  thermocouple  tip 
measurements  (TC09  &  TC10).  The  following  equation  was  then  fonnulated  for  Qcond. 
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Equation  4-4 


Qcond  —  ^ 


AT 


Side  HS 


+ 


AT 


Bottom  HS 


Rt,Housing_Side  R t, Housing  _B  ottom 


Because  of  a  limitation  of  available  thermocouples  only  two  TCs  were  used  for 
the  external  temperature  measurements.  For  the  sides,  TCI 3  was  placed  at  the  center  of 
the  conduction  area.  Because  of  symmetry  the  same  temperature  was  expected  to  be  on 
the  alternate  side  of  the  heat  sink.  The  same  methodology  was  used  for  the  bottom,  by 
placing  TC15  in  the  center  of  the  base  conduction  area  (Figure  4.1).  Next,  the  thermal 
resistances  for  the  bottom  and  sides  will  need  to  be  computed. 

Equation  4-5  and  Equation  4-6  were  used  to  calculate  the  thennal  resistances  of 
the  housing  side  and  bottom.  The  resistances  were  both  sums  of  the  material  conduction 
resistance  and  conduction  interface  resistance. 

Equation  4-5 

n  _  Aside  1 

^t, Housing  _Side  ~  Z  Z  '  ZZT  Z 

KPC/iHS_Side  iLrLc,SideRHS_Side 


Equation  4-6 


R 


■‘FF 


t, Housing  JSottom 


k FF^HS_Bottom 


+ 


-‘Bottom 


^PC^HS_B  ottom 


+  3 


2  hr 


c,Bottom^HS_B  ottom 


where  Lside  is  the  thickness  of  the  housing,  kPC  is  the  thermal  conductivity  of  the 
polycarbonate,  and  AHSSide  is  the  contact  area  of  the  heat  sink  side  to  the  housing.  For 
the  bottom  of  the  heat  sink,  LFF  is  the  thickness  of  the  Fiberfrax  insulation,  LBottom  is  the 
thickness  bottom  of  the  housing,  AHSBottom  is  the  contact  area  for  the  bottom  of  the  heat 
sink,  and  kFF  is  the  thennal  conductivity  of  the  Fiberfrax.  To  estimate  the  interface 
resistance  a  value  for  the  interface  conductance  was  needed.  Of  the  interface  between  the 
aluminum  heat  sink  and  polycarbonate  housing  an  aluminum  to  aluminum  conductance 
[18]  value  was  used  for  hc  Side.  Of  the  interface  between  the  Kapton  heater  to  Fiberfrax, 
Fiberfrax  to  Fiberfrax,  and  Fiberfrax  to  polycarbonate  a  single  value  of  a  ceramic  to 
ceramic  interface  conductance  [18]  was  used  for  hc  Bottom. 

Table  4-1  and  Table  4-2  show  the  variables  and  results  of  the  thennal  resistance 
calculations.  The  housing  side  had  a  resistance  of  5.98  +  1.60  (m2K/W)  and  the 
housing  bottom  was  found  to  be  72.36  +  12.26  (m2K/W).  The  uncertainty  of  both 
calculations  comes  from  variability  of  the  values  in  Table  4-1  and  Table  4-2  and  was 
found  using  the  methods  detailed  in  APPENDIX  H. 
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Table  4-1  Thermal  Resistance  Calculation  for  Housing  Side 


Material 

Resistance 

Type 

Thermal 

Conductivity 

(W/m-K) 

Length 

(m) 

Area 

(m2) 

Interface 

Conductance 

(W/m2K) 

Resistance 

(K/W) 

k 

L 

f^HS_Side 

^c,Side 

R« 

Aluminum  (AL) 

Conduction 

209 

N/A 

0.003 

N/A 

AL/PC 

Interface 

0.003 

7100 

0.02 

Polycarbonate  (PC) 

Conduction 

0.205 

0.00787 

0.003 

5.95 

Total  Resistance  (K/W), 

Rt,Housing_Side 

5.98 

Table  4-2  Thermal  Resistance  Calculation  for  Housing  Bottom 


Material 

Resistance 

Type 

Thermal 

Conductivity 

(W/m-K) 

Length 

(m) 

Area 

(m2) 

Conductance 

(W/m2K) 

Resistance 

(K/W) 

k 

L 

^HS_Bottom 

^c,  Bottom 

R. 

Aluminum  (AL) 

Conduction 

209 

N/A 

0.005 

N/A 

Al/Fiberfrax 

Interface 

0.005 

1750 

0.12 

Fiberfrax  (FF) 

Conduction 

0.0462 

0.00635 

0.005 

29.38 

Fiberfrax  /Fiberfrax 

Interface 

0.005 

1750 

0.12 

Fiberfrax  (FF) 

Conduction 

0.0462 

0.00635 

0.005 

29.38 

Fiberfrax  /PC 

Interface 

0.005 

1750 

0.12 

Polycarbonate  (PC) 

Conduction 

0.205 

0.01270 

0.005 

13.24 

Total  Resistance  (K/W), 

^t, Housing  _Bottom 

72.36 

The  next  loss  contribution  considered  was  the  losses  due  to  radiation,  QRad-  As 
mentioned  in  Section  2.1,  at  high  altitudes,  radiation  should  not  be  ignored.  The  analysis 
of  radiation  heat  transfer  was  carried  out  by  dividing  the  surface  of  heat  sink  into  two 
different  parts,  which  could  be  easily  computed.  The  first  part  of  radiation  exchange  was 
the  front  and  rear  profiles  of  the  heat  sink  (left  side  of  Figure  4.2)  with  the  ambient.  The 
second  part  was  surface  to  surface  exchange  from  the  area  within  the  heat  sink  channels 
to  the  ceiling  of  the  housing  (right  side  of  Figure  4.2). 
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Figure  4.2  Heat  Sink  Radiation  Surfaces 


The  radiation  from  the  front  and  rear  profiles  was  considered  to  be  a  radiation 
emitted  to  a  blackbody  (i.e.,  100%  absorption  by  the  environment),  since  there  was 
sufficient  distance  to  any  objects,  reflection  would  be  minimal.  Equation  4-7  is  for 
radiation  thermal  energy  emission  from  the  front  and  rear  profiles. 

Equation  4-7 

Qr,P  =  2 Apas(Ts4  -  Tamb 4) 

Here,  Ap  is  one  of  the  transmitting  surface  area,  a  is  Stefan-Boltzmann  Constant 
(5. 6704x1  O'8),  £  is  the  emissivity  of  unfinished  aluminum  (0.04-0.07)  [18],  Ts  is  the  heat 
sink  average  surface  temperature  and  Tamb  is  the  temperature  of  the  ambient 
environment.  Tamb  was  estimated  using  the  inlet  air  temperature  measurement. 

Figure  4.2  (right)  shows  the  emitting  area  of  the  heat  sink  channels.  This 
calculation  was  slightly  more  complex  than  that  of  the  profile  radiation.  Because  of  the 
proximity  of  the  channels  to  the  ceiling  of  the  housing  this  needed  treated  as  a  radiation 
exchange  between  two  surfaces.  To  simplify  the  calculation,  the  exchange  from  one  heat 
sink  channel  surface  was  evaluated  and  then  multiplied  by  the  five  channels  to  get  the 
total  exchange.  Equation  4-8  is  for  radiation  exchange  between  two  surfaces. 

Equation  4-8 

o(T^-T2*) 

Ql~2  l-£l  1  1  —  e2 

2l1F1_2 


The  subscript  1  denotes  a  property  of  the  aluminum  heat  sink  single  channel  and 
the  subscript  2  denotes  the  properties  of  the  polycarbonate  housing  viewed  by  the 
channel.  Here,  the  heat  sink  surface  temperature,  7j,  has  already  been  computed.  The 
housing  temperature  was  estimated  by  an  average  of  the  inlet  air  (TC01  &  TC02)  and  the 
exit  air  temperature  (TC11,  TC  12  &  TC16).  In  the  denominator  of  Equation  4-8,  there 
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are  three  terms  representing  the  radiation  resistance.  The  areas  A1  &  A2  were  measured 
using  the  CAD  model  of  the  heat  sink.  The  emissivity  of  the  polycarbonate  was  estimated 
to  be  0.84.  The  view  factor  of  surface  1  to  2,  y41F1_2»was  simplified.  Working 
backwards,  the  view  factor  from  the  flat  surface  of  the  housing  to  the  aluminum  channel 
(F2-i)  only  sees  the  aluminum  and  therefore  100%  of  its  radiation  is  transmitted  to  the 
aluminum  (i.e.,  F2_t  —  1).  Because  of  the  relationship,  ^iFi-2  =  A2F2_1,  ^iF1_2  was 
simplified  to  be  represented  as  A2.  Therefore,  the  resistance  simplifies  to  Equation  4-9. 

Equation  4-9 

1  —  £i  1  1  —  £  2 

^  —  _ _  _| _ | _ _ 

£\A-l  A2  ^2^2 


Throughout  the  entire  experiment,  the  losses  from  any  given  test  averaged  -2.18 
±  0.96  W.  This  varied  with  altitude  as  well,  with  more  losses  at  higher  altitude  and  less  at 
lower  altitude. 


4.1.2  Nusselt  Number 

After  calculating  the  convective  heat  transfer  from  the  heater  to  the  passing  fluid 
the  Nusselt  number,  as  a  function  of  Reynolds  number,  could  found  through  Equation 
4-2.  Figure  4.3  shows  the  experimental  results.  Error  bars  were  added  based  on  the 
uncertainty  of  measurements.  For  details  of  how  the  uncertainty  was  calculated  refer  to 
APPENDIX  H. 


Average  Nusselt  Number  vs.  Reynolds  Number 
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Figure  4.3  Measured  Nusselt  Number  vs.  Reynolds  Number  at  Various  Altitude  Pressure 
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Figure  4.3  shows  each  individual  test  taken  for  different  atmospheric  pressures 
denoted  by  the  different  colors  in  the  legend.  It  is  clear  that  although  multiple  tests  were 
run  at  various  pressures,  ranging  from  1.0  -  0.1  atm,  the  Nusselt  number  for  each  pressure 
falls  within  a  single  curve.  Therefore,  the  conclusion  can  be  made  that  the  Nusselt 
number  is  independent  of  fluid  pressure  or  density,  and  it  is  only  a  function  of  Reynolds 
number. 


4.1.3  Nusselt  Number  Flat  Plate  Correlation 

It  was  useful  to  evaluate  accepted  methods  for  predicting  the  Nusselt  Number  and 
compare  them  against  the  measurements  taken.  The  first  comparison  was  done  on  flat 
plate  laminar  flow  correlation.  Predicting  a  heat  sink’s  perfonnance  based  on  this 
correlation  is  a  widely  accepted  practice  in  industry. 

Equation  4-10  [18] 

Nul  =  0.664 Re^^Pr1/3-,  Laminar,  Aver  age,  Tf,  0.6  <  Pr  <  50 


This  correlation  is  for  the  average  Nusselt  number  of  a  flat  plate  with  parallel 
laminar  flow.  The  Reynolds  number  ( ReL )  is  based  on  the  length  scale  of  the  flat  plate. 
To  extend  this  correlation  to  a  parallel  plate  heat  sink  surface  with  heat  transfer 
enhancement  feature,  an  efficiency  (e0)  was  applied  to  the  thermal  resistance  of  the  heat 
sink  surface,  to  represent  the  surface  enhancement. 

Equation  4-11 

1 

Rt,HS  —  J— 

^0  ™basen-L 


Therefore,  the  correlation  for  the  enhanced  surface  becomes: 


Equation  4-12 

Nuava  ~  Nu,60  ( — 
avg  L  °\WHSl*) 

where  Awet  is  the  total  area  of  the  heat  sink  surface  in  contact  with  the  air  flow  and  WHS 
is  the  length  of  the  heat  sink  along  the  direction  of  the  flow.  The  efficiency  e0  was 
detennined  from  a  correlation  for  predicting  the  surface  efficiency  of  a  fin  array  [18] 
shown  in  the  following  equation. 


Equation  4-13  [18] 


1  - 


NAfjn 

^wet 


(i  -«,) 


where  Afin  is  the  area  of  a  single  fin,  N  is  the  number  of  fins,  and  £f  is  the  efficiency  of  a 
single  fin.  6f  is  estimated  by  the  following  equation. 
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Equation  4-14  [18] 

tanh  mLc 

6f  - - 

‘  mLc 


where  mLc  is  a  correlation  for  estimating  the  temperature  distribution  within  the  heat  sink 
cross  section  (AHS  Profile). 


mLc 


Equation  4-15  [18] 


2  h 


1/2 


AL 


JiALAHS_Profile) 


3/2 


This  helps  tailor  the  correlation  for  different  surface  enhancements.  Figure  4.4  is 
Equation  4-12  against  the  experimental  data.  It  is  clear  that  this  predictive  method  falls 
short  of  characterizing  the  experimental  observations.  There  are  several  reasons  for  this. 
The  first  being  disagreement  of  the  actual  intent  of  this  correlation  and  the  experimental 
setup. 


Comparison  of  the  Average  Nusselt  Number  to  the  Flat  Plate  Correlation 
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Figure  4.4  Comparison  of  the  Average  Nusselt  Number  of  Flat  Plate  Correlation  vs.  Test 

Data 


The  application  of  this  correlation  is  really  meant  for  a  heat  sink  within  a  free 
flowing  fluid.  The  experiment,  by  design,  did  not  allow  free  flow  over  the  heat  sink  but 
instead  forces  the  fluid  through  the  FIS  channels,  eliminating  any  bypass  air.  In  a  free 
stream  a  large  portion  of  the  air  passing  by  the  heat  sink  will  take  the  path  of  least 
resistance  passing  over  the  heat  sink  geometry.  Although  this  correlation  is  not  applicable 
to  this  setup  it  does  offer  insight  into  the  performance  difference  between  external  flow 
forced  convection  and  confined  no  bypass,  or  internal  flow. 
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It  is  clear  from  the  comparison  of  this  correlation  to  the  observed  test  results  in 
Figure  4.4  that  the  channeled  flow  shown  as  test  data  has  much  higher  heat  transfer 
perfonnance  than  external  flow  shown  in  the  correlation  curve.  This  is  because  for  a 
bypass  air  TMS  the  air  will  take  the  path  of  least  resistance  and  tend  to  pass  over  the  heat 
sink  rather  than  through  the  channels.  Another  observation  is  that  higher  Reynolds 
number  flow  will  have  a  higher  percentage  of  bypassed  air  resulting  in  larger  difference 
in  Nusselt  number. 


4.1.4  Nusselt  Number  Internal  Flow  Correlation 

The  next  method  of  prediction  was  internal  flow  correlation  which  is  more 
applicable  to  the  experiments.  Several  correlations  exist  for  describing  the  heat  transfer 
between  the  channel  wall  and  the  fluid  flowing  through  it  and  the  appropriate  form 
depends  on  the  Reynolds  number  of  the  flow  within  the  channel. 

If  the  flow  is  laminar,  then  its  Reynolds  number  ( ReD )  is  below  2,300.  In  the 
instance  of  laminar  flow,  special  care  must  be  taken  to  account  for  the  thermal  boundary 
layer  growth  as  a  result  of  heat  transfer  at  the  surface  of  the  heat  sink,  decreasing  the 
HTC.  Figure  4.5  shows  dependency  of  the  local  Nusselt  number  to  thermodynamically 
developing  laminar  flow. 


4.36 

3.66 


Figure  4.5  Local  Nusselt  Number  Thermal  Development  [18] 


At  the  thermal  entry  region  of  a  flow  within  a  channel  the  Nusselt  number  can  be 
4-5  times  higher  than  it  is  after  it  has  become  fully  developed.  The  length  of  the  thermal 
developing  region  is  a  function  of  the  Reynolds  number,  Prandtl  number  and,  diameter 

f  x  RePr\ 

\d~  ~  ~~2o — )  *n  Edition  to  the  thermal  boundary  layer,  there  also  exists  a  hydrodynamic 
boundary  layer.  The  length  of  the  hydrodynamic  entry  region  is  a  function  of  the 

(x  Rb\ 

—  =  — ).  The  hydrodynamic  boundary  layer 

D#  20/ 

propagates  from  the  channel  walls  and  grows  larger  within  the  duct  cross-section,  for 
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laminar  flow.  This  further  decreases  HTC  by  slowing  done  the  motion  of  the  fluid  at  the 
channel  surface.  Both  the  thermal  and  hydrodynamic  boundary  layers  were  developing 
within  the  channels  but  their  impact  on  the  overall  heat  transfer  was  less  at  lower 
Reynolds  number  due  to  the  decrease  length  of  the  developing  region.  The  correlation 
detailed  in  Equation  4-16  [18]  is  for  internal  laminar  flow  with  a  combined  entry  region 
both  thermal  and  hydrodynamic  developing. 


Equation  4-16  [18] 


ReD 


DH  '  U  oo 
V 


<  2,300 


For  flow  that  was  not  laminar  (Re>  4,000)  a  turbulent  correlation  was  needed.  For 
a  fully  developed  turbulent  flow  within  a  channel  the  presence  of  a  thermal  boundary 
layer  has  less  of  impact  on  the  heat  transfer  than  with  laminar  flow.  This  is  due  to  the 
decreased  size  of  the  boundary  layer  for  a  fully  developed  turbulent  flow.  The  boundary 
layer  is  confined  to  the  area  at  the  surface  of  the  duct  walls.  However,  just  like  with 
laminar  flow,  at  the  entry  region  of  the  heat  sink  the  fluid  is  not  hydrodynamically  fully 
developed.  The  length  of  the  hydrodynamic  entry  region  is  a  function  of  the  hydraulic 

(X  \  .  X 

—  w  10).  With  —  =  9.5  the  flow  through  the  channels  never  reached 
dh  /  dh 

fully  developed  turbulent  during  this  experiment.  Unfortunately,  correlations  for 
predicting  the  Nusselt  number  within  the  entry  region  of  a  developing  turbulent  flow  are 
not  as  prevalent.  In  the  absence  of  the  desired  correlation,  the  Nusselt  number  for  a  fully 
developed  turbulent  flow  was  used.  The  correlation  expressed  in  Equation  4-17  is  a  well- 
supported  correlation  for  the  average  Nusselt  number  of  turbulent,  fully  developed  flow. 

Equation  4-17  [18] 

Nud  =  0.023 Re4/5Pr03;  ReD  >  10,000 


Figure  4.6  is  a  comparison  of  the  measured  average  Nusselt  number  with  that  of 
correlative  predictions  for  the  Nusselt  number  in  internal  laminar  flow  having  combined 
entry  region  (Equation  4-16  [18])  and  in  fully  developed  internal  turbulent  flow 
(Equation  4-17  [18]). 
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Comparison  of  the  Average  Nusselt  Number  to  Correlations 
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Figure  4.6  Comparison  of  the  Average  Nusselt  Number  and  Internal  Flow  Correlations 


Figure  4.6  shows  that  Equation  4-17  [18]  has  better  agreement  with  the  measured 
date  than  that  of  Equation  4-16,  for  their  applicable  range.  Equation  4-17  follows  the 
measured  values  reasonably  well  with  exception  of  Reynolds  number  below  10,000,  as 
expected.  Below  Reynolds  number  1,500  Equation  4-16  describes  the  measurements 
well. 


4.2  Pressure  Loss 

The  next  focus  of  the  results  and  discussion  is  the  pressure  loss  over  the  heat  sink. 
Theoretically  the  pressure  loss  is  related  to  HTC  and  may  be  used  to  describe  the 
performance  of  a  heat  sink.  However,  knowledge  of  the  pressure  loss  can  be  used  to  size 
both  heat  sink  and  fan. 

For  this  experiment  the  measurement  of  the  pressure  loss  was  made  directly,  by 
two  static  pressure  taps  before  and  after  the  heat  sink  and  a  single  differential  pressure 
sensor.  It  is  useful  to  study  this  measurement  using  the  non-  dimensional  pressure  loss 
coefficient  77.  This  is  done  by  dividing  the  static  pressure  loss  by  the  average  fluid 
momentum  term,  -j)  UHS2 . 

Equation  4-18 

_  APstc 

”  =  'P77) 
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Figure  4.7  shows  the  independence  of  the  pressure  loss  coefficient  and  density 
when  plotted  against  Reynolds  number.  An  uncertainty  analysis  was  conducted  and  error 
bars  were  added  to  Figure  4.7.  The  error  of  this  measurement  is  relatively  small  except 
for  at  low  Reynolds  numbers.  A  more  complete  description  of  the  uncertainty  analysis  is 
in  APPENDIX  H. 
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Figure  4.7.  Comparison  of  the  Pressure  Loss  Coefficient  at  Different  Pressure 


It  is  useful  to  compare  the  measured  results  with  correlations  for  the  pressure  loss 
coefficient. 


4.2.1  Internal  Flow  Friction  Factor  Correlation 

The  typical  approach  to  predicting  the  pressure  loss  over  a  given  geometry,  such 
as  this,  is  to  use  a  correlation  that  has  been  established  to  describe  the  loss  given  the 
course  of  the  airflow  and  mode  of  pressure  loss.  For  instance  there  are  two  types  of  losses 
taking  place  that  contribute  to  the  total  static  pressure  loss  across  the  heat  sink,  the  major 
and  minor  losses.  The  minor  losses  are  from  changes  to  the  cross-sectional  area  of  the 
flow,  such  as  enlargements  or  contractions.  In  the  case  of  the  heat  sink,  there  was  a  minor 
loss  for  the  abrupt  contraction  in  cross-sectional  area  at  the  entrance  to  the  heat  sink, 
shown  in  Equation  4-19  and  a  loss  associated  with  the  rapid  expansion  in  cross-sectional 
area  at  the  exit  of  the  heat  sink,  shown  in  Equation  4-20. 

Equation  4-19 

(  DlV 

Sudden  Contraction-,  riMinor  —  0.42  II  —  —  0.221 
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Equation  4-20 


2 


Sudden  Expansion; 


dMinor  — 


=  0.276 


The  minor  loss  coefficient  from  the  sudden  contraction,  Equation  4-19,  and 
sudden  expansion,  Equation  4-20,  is  entirely  geometry  dependent  and  therefore  will 
remain  constant  throughout  the  entire  range  of  air  velocity  tested. 

Unlike  the  minor  losses,  the  major  losses  are  dependent  on  the  velocity  of  the 
flow.  Correlations  for  predicting  the  major  loss  is  shown  in  Equation  4-21  where  the 
friction  factor,  /,  is  Reynolds  number  dependent. 

Equation  4-21 

L 

V Major  ~  f  71 

uHy 


Two  correlations  for  the  friction  factor  were  chosen.  Both  are  well  established 
correlations.  Equation  4-22  [18]  [17]  is  for  laminar  fully  developed  flow  and  Equation 
4-23  [18]  is  for  fully  developed  turbulent  flow. 

Equation  4-22  [18]  [17] 

64 

/  =  - — ;  ReD  <  3,000 

tieD 

Equation  4-23  [18] 

/  =  0.316fleD~1/4;  2,300  <  ReD  <  20,000 


A  comparison  of  these  two  correlations  to  the  measured  results  is  presented  in  Figure  4.8. 
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Comparison  of  the  Pressure  Loss  Coefficient  to  Correlations 
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Figure  4.8  Comparison  of  the  Pressure  Loss  Coefficient:  Measurement  vs  Correlation 


The  most  obvious  observation  from  this  chart  is  the  apparent  failure  of  the  fully 
developed  laminar  flow  correlation  to  predict  the  loss  coefficient.  The  disagreement  here 
is  profound  enough  to  be  considered  unusable  for  this  application.  There  are  two  possible 
contributions  for  the  disagreements  observed:  the  major  loss  correlation  under-predicts 
the  loss  coefficient  or  the  approximations  for  the  minor  losses  under-predict  the  losses. 
Since  the  minor  loss  prediction  method  remains  constant  throughout  all  cases,  an  under¬ 
prediction  of  this  would  simply  mean  a  negative  offset  of  the  curve.  This  could  be  the 
case  since  the  curve  for  the  measured  losses  in  the  laminar  region  fit  the  same  arc  as  the 
correlation.  Another  explanation  for  this  can  be  found  when  comparing  the  intended  use 
of  the  minor  loss  approximation  to  the  case  of  this  experiment.  The  intent  of  the 
approximation  is  meant  for  the  simpler  situation  of  a  circular  pipe  of  certain  diameter 
suddenly  experiencing  a  change  in  that  diameter.  It  does  not  necessarily  account  for 
complex  geometry  within  the  flow  path  (i.e.,  heat  sink).  The  other  explanation  for  the 
disagreement  would  be  the  correlation  for  the  friction  factor.  This  possibility  is  explored 
further  in  the  next  section,  where  an  enhanced  method  of  adapting  the  friction  factor  for  a 
parallel  plate  heat  sink  is  explored. 


4.2.2  Friction  Factor  Simons  Method 

A  paper  by  Simons  [20],  details  an  analytical  method  for  predicting  the  pressure 
drop  across  a  parallel  plate  heat  sink.  The  method  is  similar  to  the  previously  outlined 
procedure  of  summing  the  major  and  minor  losses.  The  computation  for  the  minor  losses 
is  identical  to  the  process  in  Section  4.2.1.  The  only  difference  is  how  the  major  loss  is 
found.  Therefore,  by  exploring  Simons  method  a  better  understanding  was  achieved  on 
the  shortcomings  of  the  previous  predictive  methods.  In  other  words,  answer  to  the 
question,  is  the  issue  with  the  minor  or  major  loss  prediction,  would  be  found. 
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Simons  method  takes  an  already  existing  friction  factor  correlation,  like  the  one 
detailed  by  Equation  4-22,  and  uses  it  to  build  an  ‘apparent  friction  factor,’  fapp  for  the 
entire  heat  sink  major  loss,  shown  in  Equation  4-24. 


fapp 


Equation  4-24 


ReD  <  2,300(est) 


The  apparent  friction  factor  uses  the  characteristic  length  of  the  heat  sink  (L*),  the 
average  Reynolds  number  through  the  channel,  and  a  friction  factor  fitted  (ffu)  to  a  fifth 
order  polynomial,  shown  in  Equation  4-25.  Note  that  the  equation  is  only  applicable  in 
the  hydrodynamically  developing  laminar  flow.  The  upper  limit  for  ReD  in  Equation 
4-24  was  estimated  based  on  the  experimental  data  within  the  paper. 

Equation  4-25 

24  -  32.5/  +  46.7/2  -  40.8/3  +  23/4  -  6.1/5 


The  apparent  friction  factor  is  related  to  the  total  pressure  loss  coefficient  through 
the  following  equation. 

Equation  4-26 

ZL 

V Minor  T  4/jpp  — 


Figure  4.9  is  the  comparison  of  Simons  correlation  against  experimental 
measurement.  It  is  observed  that  the  Simons  method  for  the  major  losses  in  a  parallel 
plate  heat  sink  yield  much  better  results  than  using  the  correlations  of  Equation  4-22  and 
Equation  4-23. 
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Figure  4.9  Comparison  of  the  Pressure  Loss  Coefficient  for  the  Simons  Method 


Figure  4.9  offers  insight  to  explaining  the  disagreement  with  other  predictive 
methods.  Since  all  three  methods  use  the  same  calculation  for  the  minor  losses,  and  the 
Simons  method  is  able  to  accurately  predict  the  measured  results,  then  the  disagreement 
must  be  due  to  the  major  loss  calculation  and  therefore  the  friction  factor.  This  further 
supports  the  conclusion  that  correlations  for  a  fully  developed  flow,  laminar  or  turbulent, 
cannot  be  used  to  describe  accurately  losses  for  a  hydrodynamically  developing  entry 
region. 

Another  observation  from  Figure  4.9  is  that  although  the  Simons  method, 
developed  only  for  hydrodynamically  developing  laminar  flow,  works  well  to  predict  the 
loss  coefficient  above  ReD  —  2,300.  In  Figure  4.9  the  curve  for  the  Simons  method 
prediction  was  extended  to  a  Reynolds  number  range  beyond  its  defined  applicability. 
This  was  done  to  illustrate  the  advantage  the  method  has  over  the  fully  developed 
turbulent  friction  factor  correlation  (Equation  4-23)  as  well  as  the  fully  developed  laminar 
correlation  (Equation  4-22). 
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CHAPTER  5 

CONCLUSION  AND  FUTURE  WORK  RECOMMENDATION 


5.1  Conclusion 

There  are  several  important  conclusions  from  this  thesis.  The  findings  from  the 
literature  review  highlighted  the  importance  of  fully  understanding  the  forced  convective 
heat  transfer  coefficient  at  high  altitudes.  This  thesis  proposed  a  method  to 
experimentally  measure  the  Nusselt  number  and  pressure  loss  coefficient  by  constructing 
a  wind  tunnel  to  operate  within  an  altitude  chamber. 

The  results  confirm  that,  for  altitudes  from  sea  level  to  52,000  feet,  the  Nusselt 
number  and  pressure  loss  coefficient  can  be  determined  from  experimental  measurement, 
numerical  simulation  and  suitable  correlations  obtained  at  sea  level.  Both  quantities  are 
functions  of  Reynolds  number  only.  Although  these  results  were  not  unexpected,  they 
now  serve  to  validate  predictive  methods  of  heat  transfer. 

In  addition  to  measuring  the  Nusselt  number  and  pressure  loss  coefficient  at  high 
altitude,  an  effort  was  made  to  evaluate  basic  methods  of  predicting  these  properties, 
through  correlation.  There  were  several  key  findings  from  this  investigation.  The  first 
was  the  entrance  effect  enhanced  heat  transfer  but  increased  pressure  drop  at  the  same 
time.  For  both  measurements  the  correlations  under-predicted  the  observed  data.  The 
primary  reason  was  due  to  the  fact  that  the  flow  within  the  heat  sink  channels  went 
through  a  thermally  and  hydrodynamically  developing  regime,  while  the  correlations 
assume  a  fully  developed  flow.  Additionally,  there  are  no  useful  correlations  that  can 
accurately  predict  the  Nusselt  number  in  the  transition  region  from  laminar  to  turbulent 
flow.  Within  the  range  of  Reynolds  numbers  tested  (500-14,000)  the  flow  exists  in  both 
laminar  and  turbulent  regimes.  This  is  of  interest  to  the  design  of  an  EMAS  thermal 
management  system  using  forced  convection,  as  it  would  likely  experience  this  transition 
within  the  flight  envelope.  Lastly,  the  traditional  laminar  flow  correlation  for  the  friction 
factor  under  predicts  the  pressure  loss  for  this  experiment  whereas  the  Simons  method  is 
more  accurate  because  of  the  entrance  effect  and  complex  geometry.  The  Simons  method 
is  also  more  accurate  in  predicting  the  pressure  drop  for  Reynolds  numbers  greater  than 
2,300.  However,  this  is  not  a  general  use  method  and  only  applies  parallel  plate  heat  sinks 
with  no  bypass  air. 

Finally,  it  was  clear  that  directing  the  air  stream  through  only  the  channels  of  the 
heat  sink,  eliminating  bypass  flow,  had  significant  higher  heat  transfer  over  an  external 
flow  situation.  Based  performance  difference  observed  between  the  data  and  flat  plate 
correlation,  a  recommendation  can  be  made  that  for  any  future  aircraft  TMS  use  a  no 
bypass  solution. 

5.2  Future  Work  Recommendation 

The  setup  chosen  for  this  experiment  was  meant  to  closely  emulate  an  aircraft 
EMAS  cooling  system.  Therefore,  any  pitfalls  uncovered  to  this  method  of  cooling  would 
directly  relate  to  problems  for  an  aircraft  system.  Additionally,  this  approach  provides  a 
good  catalyst  to  direct  future  work  within  the  field. 
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The  first  recommendation  to  come  from  this  research  would  be  the  development  of  a 
better  correlation  that  takes  into  account  the  entrance  effect.  For  any  TMS  of  this  type  the 
airflow  will  go  through  both  a  thennal  and  hydrodynamic  developing  state.  In  order  for  a 
correlation  to  be  useful,  in  predicting  performance,  it  needs  to  take  this  effect  into 
account. 

Another  recommendation  would  be  to  focus  research  on  the  development  of 
optimized  thermal  management  solution  for  the  range  of  Reynolds  numbers  across  the 
flight  envelope.  Further  research  can  be  conducted  to  verify  if  any  methods  for  heat 
transfer  surface  enhancement  are  equally  effective  across  the  flight  envelope. 
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APPENDIX  A 
FAN  SPEC  SHEET 


Rotron  /  Airscrew 


Tubeaxial  Fans  -  Propimax  2L  3715SF  MADE  IN  THE  USA 


Part  Number 

011858000 

Unit  Description 

PX2L.Q28MCDC  .S.  37 1 5SF 

Arflow  |CFM) 

110 

Pressure  Shutoff  (IWG) 

3.29 

Line  Voftage  (DC) 

28 

PhaseFreauency  (Hz) 

N'A 

Capacitor  ipF) 

N'A 

Power  (Watts) 

35.672 

Run  Cirrent  (Amps) 

1.30 

Locked  Rotor  (Amps) 

0 

RPM 

19.500 

Weight  (lbsikoi 

0.6/ 0.2 

DBA  at  3  feet 

75 

Max  Ambient  Temp  (*C) 

85 

Description 


Accessories: 


The  Propimax  2  brushless  DC  fan  is  an  axially  compact 
tubeaxial  fan  and  is  equipped  with  an  aerodynamic  venturi 
configured  case  Utilizing  the  unique  Rotron  inverted  motor 
construction  the  axial  length  of  the  fan  is  limited  to  the 
minimum  measurement  required  by  the  highly  efficient 
2-blade  propeller,  and  the  internal  or  external  drive 
electronics. 


Low  Speed  Warning  Detector 
Air  Flow  Switch 
Finger  Guard 


Construction 

The  venturi  case  and  propeller  are  cast  aluminum.  Rotron’s 
unique  inverted  motor  rotorprop  impeller  stainless  steel  shaft 
runs  on  two  high-precision  (ABEC  5),  double-shielded, 
stainless  steel  ball  bearings  with  an  integral  lubricant 
reservoir  to  give  longer  maintenance-free  life.  Brushless  DC 
drive  electronics  circuitry  is  mounted  in  an  external  module. 


Life  Expectancy 

Contact  AMETEK  Rotron  for  specific  ratings.  Rotron  rates  life 
expectancy  in  hours  L10  or  MTBF  at  a  given  ambient 
temperature.  This  is  the  time  at  which  90%,  for  L10,  and 
50%,  for  MTBF,  of  an  available  test  lot  are  still  operating. 


Environmental  Resistance  Guide 

Designed  to  meet  MIL-STD-810C 


Altitude 

Method  500.1,  pr  1 

Acceleration 

Method  513.2 

Vibration 

Method  514.2 

Shock 

Method  516.2 

Salt  Fog 

Method  509.1 

Rain 

Method  506.1 

Fungus 

Method  508.1 

Humidity 

Method  507.1,  prl.  II.  Ill 

Explos.  Atm. 

Method  511.1,  pr  1 

For  specifics,  please  contact  the  factory. 


Woodstock,  N.Y.  USA  (P)  845-679-2401  (F)  -  845-679-1870  (email)  milinquiry@ametek.com 
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a  [Measured 

1  B5IV.  1  iAp 

ROTRON  TECHNICAL  MOTOR  DIVISION  1 

Project:  iTest: 

TYP.  CURVE  1 A799-16E 

Model: 

PROPIMAX 

Type: 

2L 

NOTES: 

EST.  NEW  LAB  REF.  FAN 

E2900-7201.  THIS  TEST  26  VDC. 

P/N  011850000,  S/N:  108092223 

Date:  3 23-09 

Senes 

371 5SF 

N3P-20-1 

Type 

ECDC 

Flow:  LD  RD 

MOTOR 

Volts 

26 

s 

o 

_J 

u. 

LU 

LU 

cc 

LL. 

RPM 

19669 

Wave 

DC 

Line  Amps 

1.18 

Phase 

DC 

Locked  Amps 

— 

Hertz 

Watts 

MFD 
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APPENDIX  B 

WIND  TUNNEL  MINOR  LOSS  CALCULATIONS 


The  losses  from  changes  in  cross-sectional  area  of  the  duct  were  looked  at  first. 
This  instance  has  two  situations  that  result  in  a  cross-sectional  area  change;  increase  or 
decrease.  The  increase  in  cross-sectional  area  is  known  as  an  expansion  and  a  decrease  in 
cross-sectional  area  is  known  as  a  contraction.  The  minor  losses  for  these  are  computed 
differently  where  the  magnitude  of  loss  is  dependent  on  the  rate  (contraction  ratio:  N  & 
expansion/contraction  angles  a  and  fl)  at  which  the  change  develops.  If  the  rate  is  high 
enough  to  be  considered  an  almost  instantaneous  change  in  cross-sectional  area  then  the 
pressure  losses  associated  with  the  abrupt  change  would  be  at  a  maximum.  Wind  tunnel 
design  avoids  these  high  pressure  loss  situations.  However,  in  some  cases  this  could  not 
be  avoided,  therefore  it  was  necessary  to  calculate  their  contributing  loss.  But  first  it  was 
necessary  to  understand  the  other  contributing  factors  to  the  minor  losses. 

Equation  B  - 1  [19] 

(puj 

Vminor  loss  V  \  ^ 

For  each  case  where  a  minor  loss  occurred,  V minor  loss  (Equation  B  -  1)  was  used. 
This  value  is  dependent  on  the  flow  velocity  Um  and  the  loss  coefficient  rj.  In  the  instance 
of  an  abrupt  change  in  duct  cross-sectional  area  the  loss  coefficient  can  be  found  for 
either  a  sudden  enlargement  (Equation  B  -  2)  or  sudden  contraction  (Equation  B  -  3). 

Equation  B  -  2  [19] 


Equation  B  -  3  [19] 

For  these  two  equations  D  represents  the  diameter,  or  hydraulic  diameter  in  the 
case  of  a  non-circular  pipe.  The  subscripts  1  &  2  represent  the  diameter  of  the  pipe  before 
and  after  the  cross-sectional  area  change. 

For  situations  involving  gradual  change,  the  equation  for  estimating  the  loss 
coefficient  is  somewhat  more  complex.  For  the  cases  of  a  controlled  enlargement, 
Equation  B  -  4  was  used  and  for  a  controlled  contraction  Equation  B  -  5  was  used. 

Equation  B  -  4 

a  4 1  a/  A ±\2 

rj  =  4.0  xtan-x  Jtan-^1  J  +  rjf 
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Equation  B  -  5 


The  loss  coefficient  is  affected  by  the  contraction  ratio  or  expansion  angle  [19]. 
Figure  B.  1  shows  the  dimensions  of  a  controlled  expansion. 


Figure  B.  1  Wind  Tunnel  Diffuser  Section  [19] 


The  controlled  enlargement  or  diffuser  section  is  best  optimized  by  keeping  the 
expansion  angle,  a,  below  3.5°.  This  ensures  that  the  flow  does  not  detach  from  the  sides 
of  the  wind  tunnel  during  expansion.  However,  in  the  instances  where  this  feature  was 
used  (See  Figure  3.10:  Sections;  4-5  &  21-22)  there  was  not  enough  length  available  to 
follow  this  guideline.  The  effect  of  not  adhering  to  this  rule  is  minimal  since  neither 
section  takes  place  within  the  flow  conditioning  or  measuring  area  (See  Figure  3.10: 
Sections;  7-21).  Instead,  the  geometry  was  shaped  to  best  fill  the  available  space  and  used 
continuous  curves  to  smooth  the  transition.  The  pressure  loss  was  computed  using 
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Equation  B  -  4.  The  other  relevant  variables  needed  were  the  area  before  transition  A l5 
and  after  A2  as  well  as  the  additional  loss  factor  iff. 

Equation  B  -  6  shows  the  calculation  for  a  controlled  contraction. 

Equation  B  -  6 


0.02 

^  8  sin  a 1 2 


2i 


In  the  case  of  a  controlled  contraction,  this  experiment  used  three  (See  Figure  3.10: 
Sections;  1-2,  7-8  &  17-18),  two  were  in  the  conditioning  and  measurement  portions  of 
the  wind  tunnel.  Figure  B.  2  shows  the  concept  of  a  controlled  contraction  section  with 
angles  a  and  />'.  Of  the  two  sections  within  the  flow  conditioning  sections  only  one  could 
be  optimized  to  have  minimal  pressure  loss. 


Figure  B.  2  Wind  Tunnel  Contraction  Section  [19] 

This  was  the  section  immediately  after  the  fan,  Section  7-8.  Since  the  contraction 
ratio,  N,  was  already  fixed  the  contraction  angle  could  be  changed  by  altering  the  length 
over  which  the  contraction  takes  place. 

Once  the  design  of  these  sections  was  complete  it  was  time  to  move  on  to  the 
other  elements  of  loss  inside  the  wind  tunnel.  Three  corners  in  total  were  used  for  this 
wind  tunnel  (See  Figure  3.10,  Sections;  5-6,  11-12  &  13-14).  Several  solutions  were 
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considered.  Figure  B.  3  shows  types  of  ducting  bends  as  classified  by  American  Society 
of  Heating,  Refrigeration,  and  Air-Conditioning  Engineers,  Inc  [22], 


INEFFICIENT 

BRANCH  LOSS  COEFFICIENT-APPX.  1.20 
(HEIGHT/WIDTH  =  1.0) 


LACK  OF  TURNING  VANES  CAUSES 
EXCESSIVE  TURBULENCE  IN  FITTING: 
RESULTS  IN  VERY  HIGH  PRESSURE  DROP 


EFFICIENT 

BRANCH  LOSS  COEFFICIENT  -  APPX.  0.21 
(FVW=1 .0,  HEIGHT/WIDTH  =  1 .0) 


SMOOTH  RADIUS  ENSURES  UNIFORMITY  OF 
AIRFLOW,  GREATLY  REDUCING  TURBULENCE: 
RESULTS  IN  VERY  LOW  PRESSURE  DROP 


EFFICIENT 


BRANCH  LOSS  COEFFICIENT-APPX.  0.15 
(HEIGHT/WIDTH  =  1.0) 

SINGLE  THICKNESS 


TURNING  VANES  ENSURE 
UNIFORMITY  OF  AIRFLOW,  GREATLY 
REDUCING  TURBULENCE:  RESULTS 
IN  VERY  LOW  PRESSURE  DROP 


Figure  B.  3  (Left)  Expanding  Corner;  (Middle)  Constant  Diameter  Corner;  (Right) 
Expanding  Corner  with  Turning  Vanes  [22] 


Two  of  these  corners  result  in  a  comparatively  low  pressure  drop.  Both  the 
constant  diameter  corner  and  the  expanding  corner  with  turning  vanes  were  used  within 
the  wind  tunnel.  However,  the  expanding  comer  offers  a  better  method  for  minimizing 
the  size  of  the  bend,  thus  allowing  for  a  longer  straight  section  before  the  Pitot  tube 
measurement.  Figure  B.  4  shows  the  concept  design  of  an  expanding  comer  with  turning 
vanes. 


Figure  B.  4  Concept  Design,  Expanding  Corner  with  Turning  Vanes. 
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The  relevant  dimensions  are  the  comer  radius  r,  entrance  width  Went,  entrance 
height  Hent,  and  the  number  of  turning  vanes  N  as  well  as  the  size  of  the  vanes.  The 
entrance  height  and  width  were  identical  to  the  exit  height  and  width,  and  since  these 
were  dependent  on  the  duct  size,  will  both  be  equal  to  2.0”.  The  number  of  vanes  ( N )  for 
each  comer  were  found  using  the  following  relationship  shown  in  Equation  B  -  7. 

Equation  B  -  7 

d 

N  =  1.4  x- 
c 

It  was  important  to  have  the  fewest  number  of  vanes  needed  to  adequately  turn  the 
flow,  as  they  added  friction  losses.  The  amount  of  vanes  was  dependent  on  the  chord 
length  of  the  vanes  c,  and  the  length  of  the  diagonal  d.  The  chord  length  is  dependent  on 
the  design  of  the  individual  vanes,  including  their  camber.  Therefore,  it  was  possible  to 
select  the  number  of  vanes  based  on  selecting  a  comer  radius.  Figure  B.  5  shows  how 
these  dimensions  were  gathered. 


Figure  B.  5  Geometry  of  Guide  Vanes  to  a  90°  Expanding  Corner  [23] 


The  camber  for  the  corner  was  chosen  to  be  0.5  inches  as  this  represented  a 
quarter  of  the  duct  width.  From  this  perimeter  the  number  of  vanes  then  becomes  four, 
which  allows  for  the  calculation  of  the  pressure  loss  for  the  expanding  comer. 

Equation  B  -  8  shows  the  calculation  for  the  loss  coefficient  associated  with  an 

T 

expanding  corner,  where  r]M  is  an  additional  loss  based  on  the  ratio  - ,  which  for  0.25, 

went 

r]M  becomes  0.17  [19]. 

Equation  B  -  8 
rj  =  rjM  +  0.02  +  0.31  x— 

vvent 
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The  loss  coefficient  for  the  third,  constant  diameter  comer,  was  found  by  a  less 
involved  process.  Not  only  was  this  portion  of  the  wind  tunnel  before  the  fan  and  not 
within  the  test  section  it  also  had  the  largest  hydraulic  diameter  and  therefore  the  lowest 
velocity.  Because  of  this,  the  loss  within  the  corner  was  low  compared  to  other  parts  of 
the  duct.  Figure  B.  6  shows  the  relationship  between  the  size  of  the  elbow  and  the 
pressure  dropped  by  the  bend. 


1.00 


0.80  - 


R 

d 


Figure  B.  6  Constant  Diameter  Corner  Loss  Coefficient  Estimation  [24] 


The  multiple  curves  in  Figure  B.  6  are  for  varying  surface  roughness  for  the  walls 
inside  of  the  corner.  A  conservative  estimate  for  the  roughness  was  chosen  (0.001)  which 
for  an  R/d  ratio  of  1.3  gives  a  loss  coefficient  of  0.35. 

The  next  elements  of  the  wind  tunnel  that  needed  to  be  designed  and  their  pressure 
losses  computed  were  the  two  honeycomb  type  flow  conditioners  (See  Figure  3.10, 
Sections;  9-10  &  15-16).  The  purpose  of  these  features  was  to  eliminate  any  large-scale 
transverse  turbulence  in  the  flow  and  align  it  with  the  desired  direction  within  the  duct  for 
development  [25].  Each  of  the  two  flow  conditioners  served  a  necessary  purpose  within 
this  wind  tunnel  design.  The  first  was  to  eliminate  any  rotational  motion  or  ‘swirl’ 
imparted  onto  the  flow  from  the  fan  before  it  enters  the  first  corner.  The  second 
honeycomb  was  to  eliminate  transverse  turbulences  in  the  flow  after  having  traveled 
through  two  90°  comers.  It  also  helped  transition  it  into  a  fully  developed  state  necessary 
for  the  dynamic  pressure  measurement  by  the  Pitot  tube.  The  sizing  of  these  elements 
was  based  on  the  diameter  of  the  duct.  The  pressure  loss  of  the  honeycomb  was  found 
using  the  computation  for  a  sudden  contraction  (Equation  B  -  2)  followed  by  a  sudden 
enlargement  (Equation  B  -  3). 
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Figure  B.  7  shows  a  page  from  the  AMCA  publication,  Laboratory  Methods  of 
Testing  Fans  for  Aerodynamic  Perfonnance  Ratings  [15].  It  shows  the  guidelines  in 
selecting  the  length  and  cell  size  for  a  honeycomb  type  flow  straightener.  Thus  a 
honeycomb  size  that  fit  the  requirements  of  this  duct,  a  1.0  in.  long  honeycomb  with  1/8 
in.  cell  size  was  purchased. 


All  dimensions  shall  be  wrthin±0.005D  except  y 
which  shall  not  exceed  0.005D 


Figure  B.  7  Flow  Straighter  Sizing  [15] 
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APPENDIX  C 

PRESSURE  DROP  EXAMPLE 

O  "5 

Pressure  drop  computation  for  the  wind  tunnel  operating  at  0.325  nr/s,  1.2  kg/m 
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APPENDIX  D 
WIRING  DIAGRAM 

Experimental  Wiring  Diagram  for  Control.  Power,  and  Data  Aquisition 


vJ 
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Applied  Pressure  (PSIA) 


APPENDIX  E 

PRESSURE  SENSOR  CALIBRATION  INFO 
Sensor  Information 


Sensor  Part  Number: 
Serial  Number: 
Sensor  Range: 
Output: 
Supply  Voltage: 
Cal  Date: 


5221015PA2F2BE1S 
F033141 
0  to  15  PSIA 
0.05-5.0  VDC 
24  VDC 
6/16/2015 


Calibration  Curve  For  Absolute  Sensor  15PSIA 


Transducer  Output  (VDC) 


Regression  Statistics 

Multiple  R 

0.999999965 

R  Square 

0.999999931 

Adjusted  R  Square 

0.999999923 

Standard  Error 

0.001381822 

Observations 

11 

Intercept 

0.038272181 

X  Variable  1 

3.002888222 
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Applied  Pressure  (INWC) 


Sensor  Information 


Sensor  Part  Number: 
Serial  Number: 
Sensor  Range: 
Output: 
Supply  Voltage: 
Cal  Date: 


26412R5WD2DT1F 
8039682 
Oto  2.5  INWC 
0.05-5.05  VDC 
24  VDC 
3/12/2015 


Calibration  Curve  For  Sensor  Range  0-2.5  INWC 


Transducer  Output  (VDC) 


Regression  Statistics 

Multiple  R 

0.999999902 

R  Square 

0.999999804 

Adjusted  R 
Square 

0.999999782 

Standard  Error 

0.000389375 

Observations 

11 

Intercept 

-0.02455293 

X  Variable  1 

0.500055114 
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Applied  Pressure  (INWC) 


Sensor  Information 


Sensor  Part  Number: 
Serial  Number: 
Sensor  Range: 
Output: 
Supply  Voltage: 
Cal  Date: 


2641001WD2DT1F 
8039683 
Oto  1.0  INWC 
0.05-5.05  VDC 
24  VDC 
3/12/2015 


Calibration  Curve  For  Sensor  Range  0-1.0  INWC 


Transducer  Output  (VDC) 


Regression  Statistics 

Multiple  R 

0.999993514 

R  Square 

0.999987028 

Adjusted  R 
Square 

0.999985586 

Standard  Error 

0.001255886 

Observations 

11 

Intercept 

-0.01045236 

X  Variable  1 

0.199952322 
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Applied  Pressure  (INWC) 


Sensor  Information 


Sensor  Part  Number: 
Serial  Number: 
Sensor  Range: 
Output: 
Supply  Voltage: 
Cal  Date: 


2641R25WD2DT1F 
8039684 
0  to  0.25  INWC 
0.05-5.05  VDC 
24  VDC 
3/12/2015 


Calibration  Curve  For  Sensor  Range  0-0.25  INWC 


Transducer  Output  (VDC) 


Regression  Statistics 

Multiple  R 

0.999996214 

R  Square 

0.999992428 

Adjusted  R  Square 

0.999991586 

Standard  Error 

0.000240106 

Observations 

11 

Intercept 

-0.00214206 

X  Variable  1 

0.05001017 

84 

Approved  for  public  release;  distribution  unlimited. 


Applied  Pressure  (INWC) 


Sensor  Information 


Sensor  Part  Number: 
Serial  Number: 
Sensor  Range: 
Output: 
Supply  Voltage: 

Cal  Date: 


26410  R1WD2DT1F 
8039685 
0  to  0.10  INWC 
0.05-5.05  VDC 
24  VDC 
3/12/2015 


Calibration  Curve  For  Sensor  Range  0-0.10  INWC 


Transducer  Output  (VDC) 


Regression  Statistics 

Multiple  R 

0.999995399 

R  Square 

0.999990797 

Adjusted  R  Square 

0.999989775 

Standard  Error 

0.000105692 

Observations 

11 

Coefficients 

Intercept 

-0.00087114 

X  Variable  1 

0.019971059 
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APPENDIX  F 

THERMOCOUPLE  CALIBRATION  PROCEDURE 


A  total  of  16  thermocouples  were  utilized  for  this  experiment.  Their  calibrations 
were  conducted  prior  to  their  instillation  within  the  experimental  setup.  The  process  for 
their  calibration  was  to  use  a  calibrated  reference  temperature  RTD  [Hart  Scientific 
1502A  with  Hart  Scientific  5628  probe  (Standard  Error,  ±0.033°C)]  and  a  silicone  oil 
bath  (Hart  Scientific  7320)  to  control  the  temperature.  The  range  of  temperatures  chosen 
for  calibration  was  0  -  150°C.  To  carry  out  the  process  of  calibration  a  Lab  VIEW  VI  was 
created  to  control  the  bath  temperature  and  record  thermocouple  measurements  via  the  NI 
cRIO  9112.  Four  NI  9211  thermocouple  modules,  each  having  four  inputs,  were  used 
with  the  cRIO. 

The  process  for  calibrating  the  thermocouples  was  to  secure  all  of  the  probe  tips 
close  to  the  RTD  tip  and  submerge  them  in  the  fluid.  This  was  to  ensure  that  the 
temperature  at  each  thermocouple  tip  was  as  close  to  the  RTD  reference  temperature  as 
possible.  The  probe  bundle  was  then  positioned  securely  in  place  within  the  oil  and 
insulation  was  placed  over  the  opening.  Figure  F.  1  shows  the  calibration  in  progress. 


Figure  F.  1  Thermocouple  Test  Setup 


Once  the  setup  was  completed  the  Lab  VIEW  program  initiated  the  calibration  by 
commanding  the  bath  to  reach  and  hold  a  steady  temperature.  With  a  range  of  0  -  150  °C, 
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measurements  were  taken  at  every  10  °C  increment  starting  at  0  °C  climbing  to  150  °C 
then  back  down  to  0  °C.  A  total  of  10  readings  were  taken  at  one  sample  per  second,  for 
each  temperature  increment  and  averaged.  The  reading  for  each  temperature  point  was 
taken  when  the  standard  deviation  of  the  RTD  fell  under  0.0005  °C.  Once  recording  is 
finished  the  VI  commands  the  bath  to  the  next  set  point  and  the  process  repeats.  A  total  of 
29  averaged  temperature  measurements  were  taken. 

It  is  worth  noting  that  at  the  time  of  calibration  the  bath  refrigeration  system  was 
not  working  properly.  For  temperature  points  below  ambient  temperature  an  ice  bath  was 
prepared  and  used  instead  of  the  oil  bath.  Also  worth  noting,  after  the  initial  calibration 
four  thermocouples  were  damaged  during  installation  and  required  a  second  calibration 
after  they  were  repaired.  The  results  of  the  calibration  are  as  follows, 


Recorded  Thermocouple  Temperature  VS.  RTD  Temperature 


y 
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Figure  F.  2  Thermocouple  Calibration  Results 


Figure  F.  2  shows  all  of  the  thermocouple  measurements  taken  for  the  calibration. 
It  is  apparent  that  the  lack  of  active  cooling  caused  a  group  of  measurements  to  be  taken 
slightly  above  room  temperature.  However,  this  does  not  adversely  affect  the  calibration. 

Once  all  the  data  was  recorded  the  next  step  was  to  create  a  regression  equation 
for  each  thermocouple.  For  this  a  linear  best  fit  regression  equation  was  chosen.  Table  F  - 
1  shows  the  results  from  the  linear  regression  analysis  for  the  thermocouple  temperature 
vs.  the  RTD  temperature. 
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Table  F  - 1  Thermocouple  Regression  Analysis 


THERMO¬ 

COUPLE 

SLOPE  (°C/°C) 

OFFSET  (°C) 

STANDARD  ERROR 
(±°C) 

TC01 

0.9945 

-1.35 

0.4165 

TC02 

0.9947 

-1.32 

0.4205 

TC03 

0.9990 

-1.55 

0.4083 

TC04 

0.9972 

-1.34 

0.3807 

TC05 

0.9986 

-1.55 

0.4465 

TC06 

0.9987 

-1.35 

0.3982 

TC07 

0.9987 

-1.16 

0.3849 

TC08 

0.9984 

-1.00 

0.3790 

TC09 

0.9986 

-1.33 

0.4239 

TC10 

0.9988 

-1.21 

0.4047 

TC11 

0.9945 

-1.13 

0.4175 

TCI  2 

0.9941 

-1.04 

0.4172 

TCI  3 

0.9954 

-0.97 

0.3947 

TC14 

0.9980 

-1.41 

0.4391 

TCI  5 

0.9978 

-1.37 

0.4705 

TCI  6 

0.9945 

-0.97 

0.4465 

The  slope  and  offset  were  programed  into  the  experimental  Lab  VIEW  VI  to  make 
the  temperature  correction.  The  Standard  Error  was  later  used  in  the  uncertainty  analysis. 
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APPENDIX  G 

HEAT  SINK  TEMPERATURE  SETTLING  TIME 


The  following  figures  were  produced  to  investigate  the  settling  time  for  steady 
state  conditions  of  the  air  velocity  and  heat  sink  temperature,  under  varying  conditions. 
The  dotted  blue  line  denotes  when  data  recording  was  triggered. 
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The  following  tests  were  run  at  extended  time  to  validate  the  ability  of  the 
program  to  determine  steady  state.  The  dotted  blue  line  denotes  the  time  at  which  normal 
triggering,  to  log  data,  would  have  occurred. 
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APPENDIX  H 

RESULTS  UNCERTAINTY  ANALYSIS 


The  uncertainty  of  the  results  of  this  thesis  was  detennined  by  the  method  of  error 
propagation  through  partial  derivatives.  The  standard  error  of  a  relationship,  f(x,y,z ),  is 
represented  by  the  sum  of  the  squares  of  the  product  of  its  partial 


derivatives 


df(x,y,z )  df(x,y,z )  „  df(x,y,z ) 


and  their  individual  measurement  uncertainties. 


7f(x,y,z ) 


df(x,  y,  z) 
dx 


df(x,  y.zy 


In  the  above  equation,  ox,  oy  &  oz  are  the  individual  uncertainties  for  the 
respective  variables  x,y,8iz.  They  are  detennined  from  measurement  uncertainty  or 
propagation  of  uncertainty  through  partial  derivation. 

The  measurement  uncertainty  was  determined  by  the  addition  of  the  calibration 
uncertainty  and  the  standard  error  of  the  regression.  The  calibration  uncertainty  was  taken 
as  the  uncertainty  of  the  device  used  for  calibration.  The  standard  error  is  the  product  of 
the  curve  fitting  the  calibration  results  to  the  best  fit  equation. 

The  regression  standard  error  is  found  by  the  method  of  Least  Squares. 

a 

w  here  £t  is  the  error  for  each  individual  measurement  and  N  was  the  total  number  of 
measurements.  The  results  for  all  the  sensors  used  in  this  experiment  are  shown  in  the 
table  below. 
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Table  H  - 1  Measurement  Uncertainty 


MEASUREMENT 

UNITS 

CALIBRATED 

RANGE 

UNCERTAINTY 

NOMENCLATURE 

STANDARD 
ERROR  (±) 

CALIBRATION 

UNCERTAINTY 

(±) 

MEASUREMENT 

UNCERTAINTY 

(±) 

SHUNT  CURRENT 

Amps 

0-5 

^l_Shunt 

0.000023 

0.0031 

0.0031 

HEATER  RESISTANCE 

Ohms 

N/A 

aR_Heater 

0.000023 

N/A 

0.000023 

DYNAMIC  PRESSURE 

lnH20 

0-0.10 

aP_Dyc:0.10 

0.00004 

0.00011 

0.00015 

DYNAMIC  PRESSURE 

lnH20 

0-0.25 

aP_Dyc:0.25 

0.00010 

0.00024 

0.00034 

STATIC  PRESSURE  DROP 

lnH20 

0-1.0 

aP_Stc:1.0 

0.00040 

0.00126 

0.00166 

STATIC  PRESSURE  DROP 

lnH20 

0-2.5 

aP„Stc:2.5 

0.00100 

0.00039 

0.00139 

ABSOLUTE  PRESSURE 

PSI 

0-15 

aP_ABS:lS 

0.00600 

0.00138 

0.00738 

AIR  AT  INLET  1 

°C 

0-150 

aTC01 

0.033 

0.4165 

0.4495 

AIR  AT  INLET  2 

°C 

0-150 

°TC02 

0.033 

0.4205 

0.4535 

HEATSINK  BASE  1 

°C 

0-150 

°TC03 

0.033 

0.4083 

0.4413 

HEATSINK  BASE  2 

°C 

0-150 

°TC04 

0.033 

0.3807 

0.4137 

HEATSINK  BASE  3 

°C 

0-150 

°TC05 

0.033 

0.4465 

0.4795 

HEATSINK  BASE  4 

°C 

0-150 

aTC06 

0.033 

0.3982 

0.4312 

HEATSINK  BASE  5 

°C 

0-150 

aTC07 

0.033 

0.3849 

0.4179 

HEAT  SINK  TIP  1 

°C 

0-150 

aTC08 

0.033 

0.3790 

0.4120 

HEAT  SINK  TIP  2 

°C 

0-150 

°TC09 

0.033 

0.4239 

0.4569 

HEAT  SINK  TIP  3 

°C 

0-150 

aTC10 

0.033 

0.4047 

0.4377 

AIR  AT  EXIT  1 

°C 

0-150 

°TC11 

0.033 

0.4175 

0.4505 

AIR  AT  EXIT  2 

°C 

0-150 

°TC12 

0.033 

0.4172 

0.4502 

HOUSING  EXTERNAL  1 

°C 

0-150 

°TC13 

0.033 

0.3947 

0.4277 

FAN 

°C 

0-150 

aTC14 

0.033 

0.4391 

0.4721 

HOUSING  EXTERNAL  2 

°C 

0-150 

aTC15 

0.033 

0.4705 

0.5035 

AIR  AT  EXIT  3 

°C 

0-150 

aTC16 

0.033 

0.4465 

0.4795 

Other  measurements,  such  as  physical  dimensions  were  also  used  in  the  derivation  of 
the  experimental  uncertainty.  Since  calibration  was  not  necessary  for  physical  dimensions 
a  standard  error  was  used  based  on  the  general  tolerance  for  each  specific  part  (±0.02  in. 
for  manufactured  parts  and  ±0.035  in.  for  3D  printed  parts) 


•  Heat  Transfer  Coefficient 

Within  this  thesis  the  HTC  was  used  to  measure  the  performance  of  the  heat  sink.  The 
uncertainty  of  this  measure  is  a  combination  of  many  various  measurements  of 
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temperature,  power,  and  pressure.  The  following  is  an  accounting  of  the  total  uncertainty 
for  the  HTC.  To  begin  the  equation  for  the  measured  HTC  of  the  base  of  the  heat  sink  is 
as  follows, 


h 


base 


Pt.HS-^base 


The  heat  transfer  coefficient  in  its  non-dimensional  fonn  is  the  Nusselt  Number,  Nuavg . 


NuaVg  hbase  ^ 


■fluid 


Where  L*  is  the  characteristic  length  that  is  computed  by, 

L*=jAt 

Using  substitution, 


ibase 


Nu 


1  vi 


base 


avg 


Rt,HS Abase  ^ fluid  Rt,HSkfiuidAbas^ 


The  thennal  resistance,  Rt,HS’  is  found  with  the  following  equation. 

n  _  A Tbase 
Rt,HS  ~  "7j 

VConv 

So,  using  substitution  again,  the  equation  for  the  Nusselt  number  becomes. 

Qconv 


Nu  avg  ~ 


NTbasekfluidAbase2 


Form  here  the  three  uncertainties  contributing  to  this  measurement  are:  Pheater> 
ATbase  &  kf  luid.  The  following  derivation  was  used  to  determine  the  total  uncertainty 
propagation  for  the  Nusselt  number  measurement. 


dNuaVg 
9  Qconv 


=  AT, 


-ii 


base 


" fluid 


A 


base 


1 

2 


dNuaVg 

dA  Tbase 

d  NllaVg 

9k fluid 
dNuaVg 

9Abase 


QconvNTbase  kfiuid  Abase 


QconvNTbase  kfiuid  Abase 


'  n  QconvNTbase  kfiuid  A 


■fluid  nbase 


3 

2 
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2 


Qconv  dPh 


;)24 


'  base  dA Th 


*  fluid  dk. 


l^]+(c 

dd )  V 


hase  dAbase. 


°Nuavg2  Qconv  {j^base  k fluid  Abase  2)j 

T  Tbase  ~QconvATbase  k.fiui(i  Abase 
”1”  i^k fluid  (  Qconv^base  k.fiui(i  A base 

”1”  (^Abase  2  Qconv^Tbase  Abase  2^ 


The  uncertainties,  0QConv,  (7&rbase  akfluid  &  GAbase  needed  to  be  individually 
detennined.  The  temperature  difference  of  the  base  to  the  air,  A Tbase,  was  found  through 
the  following. 

ArT,  .  trco3+rco4+rco5+rco6+rco7)  (rcoi+rco2) 

A Tbase  =  ( Tbase  ~  ^ air-in )  =  “ - - - “  ~  " - “ - " 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

dAThnQO  1  dAThncp  1 


dA Tbase  _ 

1 

dA  Tbase 

dTCOl  “ 

~2’ 

dTC02 

dA  Tbase 

1 

dATbase 

dTC03 

~5: 

dTC  04 

dATbase 

1 

dATbase 

dTCOS 

~  5; 

dTC06 

dTC 07  5 

2  _  (  dAThasA2  /  dATbase\2  (  dA Tbase 

°^Tbase  -  ^rcoi  dTC01  )  +  [aTC02  dTCQ2  J  +  [aTC03  gTC03 

,  l  dA Tbase\2  (  dA TbasA2  /  dATbase\2 

+  [^TCOA  QTCQ4  )  +  [aTC05  dTCQS  )  +  \°TCOe  QTCQb  ) 

(  dA Those)2 

+  \0TCO7  gTCQ7  J 


(  dATbase' 

+  \aTC07  gTCQ7  ^ 


°TCOl[-^ 
+  i^TCOS 


(“2))  +  (ffTC03  (5))  ' 

°rC05  (5)^  +  (^aTC06  (5)^  +  (^aTC07 


+  (  aTC04  (  r 
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For  the  uncertainty  of  the  power  being  dissipated  through  convection  by  the  heat 
sink,  Qconvn  the  following  equation  was  used. 

Qconv  ~  ( J Shunt  Re, Heater)  R losses 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

9  Qconv 


9 1  Shunt 
d  Qconv 

dRe, Heater 
d  Q  Conv 


2 / Shunt  x  Re, Heater 


=  l 


Shunt 


dQ 


=  -1 


losses 


® Qconv  l  A 


dS^\  +  ( a„ 


Shunt  I 

u 1  Shunt  ' 


dR 


dQconv  Y 


e, Heater , 


Q  losses 


9  Qconv  \ 
dQi  osses ' 


° Qconv2  -  (aIShunt(21  Shunt  X  Re, Heater))  +  (aRe, Heater  0  Shunt2))  +  (CTQlosses  (  !)) 


Of  these  the  uncertainties,  &ishunt  &  cr/?e  Heater  are  currently  known  from 
calibration.  However,  the  uncertainty  of  the  losses  due  to  conduction,  QCiiosses,  and 
radiation  Qr, losses  will  need  to  be  found. 

Qiosses ,  is  the  sum  of  the  two  following  variables. 

Qlosses  —  Qc.losses  T  Qr,  losses 

Therefore,  an,  becomes, 

7  v  losses  7 

® Qlosses  (^Qc.losses  y +(<  TQr,  losses  ) 

The  losses  through  conduction  are  found  by  the  following  equation. 

„  _  n  hi  side  _HS  hT  Bottom  HS 

Qc, losses  ~  ^  p  '  7J 

^t, Housing  _Side  rl  t,Housing_Bottom 

Breaking  this  down,  the  uncertainty  of  this  measurement  was  determined  from  the 
following  derivation. 

d Qc.losses  _  2  d Qc.losses  _  OAT,  n  -2 

a  » rp  —  p  >  ci  p  —  —£l±L side_HS^t,Housing_Side 

Otsl  side_HS  rtf.  Housing  _Side  ®  ^t,Housing_Side 

d  Qc.losses  _  1 

dATBottorilHS  Rt, Housing  Bottom 


96 

Approved  for  public  release;  distribution  unlimited. 


dQ 


c,  losses 


dR 


=  -AT, 


Bottom_HS^t, Housing  _Bottom 


-2 


t, Housing  _Bottom 


(J, 


2  _  (  d  Qc, losses  \  (  d  Qc, losses 

Qc, losses  1  SideJiS  SAT  /  l  Housing  Side  HP 

\  °  ***  Side_HS  /  V  ° K  t, Housing  _Side  / 

O 


*Side 

d  Qc, losses 


9  Qc, losses 


I  v-  1  I  V.  LjlUJJCJ 

I  Bottom  HS  HAT1  I  \  Rt, Housing  Bottom  HD 

\  " Bottom_HS /  \  ‘  ” t't,Housing_Bottor 


O, 


Qc, losses 


_Side 


-2 


(°A TsideHS(2Rt, Housing 

+  {°Rt, Housing  side  {~^-^Side_HS^t, Housing _Side 
(^A T Bottom_HS  ( Rt,Housing_Bottom  )) 

(° Rt,H0using_Bottom  (_ ^B ottom_HS ^t,Housing_B ottom 


-2 


Of  these  uncertainties,  aR 
needed  to  be  detennined. 


t,Housing_Side 


'  ^AT side_HSr  Housing _Bottom  Bottom_HS 


The  calculation  for  Rt,Housing_side  is  as  follows, 

Aside 


R 


t, Housing  _Side 


+ 


K PC^HS_Side  h c,Side^HS_Side 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

dRt 


vt,Housing_Side  _  -1a  -1 

—  Kpr  Al 


9Lc; 


'-PC  siHS  Side 


Side 


dR 


t, Housing _Side  _  -2  .  -1 

—  —^sideRpc  AHs_side 


dK, 


PC 


dR 


t, Housing  side 


dA 


HSSide 

dR 


—  A HSside  2{jJSide^PC  1  +  h c.Side  ) 


t, Housing _Side  _  ,  -2  .  -1 

—  ~llc,Side  Afisside 


dhc,Side 


dRt 


°Rt, 


2  _  1  _  “ -'t,Housing_Side  \  ,  /  _  u nt, Housing _Side 


dRt 


=  °L 


,Housing_Side  y  LSide 

(  dR 


+ 


d  Aside 

\2 

t, Housing _Side  \ 


+  °K 


PC 


dK, 


PC 


HS_Side  Q ^ 


+  <Th 


dR 


HS  Side 


c.Side 


t,Housing_Side 

dh-c.side  / 
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aH 


,Housing_Side 


+ 

+ 

+ 


iHS_Side 
-2 


A) 


(' aLside{KPC  lA 
(^KPC  (  Lside^pc  ^  AHS_Side 
( °AHS_Side  ~AHSside  2  iLSide  KpC 

{^C,Side(.-h 


lc,Side 


A 


HS_Side 


1))‘ 


+  h, 


-i 


c.Side 


The  uncertainty  for  the  area  of  the  heat  sink  side,  AHSSide,  needed  to  be  found. 
Using  a  standard  manufacturing  tolerance  of  ±  0.02  in.  (  +0.00005  m),  the  uncertainty 
aA  side  was  found  through  the  following  derivation. 

2  2 
2  _  ( _  9Ahs  side\  ,  ( _  dAHS jide\ 

aAHS  Side  -  \ Vhhs  9Hhs  )  +  ^u/HS  dw^  j 

c tHhs  =  °wHS  =  0.000508m 

aAHS_side2  =  (0-00005 (kTHS))2+ (0.00005 (Hhs)) 2  =  2.92  x  10-9m2 
This  is  also  the  same  method  used  to  find  the  uncertainty  of  the  Abase  ( AHS  Bottom ) 

<*Abase 2  =  (0.00005  (kk^))  2  +  (0.00005  (SHS))2  =  3.21  x  10-9m2 


The  calculation  for  A Tside  HS  was  made  using  the  following  equation. 


ATside_HS  —  A(T aVg  side  TC 13)  — 


(  f  A7’C04+ATC06'\  /ArC09+ATC10\\ 


It 


-J  +  V 


-  ATC13 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


dATside _HS  _  1 

dTCOA  ~  4' 

dATside  HS  _  1 

3TC09  ~  4' 

9ATside  HS  _ 

3TC  13 


dATside _HS  _  1 

arco6  ~~  4 

dATside  HS  _  1 

arano  ~~  4 


/  dATside  hs\2  (  9ATside  HS^2  /  dA Tside  HS^ 

V  TC04  3TCQA  )  +  V  rc°6  3TC06  )  +  V  rC09  3TC09  ) 
(  dATside  Hs\2  (  dATside  HS^2 

+  vrC10  arcio  )  +  vTrC13  3tc  13  / 
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aA  TsideHS2  -  (  °VC04  (4)  j  +  (  CTrcoe  (4)  j  +  ( aTC09  (4)  j  +  ( °TC10 


+  (otci3(.  1)) 


To  estimate  the  thermal  resistance  of  the  bottom  of  the  housing,  Rt, Housing -Bottoms 
the  following  equation  was  used. 

p  _  n  Lpp  L Bottom  0  1 

^t, Housing -Bottom  ~  ^  is  a  'is  a  '  ^  u  A 

^FF^HS -Bottom  ^PC^HS -Bottom  ^c.Bottom^HS -Bottom 


Next,  the  uncertainty  becomes, 

dR 


t, Housing_Bottom  _  -1  A  -1 

—  ZKff  Ahs  Bottom 


dL 


FF 


dR 


t, Housing  Bottom  _ .  ,  -1 

- =  -2 LFFKFF  zAhs  Bottom 


dK, 


FF 


dR 


t, Housing -Bottom 


dA 


HS_Bottom 


—  A HSside  2 (2 LppKpp  1  +  LBottomKpC  1  +  3hCj Bottom  ) 


9  Rt, Housing -Bottom  _  T/r  -1  A  -1 

—  KPC  Ahs  Bottom 


dL 


Bottom 


dR 


t, Housing -Bottom  _  T  ly.  -2  A  -1 

—  f BottomRpc  Apis  Bottom 


dK ] 


PC 


dR 


t, Housing -Bottom  _  ou  -2  A  -1 
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u  ‘‘■c, Bottom 


alh, 


Housing  _Bottom 


dR 


+  [<Ja 


t , Housing  -Bottom 


dLpp 
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t, Housing  _Bottom 
HS_Bottom  Q 


+  Ok 


dR 


1HS_Bottom 
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°Kff 
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t, Housing -Bottom^ 


dKpp 

dR 

L Bottom 


t, Housing  -Bottom^ 


dL 


t, Housing  -Bottom 


PC 


dK 


PC 
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Bottom 

d  Rt, Housing  _Bottom\ 

c,  Bottom  HU  I 

u  nc, Bottom  / 
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CJRt, 


Housing  _Bottom 


—  (°Lff(2^fF  1  AHs_Bottom  *))  +  iyaKFF{~ ^tFFKFF  2 AHS  Bottom  x)) 
+  'Ahs  Bottom  (_ AHS_Bottom  2  (2 LppKpp  1  +  LBottomKpC 


+  3h, 


lHS_Bottom 
-1 

c, Bottom 


(c 


+  (Oi. 


L Bottom 


T"  (^Kpc[  LrtnttrnnKpr  A 

3”  Bottom  (  3hc, Bottom  AHS  Side 


[Kpc  1  Apis  Bottom  X)) 
2 


■‘Bottoml'-PC  rlHS_Bottom 
-2 


-or 


The  calculation  for  the  difference  in  temperature  of  the  heat  sink  base  and  outside 
the  housing  bottom  is  as  follows. 


AT, 


Bottom  HS 


avg.Base 


TC 15) 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  equation. 

2  _  (  dATBottom_Hs\  (  9 ATBottom  HS\ 

TBottomHS  ~  \^TavgiBase  j  +  \aTC15  (IT CIS  ) 

Since  <JTavg  Base  has  been  previously  calculated,  the  equation  is  simplified. 

Bottom_HS  avg.Base  3~  (^Vci5( 


The  next  step  was  to  compute  the  uncertainty  of  the  losses  from  radiation.  The 
calculation  for  Qr  ioss  is  shown  in  the  following  equation. 

Qr.loss  =  Nf_chQl-2  +  Adas  (Ts^  —  Tamb 3") 

For  simplicity  this  equation  was  broken  into  two  parts;  the  front  and  rear  profile  radiation 
(Qr  p)  and  radiation  exchange  between  the  surface  of  the  heat  sink  channels  and  PC 
housing  (<3r,i-2)- 


Qr,loss  —  Qr, 1—2  3"  Qr,P 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

2  _  l  dQril0SS\2  (  dQr,l0SS\2 

°Qr,loss  \°Qr, 1-2  QQrl_ 2)  y  dQrp  J 

^Qr.loss2  =  (aQr, 1-2  W)  +  (^r.pC1)) 

The  first  uncertainty  to  be  investigated  was  the  profile  radiation,  aQrp .  The  equation  for 
Qr  p  is  as  follows. 
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Qr,p  ~  2i4pCTf(Ts^  T'amb^') 


The  variables  in  this  equation  are  as  follows. 

•  Ap  =  Surface  area  of  heat  sink  profile 

•  s  =  Emissivity  of  unfinished  aluminum  -0.055  (0.04-0.07) 

•  g  =  Stefan-Boltzmann  Constant  (no  uncertainty) 

•  Ts  =  Average  heat  sink  surface  temperature 

•  Tamb  =  Average  of  air  at  inlet  temperature. 

The  uncertainty  of  Qr  p  was  detennined  from  the  following  derivation. 

9Qr 


cr,p 


dAp 

dqQr, 


=  2 os(T*  -  Tamb*y, 


v-=2Avo(J^-Tam^y, 


de  -  p 

dQr.p 


—  2  A-nGS 


d(T*-Tamb*)  ~  ^ 


2  _  (  Wr.p' 

aQr,p  -  (  % 

p 


^  dA 


,  (  dQr,p\  ,  (  9Qr,p  ' 

+  r*  —  +  aoS-T^*), 


O. 


Qr,p 


=  (%(2 -  Tamb 4}))2  +  [oa  {2Ape(Ts^  -  Tamb*))f  +  (o£  {2Apo(Ts^ 
—  7’amfe^))  j  +  (a(Ts4-Tamb4)(2ApG£yj 


From  this  equation,  two  uncertainties  need  to  be  further  described,  oA  &  &(TS4-Tamb4y 

First,  Ap  was  found  by  a  direct  measurement  with  the  model.  Therefore,  it  was 
understood  as  having  the  same  general  tolerance  applied  to  other  manufactured  parts, 
±0.02  in.  Its  total  uncertainty  was  approximated  by  the  following. 


°Av  gbhs 


dAr 


dH 


+  Op 


HS  > 


My  V 

dPweJ 


oBu.  —  °p  r  —  0.000508m 

DHS  Mvet 

oAp2  =  (0.00005(BHS))2  +  (0.00005(Pwet))2  =  2.25  x  10"8m2 
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Secondly,  the  uncertainty  0(Ts4-Tamb*)  is  found  from  the  measurements  for  Ts  and 
Tamb,  which  are  the  average  temperatures  of  the  heat  sink  surface  and  air  at  inlet,  in  K. 

(  (d(Ts4  -  Tamb4)\\2  _  (  ( d(T4-Tamb4)\\ 2 

a(Ts4-Tamh4)  =  o-rs  - ~ - II  +  I  aTn 


dTQ 


■  amb 


%S4-U4)  =  (^rs(47s3))  +  (crr  amb  {—4Tamb  3)  ) 


dTamb 
2 


Next  the  uncertainties,  <jTs8i  crTamb  ■>  needed  to  be  individually  detennined. 

The  equation  for  Ts  is  as  follows. 

1  ((TC03  +  TC04  +  TCOS  +  TC06  +  TC07 )  (TC09  +  TCIO)^ 

Ts=2\  5  +  2  j 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


dTQ 


dTQ 


dTQ 


dTQ 


dTC 09  4'  dTC 10  4'  dTC 03  10'  3TC04  10' 


an 


an 


an 


arcos  io'  arco6  10'  arc 07  10 


nr,2  -  I  n’rcos 


an 


arco3 


)  +  (aTco4  dTC04)  +  {aTC05  dTC  os)  +  {aTC06  dTeoe) 


an 


+  ( °rco7 


V 


dTC07j 


an 


+  ( °TC09 


V 


arco9/ 


an 


+  l  °rcio 


V 


arcio/ 


°T  \  = 


(lo))  +(<7r“4(io))  +  (io))  +('7™6(io)) 

+  (<r,'C07  (lo))  +  (ffTC09  (4))  +  (aTC10  (j)) 


The  equation  for  Tamb  is  as  follows. 


Tamb 


(TC  01  +  TC02) 


The  uncertainty  of  this  measurement  was  determined  from  the  following  derivation. 

2  _  (  dTamb\ 2  /  dTamb\2 

aTamb  -  [aTC01  dTC01)  +  \°TC02  QTCQ2) 


(7^ 


amb 


—  aTC01 


+  °TC02 
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Lastly,  the  uncertanity  for  the  second  term  of  the  radiation  losses  was  found 
through  the  following.  The  equation  for  estimating  the  radiation  exchange  between  the 
heat  sink  channels  and  the  PC  housing  is  shown  below. 

Qr,  1-2  —  Nf_chQl-2 

Where  the  total  radiation  loss  (Qr  1-2)  is  a  product  of  the  loss  for  a  single  channel  ( Q1-2 ) 
multiplied  by  the  total  number  of  channels,  Nf  _ch.  Nf_ch  is  a  constant  value  representing 
the  total  number  of  heat  sink  fin  channels  (5)  and  has  no  uncertainty.  Therefore,  the 
uncertainty  for  the  total  radiation  loss  from  all  the  heat  sink  channels  simplifies  to  the 
following  equation. 

^Ws-s2  = 

®QrS-S  ~  ^<7(3l-2 

The  equation  for  Qi_2  is  as  follows. 

„  _  <Kr14-r24) 

Ql~2  l-e.  1  l-s2 
^1^1-2  £2^2 

It  was  deirseable  to  break  this  equation  up  into  two  tenns;  the  temperature  difference 
between  the  two  surfaces  (T^  —  r24),  and  the  surface  emissive  resistance,  Rs.  The 
simplified  equation  for  Q i_2  then  becomes, 


Q1-2  — 


<j(t^  -  r24) 


Rc 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


o, 


Q 1 


9Ql-2 

-2  -  {  a(Ti4-T2*)  _  724) 

2 


H^) 


aQi -22  =  (°(TS-TMaRs  '))  +(^s(-<7’14-T24)«s  2)) 


Next  the  uncertanities  <7(Ti4_r24-)  and  aRs  needed  to  be  computed.  The  uncertainty 
(jp1  4^  is  found  from  the  following  derivation. 


cr(Ti4-r24)  -  ^ 


'd{7^  -  T24y 
,  3Ti  ' 


+  Ot 


'a(r!4-r24)' 


dT7 


cr(r14-r24)  =  (^C4^3))  +  (^C-4^3)) 


Since  T±  is  the  average  heat  sink  surface  temperature,  its  uncertainty  (er7i )  has  already 
been  detennined,  oTi  —  oTs. The  calculation  for  the  uncertainty  of  the  second  surface,  T2, 
still  needed  to  be  computed.  T2 was  found  from  the  following  equation. 
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_  1  f(TC01  +  TC02)  (TC11  +  TC12  +  TC16)\ 

T2=2\  2  +  3 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

°T^  =  (ffTC01d¥an)  +  (aTC02^k)  +  (CTrC11  a?Hl)  +  (aTC12d¥ Si)  +(aTC16d¥^) 
-  '  ^ 


-  Urcoi  (4)  j  +  UtC02  (4)  )  +  [°TC11 


The  surface  emissive  resistance  was  found  through  the  following  equation. 


Rs  = 


1  -  £1  1  1  -  £2 

+  — + 


£l-dl  -42  £2^2 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


dRs 

ds1 

dR 


-  -F  _1- 

—  tl  , 


dR c 


ds 


°R2  =  l  <*E. 


=  24r2(i-£l-1) 

—  ^2  _ 2  (2  —  £2_1) 


S  -  _ c  -2  a  -1. 

—  ^2  n2  > 


dA9 


a/? 


°RS2  =  (^bX~£1  ^1  X))  +  (CTA1  (f>l  2(1-^1  1))) 

+  (o'£2(-£2_2^2”1))  +  (oh2  (^2_2(2  -  £2_1))) 


Using  the  same  concept  for  using  a  general  tolerance  of  ±0.02”  to  all  physical 
dimensions  the  following  uncertainties,  aAl  and  oAz  were  found  from  the  following 
calculation. 


aAl2  =  (0.00005 (J/1/hs))2+  ^0.00005  j  =  2.67  x  10 
er^22  =  (0. 0000 5 (M^))2 +(0.0000 5(h))2  =  2.68  x  10~9m2 


—9  2 

m 
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Pressure  Drop 


The  pressure  loss  across  the  heat  sink  was  represented  in  the  non-dimensional  form  as 
the  pressure  loss  coefficient,  p.  This  measurement  was  determined  from  the  following 
equation. 

A Pstc 

The  uncertainty  of  this  measurement  was  determined  from  the  following  derivation. 

dp  , 


d^Vst 


=  2  x  p  1  x  Ut 


2  x  CApstc  x  p  2  x  Uh 


=  -4  x  CApstc  x  p  1  x  Uh 


°^VStc  ' 


2  2 
V  =  (crApstc(2  x  P”1  x  +  (op(-2  X  C Apstc  x  p“2  x 

+  (^UhsC- 4  x  CA Pstc  x  p  1  X  t/HS  3)) 

The  three  variables  in  this  derivation  are  Apstc,p  &  UHS.  Each  of  these  have 
propagated  uncertainties  that  will  need  to  be  derived,  with  the  exception  of  A pstc.  A pstc 
is  the  measurement  for  the  static  pressure  difference  across  the  heat  sink.  It  only  needed 
to  be  multiplied  by  a  conversion  factor  (C  =  248.84  Pa/inH20)  to  convert  it  from  inH20 
to  Pa.  The  other  uncertainties,  p  8lUhs  were  then  found. 

The  air  density  (p)  was  calculated  through  the  following  equation. 

_  CPabs 

P  —  'DT 

Jvi  amb 

where,  JZ  is  the  gas  constant  for  air  and  contributes  no  uncertainty.  However,  Pabs  & 
Tamb  are  measurements  and  do  have  added  uncertainty.  The  uncertainty  of  p  was 
determined  from  the  following  derivation. 


dp 

dPabs 

dp 


=  c  xJZ  1  x  Tn 


CPabs  X  JZ  X  Ta 
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on 


2 


°P2  l  GPabs 


(°PaJC  X  X-1  X  Tc 


amb 


dP  A 2  /  dP  A 2 

dPabs'  \Tamb  dTamb) 

_1))  +(ffram6(-CPabsX^-1 


x  T, 


amb 


In  this  derivation,  c  is  the  conversion  from  PSI  to  Pa  (c  =  6894.757  Pa/PSI). 

The  uncertainty  for  Tamb  has  been  previously  computed  during  the  radiation  loss 
uncertainty  derivation.  <JPabs  was  found  from  the  calibration  of  the  absolute  sensor. 

The  next  step  was  to  find  the  uncertainty  of  the  heat  sink  velocity.  The  equation 
for  measuring  the  heat  sink  air  velocity  is  shown  below. 


Uhs  — 


A 


Duct 


I  2  CpDyC 


A 


HS 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


dUfis  _  V2  1 

~YcPoyc 


-1/2 


1-1/2 


A 


Duct 


dPDyc 

9Uhs  _  V2 
~g^~  ~  ~  ~  CPoyc  P 

dUfis 


A 


HS 


At 


l/2n-3/2  “Duct 


d^Duct 

dUHs 
dA 


2  -ruy  c  r-  4 

^2  Cpoyc1^2  P  1^2A-hs 


HS 

-1 


—  V2Cp£)yc1^2p  x!2 ADuct^ A 


HS 


HS 


dUJ 


°uHS  =  | 


HS 


°UHS  =  I 


VDyc 


(f(CPoyc) 


dpD 


yc/ 


VDyc 


(  dUHS ^  /  dUHS  \  /  dUHS \ 


— 1/2  „  — 1/2  m>uct 
AhS 


+ 


aADuct  (V2(CpDyc)2p  zi4W5  J  +  f°hHS  (  V2(CpDyc)2p  2) 


The  remaining  uncertainties  from  this  equation  are  &aDuc[  &  OhHS.  Using  the  same 
approach  as  previously  outlined,  the  uncertainties  for  these  two  measurements  is  shown 
below. 

aADuct2  =  (0.00005 (0.05m))2 +(0.00005(0. 05m))2  =  1.29  x  10"9m2 
GAhs2  =  ((5)  x  0.00005(b))2+((5)  x  0.00005(h))2  =  4.03  x  10-9m2 
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•  Reynolds  Number 

The  last  uncertainty  needing  computation  was  that  of  the  Reynolds  number.  Both 
the  Nusselt  number  and  pressure  loss  coefficient  are  dependent  on  Reynolds  number. 
Because  the  Reynolds  number  was  not  an  independent  variable  its  uncertainty  was  not 
needed  when  computing  the  error  of  the  experiment.  However,  it  was  computed  for 
reference. 


The  equation  for  the  Reynolds  number  within  a  heat  sink  channel  is  as  follows. 

p  x  UHS  x  Dh  hs 
Re„s  = - - - - 

The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 

dReHS 


c,  Uhs  x  Dh  hs  x  P 
dp 

dRej-is 

— —  =  p  x  Dh  hs  x  p 

dUHS 

dReHS 


-l 


°ReHs2  ~  \°P 


dD 

dRejis 
dp 


—  P  x  UHS  x  p  1 


H_HS 
2 


)  +(aa^^)  +(<7d 


dRe 


HS 


HyHS  dD 


H  HS  ) 


JReHS 


—  {op{Uhs  x  Dh_hs  X  M  X))  +  {^Uusip  x  Rh_hs  x  M  X))  +  dh  hs(P  x  ^ hs  x  M  X)) 


Of  the  four  uncertainties  in  this  derivation,  op,  &uHS,  <JDh  hs  &  two  have  been 
found  previously,  while  the  uncertainties  of  the  DH _HS  and  p  still  need  to  be  found.  The 
dynamic  viscosity,  p,  was  determined  from  the  thermo-physical  properties  of  air.  For  this 
reason  it  was  constant  and  considered  to  have  no  uncertainty.  The  equation  for  the 
hydraulic  diameter  of  the  heat  sink  channel  is  shown  below. 

n  _  4  x  AHs_ch 
uh_hs  —  7, 

'wet 


The  uncertainty  of  this  measurement  was  detennined  from  the  following  derivation. 


3  DHhs 
dAHs_ch 


4  x  Pwet 


dD 


dR 


H-HS  =  _4  x  Am_ch  x  Pu  -2 


wet 


wet 


2  2 
°DHHs 2  =  (°h/fS_Ch(4  X  Pwet  X))  +  (^PwetC-4  X  AHs_ch  X  Pwet  2)) 


The  uncertainty  of  the  heat  sink  wetted  area,  AHS,  has  already  been  derived,  but 
needed  to  be  modified  so  that  it  would  apply  to  a  single  channel,  shown  below. 
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(Tahsch 2  =  (0.0000  5  (b))2+(0. 00005(f))2  =  1.61  x  10 -10m2 


The  last  uncertainty  is  for  the  wetted  perimeter  of  the  heat  sink.  This  was  simply 
the  perimeter  of  a  single  heat  sink  channel.  The  equation  for  Pwet  is  shown  below. 


wet 


21  +  2b 


By  using  the  same  principal  of  applying  a  general  tolerance  to  its  dimensions  the 
following  uncertainty  were  derived. 


2  /  dPwet\ 


2 

)  +(°i 


dP 


wet 


V 


db  ) 


aPwet2  =  (0.00005(2))2  +  (0.00005(2))2  =  1.41  x  10“4m 


108 

Approved  for  public  release;  distribution  unlimited. 


•  Results 


Table  H  -  2  Experimental  Results  and  Uncertainty 


Testing  Parameters 

Experimental  Results  and  Uncertainties 

Target 

Chamber 

Pressure 

(atm) 

Target 
Flow  (% 
of  Max) 

Heater 

Set  point 
(W) 

Heat  Sink 
Reynolds 
Number  Re 

Uncertainty 
of  the 
Reynolds 
Number  (±) 

Average 

Nusselt 

Number 

Uncertainty 
of  the 
Nusselt 
Number  (±) 

Pressure 

Loss 

Coefficient 

Uncertainty 
of  the  Loss 
Coefficient 

(±) 

ReD 

±oRe 

Nu 

n 

±on 

0.1 

100% 

20 

1223 

97 

349.8 

21.1 

2.043 

0.025 

0.1 

90% 

20 

1072 

85 

335.7 

21.3 

2.164 

0.031 

0.1 

75% 

20 

872 

69 

302.8 

19.9 

2.424 

0.046 

0.1 

60% 

20 

709 

56 

284.6 

18.9 

2.731 

0.071 

0.1 

45% 

20 

472 

37 

243.6 

18.4 

3.561 

0.150 

0.2 

100% 

20 

2815 

222 

571.1 

29.6 

1.528 

0.009 

0.2 

80% 

20 

2193 

173 

489.7 

25.6 

1.663 

0.014 

0.2 

60% 

20 

1636 

129 

410.8 

22.0 

1.854 

0.025 

0.2 

40% 

20 

1053 

83 

333.3 

19.3 

2.141 

0.060 

0.2 

30% 

20 

759 

60 

294.5 

18.7 

2.712 

0.115 

0.3 

100% 

20 

4226 

334 

734.9 

33.5 

1.282 

0.006 

0.3 

80% 

20 

3349 

264 

598.1 

26.9 

1.406 

0.009 

0.3 

60% 

20 

2490 

197 

506.2 

23.2 

1.540 

0.016 

0.3 

40% 

20 

1603 

127 

396.1 

18.9 

1.808 

0.039 

0.3 

20% 

20 

756 

60 

296.5 

16.9 

3.878 

0.187 

0.4 

90% 

40 

6185 

488 

955.4 

30.4 

1.202 

0.004 

0.4 

80% 

40 

5255 

415 

833.8 

26.0 

1.218 

0.005 

0.4 

60% 

20 

4001 

316 

627.7 

22.9 

1.312 

0.009 

0.4 

40% 

20 

2647 

209 

522.6 

24.6 

1.529 

0.021 

0.4 

20% 

20 

1186 

94 

329.1 

19.4 

1.937 

0.096 

0.5 

80% 

40 

7210 

569 

1063.5 

33.1 

1.179 

0.004 

0.5 

60% 

20 

5465 

431 

808.9 

31.5 

1.191 

0.006 

0.5 

40% 

20 

3619 

286 

645.3 

28.8 

1.331 

0.013 

0.5 

20% 

20 

1742 

138 

419.4 

21.3 

1.719 

0.054 

0.6 

100% 

40 

8651 

683 

1205.0 

40.7 

1.155 

0.003 

0.6 

80% 

40 

6950 

549 

1015.2 

34.9 

1.177 

0.005 

0.6 

60% 

40 

5216 

412 

810.0 

28.8 

1.212 

0.008 

0.6 

40% 

20 

3402 

269 

570.0 

25.7 

1.353 

0.019 

0.6 

20% 

20 

1646 

130 

412.8 

25.8 

1.803 

0.080 

0.7 

100% 

40 

9742 

769 

1315.9 

46.2 

1.139 

0.003 

0.7 

80% 

40 

7805 

616 

1098.1 

38.2 

1.164 

0.004 

0.7 

60% 

40 

5855 

462 

892.4 

31.4 

1.193 

0.008 

109 

Approved  for  public  release;  distribution  unlimited. 


0.7 

40% 

40 

3849 

304 

666.4 

26.8 
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